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ABSTRACT 


Several  reports  in  the  literature  imply  extremely  large  number  densities 
of  metastable  nitrogen  molecules  can  be  found  in  the  effluents  of  dielectric- 
barrier  or  ozonizer  discharges  in  flowing  nitrogen.  Ve  have  constructed  a 
dielectric-barrier  discharge  which  operates  at  1  to  10  kV  in  flowing  nitrogen 
or  mixtures  of  nitrogen  in  either  argon  or  helium  at  total  pressures  between 
0.3  and  6  Torr.  Spectroscopic  observations  in  the  afterglow  of  this  discharge 
indicate  that  virtually  all  observed  excitation  between  200  and  850  nm  results 
from  discharge  streamers.  All  emissions  are  quenched  upon  adding  a  grounded 
loop  down  stream  from  the  discharge.  The  absence  of  NO  Gamma-band  emission 
when  NO  is  added  downstream  from  the  grounded  loop  indicates  that  N2(A)  number 
densities  in  the  afterglow  are  below  10^  molecules  cm~3f  some  ten  orders  of 
magnitude  below  previous  estimates. 


I .  INTRODUCTION 


The  dielectric-barrier  (ozonizer)  discharge  is  a  high  voltage  ac  discharge. 
A  number  of  groups  have  investigated  the  afterglow  of  dielectric-barrier  dis¬ 
charges  through  flowing  nitrogen  for  application  to  spectroscopy3""**  and  trace 
analysis. 5-12  They  have  reported  strong  nitrogen  second-positive  emission 
(N2(C3nu)  -»  N2(B3Ilg))  even  at  times  as  long  as  100  ms  downstream  from  the 
discharge.  Since  the  radiative  lifetime  of  N2(C3nu)  is  37  ns,33  these  results 
suggest  that  N 2 ( C 3 IIU )  is  being  produced  somehow  in  the  discharge  afterglow. 

Ung**  postulated  that  the  strong  second-positive  emission  resulted  from  energy 
pooling  reactions  of  N2(A3EU+)  either  with  itself  or  with  highly  vibrationally 
excited,  ground-electronic  state  nitrogen: 

2N2(A3V)  N2(C3nu)  +  N2(X3  Eg+ )  (1) 

N2(A3V)  +  N2(X3Ig+,v'  >  20)  -»  N2(C3nu)  +  N2(X1Eg+)  .  (2) 

Dodge  and  Allen, 6  noting  the  fall  off  from  linearity  of  atomic  zinc  emissions 
as  a  function  of  added  [ Zn ]  estimated  N2(A)  number  densities  to  be 
1033  molecules  cm-3  or  about  five  to  six  orders  of  magnitude  greater  than  is 
produced  by  microwave3^  or  hollow-cathode33  discharges.  Jurgensen  and 
Winefordner  have  estimated  similarly  large  N2(A)  number  densities, 

1036  molecules  cm  3. 

Intrigued  by  these  reports  of  large  N2(A)  number  densities,  we  undertook  a 
spectroscopic  study  of  the  dielectric-barrier  discharge  afterglow  in  order  to 
determine  number  densities  of  metastable  nitrogen  quantitatively.  Our  pro¬ 
cedure  involves  determining  absolute  photon  emission  rates  of  various 
spectroscopic  features.  In  situ  calibration  of  our  detection  system  relies  on 
observations  of  the  air  afterglow  chemiluminescence  under  carefully  controlled 
conditions.  The  absolute  photon  yields  of  this  reaction  at  visible  wavelengths 
are  well  established.3'7-3^ 
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II.  EXPERIMENTAL 


A.  Apparatus 

The  afterglow  of  the  discharge  was  monitored  as  it  passed  through  a  1  in. 
diameter  glass  flow-tube  with  quartz  windows  at  various  points  along  its 
length  (see  Fig.  1).  A  wire  loop  was  inserted  through  a  1/4  in.  port  at  the 
entrance  to  the  flow-tube.  The  loop  was  grounded  through  a  1  kfi  resistor  and 
served  as  a  groundpoint  to  ensure  the  discharge  did  not  extend  into  the  obser¬ 
vation  region.  A  Sargent-Welch  No.  1402  (160  1pm)  mechanical  pump  exhausted 
the  flow  tube,  maintaining  flow  velocities  on  the  order  to  300  cm/s~l. 

Nitrogen  flow  rates  ranged  from  20  to  200  umoles-^  and  total  pressures 
varied  between  0.5  and  6  Torr.  Most  spectra  were  taken  28  cm  (~  100  ms) 
downstream  from  the  discharge  terminus.  These  conditions  were  picked  to 
reproduce  those  of  Dodge  and  Allen^  and  of  Jurgensen  and  Winefordner . ^ 

The  discharge,  which  is  based  on  Dodge  and  Allen's  design, 6  consists  of 
two  coaxial  glass  tubes  which  form  a  2.8  mm  annulus  (see  Figure  2).  One 
electrode  is  formed  by  pinching  off  the  bottom  of  the  inner  tube,  filling  it 
with  saturated  NaCl  solution  and  immersing  a  wire  from  the  power  supply  in  the 
solution.  The  second  electrode  is  made  by  wrapping  copper  foil  tape  around  the 
outer  tube.  The  electrodes  form  a  disci arge  region  that  is  30  cm  long. 

The  ac  power  supply  used  with  the  discharge  consisted  of  a  15  kV  neon  sign 
transformer  (center  tap  of  secondary  grounded)  with  a  variac  in  the  primary 
circuit  to  control  the  output  voltage.  To  reduce  interference  with  nearby 
electronics,  only  half  of  the  transformer  secondary  was  used  during  most  of  our 
work.  Two  100  k 52  power  resistors,  placed  in  the  secondary,  limited  the 
current . 

Spectra  were  taken  using  a  Spex  No.  1269,  1.26m.  monochromator  (1  to  3  mm 
slits)  and  a  Hamamatsu  No.  R955  photomultiplier  tube.  The  photomultiplier  tube 
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was  cooled  to  -40°C  and  operated  at  a  gain  of  1.7  x  10^.  The  output  of  the 
photomultiplier  tube  was  measured  using  a  Keithley  No.  417  picoameter  which  fed 
into  a  Soltec  No.  1242  strip  chart  recorder. 

B.  Photometric  Calibrations 


Calibrating  the  apparatus  for  absolute  photon  emission  rates  was  done  in 
two  steps.  First,  the  relative  response  function  (R\)  of  the  monochromator  was 
measured  using  an  ultraviolet  irradiance  standard  between  200  and  400  nm  and 
the  relative  intensities  of  the  continuum  produced  by  the  0/N0  air  afterglow 
reaction  between  400  and  800  nm  (react: on  (3)).^~20  This  relative  response 
function  corrects  for  wavelength  variations  in  detector  efficiency  and  the 
transmission  of  the  optical  system.  Second,  the  absolute  response  at  580  nm 
was  measured  using  the  0/N0  air  afterglow  titration. 

When  oxygen  atoms  and  nitric  oxide  are  mixed,  the  following  reaction  acts 
to  produce  a  continuum  that  extends  from  3/5  nm  to  beyond  3000  nm: 

0  +  NO  -»  N02*  +  M  (3a) 

N02*  ->  N02  +  h\>  (a  continuum)  .  (3b) 

The  intensity  of  the  emission  is 

i°/N0  =  kX(0,lN0,*X  ’  (4) 

where  k^  is  the  air  afterglow  rate  coefficient  (k\  =  1.25  x  10"^  photon 
cm"-*  s"1  nm~^  at  580  nm)^  and  AX  is  the  monochromator  bandwidth.  Measuring 
the  air  afterglow  intensity  at  a  fixed  wavelength  under  conditions  of  known  [Oj 
and  [N0J  place  the  relative  response  function  on  an  absolute  basis. 
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In  the  titration  known  number  densities  of  0  atoms  are  formed  by  mixing  NO 
with  a  constant  flow  of  N  atoms.  The  N  atoms  are  produced  by  passing  nitrogen 
gas  through  a  microwave  discharge.  The  NO  quantitatively  converts  the  N  to  0 
through  reaction  (5), 21 

N  -■  NO  -*  0  +  N2  (5) 

Veil  below  the  titration  point  N2  first-positive  emission  is  observed  from 
the  recombination  of  N  atoms: 


N  +  N  +  M  ->  N2(B3ng)  +  M  ,  (6) 

N 2 ( B 3 Ilg )  -*  N2(A3£u  +  )  +  hv  (first  positive)  .  (7) 


As  more  NO  is  added  the  intensity  of  the  first-positive  emission  drops.  At  the 
titration  point  no  emission  from  the  flow  is  observed  and  the  0-atom  number 
density  is  equal  to  the  initial  N-atom  number  density.  As  still  more  NO  is 
added,  the  0/N0  air  afterglow  is  observed  and  its  intensity  varies  linearly 
with  the  number  density  of  NO  added: 

I°/N°  =  K[N]q  ( [ NO ] — [ N I o )  ,  (8) 


where  subscript  0  refers  to  N-atom  number  density  determined  from  the  titration 
and  point.  Equation  (8)  shows  that  K  will  be  determined  by  the  ratio  of  the 
square  of  the  slope  to  the  intercept  of  the  line  determined  by  the  variation  in 
intensity  of  the  air  afterglow  as  a  function  of  added  NO  number  density. 

Combining  Eqs.  (h)  and  (8)  converts  the  observed  intensities  (in  nA)  to 
absolute  photon  emission  rates: 

It  rue  *  xob~  k580  AX  R580/(K  RXobs>  ■ 
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Where  I0bs  *s  the  total  integrated  band  intensity  of  the  spectral  feature  of 
interes  t . 

The  air  afterglow  intensities  must  be  corrected  for  the  slow  removal  of  0 
atoms  in  the  following  reactions: 

0  +  NO  +  M  -*  NO2  +  M  (k  =  7  x  10--^^  cm^  molecules_^s ) ,  22  (10) 

and 

NO2  +  0  -»  NO  +  O2  (k  =  9.5  x  10“12  cm3  molecule--*-  s-*-,^^  (11) 

where  k  is  the  rate  coefficient  of  the  corresponding  reaction. 
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III.  RESULTS 


A.  General  Observations 


Figures  3  and  4  show  typical  spectra  in  the  afterglow  region  with  the  loop 
floating  and  with  it  grounded.  The  second-positive  bands  (N2(C^nu)  ->N2(B3[ig)) 
were  by  far  the  strongest  system.  Other  observed  systems  included  the 
first-positive  system  (^(B-^Ilg)  ->  N2<  A^  Eu  + ) ,  N2+  first-negative  bands  at  391  and 
428  nm  (v',v"  =  0,  0  and  0,  1  of  N2  +  (B^EU+)  -»  N2  +  (X^EU+)),  and  a  relatively  weak. 
NOy  system  (N0(A^Z+)  -*  NO(X^II)).  The  presence  of  the  NOy  bands  -  caused  mostly 
*  by  trace  oxygen  in  the  nitrogen  -  made  it  impossible  to  detect  the  Vegard-Kaplan 

bands  that  would  have  signaled  the  presence  of  N2(A)  in  the  discharge  effluent. 
The  output  of  the  discharge  was  observed  for  pressures  ranging  from  0.3  to 
6  torr  and  supply  voltages  from  1  to  10  kV.  The  observed  intensities  were 
roughly  proportional  to  the  voltage  but  pressure  seemed  to  have  the  most  effect 
on  the  discharge  output.  At  higher  pressures  the  output  was  a  purple-pink 
color.  The  flow  was  also  very  unsteady  in  appearance.  At  0.3  torr  the 
discharge  output  was  bright  pink  and  uniform.  At  lower  pressures  the 
intensities  of  the  first-  and  second-positive  systems  were  greater  than  at  high 
pressures  and  a  few  bands  in  each  system  emitted  nearly  all  of  the  signal. 

With  the  loop  grounded  the  afterglow  was  non-existent  downstream  from  the 
groundpoint.  Spectra  taken  under  these  conditions  showed  nothing  except 
scattered  light  -  as  was  verified  when  the  flow  tube  was  modified  to  position 
the  ground  point  upstream  from  the  light  trap. 

Measurements  of  the  current  flowing  out  of  the  groundpoint  were  made  by 
monitoring  the  voltage  drop  across  a  1  kS?  resistor  between  the  groundpoint  and 
ground.  At  6  torr  the  current  was  about  60  mA  (with  5400V  across  the 
discharge).  At  lower  pressures  the  figure  rose  to  around  400  mA. 
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B.  Estimate  of  N?(A)  Number  Density 

We  looked  for  the  excitation  of  N0(A)  by  N2(A)  to  determine  if  N2(A)  were 
present  in  the  afterglow: 

N2(A)  +  NO  ->NO(A2E+)  +  N2(X1Eg+)  ,  (12) 

NO(A2E+)  ->NO(X2II)  +  hv  (y  bands)  .  (13) 

The  differential  equation  governing  the  number  density  of  NO(A2E+)  is 

~  [m^-  =  kex  [N2(A)  ]  [NO]  -  krad  [N0(A)  ]  .  (14) 

kex  (2  x  10_H  cm^  molecule-^  s~3)23  the  excitation  rate  coefficient  and 
krad  is  the  radiative  decay  rate.2^  N0(A)  is  in  steady  state  in  our  reactor 
because  the  radiative  lifetime  of  N0(A)  is  short  (-  200  ns)  compared  to  the 
time  the  gas  resides  within  the  monochromator  field  of  view  (-1  ms).  This 
means  that 


so 


lN0(A)  =  krad  [N0(A) ]  =  kex  [N2(A)][N0] 


IN2(A)) 


dIN0( A) 
d  [  NO  ] 


(15) 


(16) 


We  looked  for  emission  from  the  0,1  y-band  at  235  nm  when  we  added  NO  to 
the  reactor.  We  could  detect  no  changes  in  emission  intensity  as  the  NO  number 
density  was  increased  to  10^  molecules  cm-3.  We  estimate  a  minimum  detectable 
signal  increase  of  2  x  10~3  nA  which  when  converted  for  absolute  system 
response  corresponds  to  an  absolute  photon  emission  rate  of  1.6  x  10-  photons 
cm'3  s~3.  Using  this  value  in  Eq.  (16)  along  with  the  NO  number  density  and 
kex  implies  an  upper  limit  to  the  N2(A)  number  density  of  2.3  x  UP  molecules 
cm-2.  This  number  is  some  ten  orders  of  magnitude  below  previous  estimates 
which  were  based  upon  more  indirect  evidence. 
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IV.  DISCUSSION 


The  facts  of  the  discharge  afterglow's  terminating  at  the  groundpoint,  the 
relatively  large  electron  currents  observed  at  the  ground  point,  and  the  lack 
of  measurable  N2<A)  in  the  afterglow,  suggest  that  the  afterglow  is  actually  an 
extension  of  the  discharge.  The  N2(C^nu)  observed  by  others  and  by  us  in  the 
afterglow  region  when  the  loop  is  ungrounded  is  continuously  formed  by  electron 
impact  excitation  of  N2  by  discharge  streamers  rather  than  by  energy-pooling 
reactions  such  as  (1)  and  (2).  This  fact  does  not  negate  the  usefulness  of  the 
dielectric-barrier  discharge  as  a  tool  for  exciting  analytically  useful 
spectra.  Useful  spectra  clearly  are  excited.  The  mechanism  of  this 
excitation,  however,  is  different  from  what  has  been  proposed. 

The  lack  of  N2(A)  under  the  conditions  employed  actually  is  not 
particularly  surprising.  N2<A)  is  well  known  to  be  quenched  with  unit 
efficiency  in  collisions  with  reactor  walls. ^  In  order  to  transport  N2(A)  for 
any  appreciable  distance,  therefore,  one  must  work  either  at  quite  high 
pressure  and,  thereby  slow  diffusion  to  the  reactor  walls,  or  else  one's 
observations  must  be  made  only  a  few  milliseconds  downstream  from  the 
discharge. 

Noxon^  was  successful  in  detecting  N2(A)  in  the  afterglow  of  a 
dielectric-barrier  discharge.  He  worked  at  pressures  on  the  order  of  an 
atmosphere  and  could  observe  Vegard-Kaplan  emission  in  the  afterglow  for  times 
on  the  order  of  a  second.  Rice  et  al.'s^  configuration  for  exciting  analyte 
spectra  in  the  afterglow  of  a  dielectric-barrier  discharge  also  employs 
pressures  on  the  order  of  an  atmosphere  and,  in  addition,  quite  high  nitrogen 
flow  rates.  Although  they  do  not  report  having  observed  Vegard-Kaplan 
emission,  their  conditions  should  be  favorable  for  the  presence  of  N2(A)  in  the 
afterglow.  One  sign  that  N2(A)  is  likely  to  be  in  their  afterglow  is  their 
observation  of  atomic  oxygen  emission  at  557.7  nm.  This  emission  is  excited  in 
the  energy-transfer  reaction  between  N2(A)  and  o  atoms. 25 
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Ve  moved  our  dielectric-barrier  discharge  to  a  2  in.  diameter  flow  reactor 
which  has  significantly  enhanced  pumping  capabilities^  to  see  if  N2 ( A )  would 
indeed  persist  in  the  discharge  afterglow.  Ve  observed  Vegard-Kaplan  emission 
quite  readily  upon  discharging  mixtures  of  nitrogen  in  argon  at  pressures  on 
the  order  of  5  to  10  Torr  and  times  on  the  order  of  10  ms  downstream  from  the 
grounded  loop.  Absolute  photometric  observations  of  Vegard-Kaplan  emission 
established  the  N2(A)  number  densities  in  the  observation  region  to  be  about 
1010  molecules  cm-^.  This  number  density  is  comparable  to  what  is  produced  by 
other  laboratory  sources.  The  dielectric-barrier  discharges  therefore,  does 
not  seem  to  offer  significant  advantages  over  other  discharge-flow  methods  of 
N2<A)  generation.  Its  one  advantage  appears  to  be  the  ability  to  generate 
N2<A)  at  pressures  on  the  order  of  an  atmosphere. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  dielectric-barrier  discharge,  as  commonly  configured  for  analytical 
applications,  is  not  a  source  of  N2<A)  for  analytically  useful  spectra. 

Rather,  the  spectra  are  excited  by  direct  electron  impact  by  discharge 
streamers.  When  configured  for  atmospheric  pressure  operation  with  short 
transit  times  between  the  discharge  and  the  observation  region,  N2(A) 
energy-transfer  reactions  may  indeed  provide  the  primary  excitation  source. 
This  high  pressure  configuration  might  be  preferred  if  the  discharge  source  is 
to  be  used  is  a  gas  chromatograph  detector.  The  dielectric-barrier  discharge 
provides  no  advantages  over  other  common  laboratory  sources  when  used  to 
generate  N2(A)  to  study  energy- transfer  reactions  unless  one  wants  to  study 
these  reactions  at  high  pressures. 
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FIGURE  CAPTIONS 


1.  Schematic  diagram  of  apparatus  used  to  study  dielectric-barrier  discharge 
afterglows. 

2.  Details  of  dielectric-barrier  discharge. 

3.  Spectrum  of  afterglow  between  300  and  450  nm  with  the  loop  flotaing  and 
with  it  grounded.  The  spectrum  has  not  been  corrected  for  variations  in 
monochromator  response  with  wavelength. 

4.  Spectrum  of  afterglow  between  550  and  800  nm  with  the  loop  floating  and 
with  it  grounded.  The  spectrum  has  not  been  corrected  for  variations 
monochromator  response  with  wavelength. 
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ABSTRACT 


We  have  extended  a  technique  for  studying  the  vibrational  distributions 
of  ground-electronic-state,  molecular  nitrogen  in  the  afterglow  of  a  microwave 
discharge  through  mixtures  of  helium  and  nitrogen.  The  technique  is  based 
upon  adding  metastable  helium  atoms  to  the  afterglow.  The  He*  ( 2 S )  excites 
the  N2(X,v)  to  N2+(B, ^Tu+)  in  a  Penning-ionization  reaction.  Since  Penning 
ionization  is  a  Franck-Condon  process,  the  vibrational  distribution  of  the 
N2+(B)  product  is  determined  by  that  of  the  N2(X,v)  from  which  it  was  pro¬ 
duced.  The  measurements  show  that  the  ground-state  nitrogen  distribution  is 
highly  non-Boltzmann,  with  vibrationally  hotter  distributions  being  produced 
with  lower  mole  fractions  of  nitrogen  in  the  discharge.  We  have  also  observed 
the  production  of  N2+(C^EU+)  from  He*  Penning  ionization  of  molecular  nitro¬ 
gen.  This  process  is  energetically  allowed  only  if  the  vibrational  energy  in 
the  ground-electronic-state  nitrogen  exceeds  3.8  eV  or  15  vibrational  quanta. 
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I.  Introduction 


Our  recent  observations  in  discharged  nitrogen  afterglows  showed  that 
IF(B3n0+)1  and  N2(B3Hg)3»2  excitation  appeared  to  result  from  V-E  transfer 
rather  than  the  more  conventional  E-E  transfer  from  one  electronically  excited 
species  to  another.  To  investigate  these  processes  more  fully,  we  needed  to 
develop  a  diagnostic  for  N2<v).  This  report  details  our  development  of  an 
N2<X,v)  diagnostic  based  upon  Penning  ionization  by  metastable  helium  atoms. 

Vibration-to-electronic  transfer  processes,  if  efficient,  could  have 
important  implications  for  chemical  laser  development  as  well  as  for  the 
modeling  of  disturbed  atmospheres  because,  unlike  the  lowest  electronically 
excited  nitrogen  metastable,  N2(A,3EU+)3,  ^(v)  does  not  destroy  itself 
efficiently  in  collisions.  Indeed,  quite  high  levels  of  vibrational  excita¬ 
tion  can  be  achieved  by  V-V  up-pumping^-^ .  Thus  N2(v)  might  act  as  an  energy 
storage  medium  in  a  laser  system,  and  N2(v)  created  in  an  atmospheric  distur¬ 
bance  could  cool  only  by  collisions  with  unlike  species.  We  note  that  Starr 
observed  V-E  excitation  of  Na  and  K  by  N2(v)®>^,  and  the  work  of  Rich  and 
co-workers  provides  evidence  of  efficient  V-E  transfer  in  systems  involving 
CO(v)10  and  NO(v)11. 

A  number  of  previous  studies  of  the  effluents  of  active-nitrogen  dis¬ 
charge  afterglows  have  shown  qualitatively  the  presence  of  vibrationally 
excited  nitrogen.  Quantitative  assessment  of  the  degree  of  ground-state 
vibrational  excitation,  however,  generally  has  proved  difficult.  Because  N2 
is  not  an  infrared  active  molecule,  the  degree  of  vibrational  excitation  in 
the  active  nitrogen  has  been  probed  primarily  by  indirect  methods.  Kaufman 
and  Kelso3^  remarked  on  the  heated  walls  of  their  flow  reactor  downstream  from 
their  discharge  and  attributed  the  source  of  this  heat  to  deactivation  of 
vibrationally  excited  nitrogen  molecules  in  collisions  with  the  walls. 

Morgan,  Phillips,  and  Schiff33*^  noted  an  excess  temperature  given  to  the 
catalytic  probe  placed  in  their  flow  reactor  which  could  not  be  attributed 
N-atom  recombination.  Again  they  suggested  that  the  source  of  this  excess 
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heat  was  vibrationally  excited  nitrogen  molecules.  Starr^  and  Starr  and  Shaw^ 
observed  excitation  of  sodium  and  potassium  electronic  transitions  when  they 
introduced  sodium  or  potassium  atoms  to  the  afterglow  of  a  nitrogen  discharge. 
They  suggested  that  the  excitation  resulted  from  a  vibrational-to-electronic 
(V-E)  excitation  from  N£(v)  to  Na  or  K.  Bass^  an(j  Tanaka  et  al.*6  used 
vacuum  ultraviolet  (VUV)  absorption  spectroscopy,  employing  the  partially 
forbidden  Lyman-Birge-Hopf ield  bands,  to  show  that  N2(v)  was  indeed  present  in 
nitrogen-discharge  afterglows.  A  number  of  other  diagnostics  have  been  devel¬ 
oped  in  more  recent  years  including  photo-ionization^ ,  photoelectron  spectros- 
copy!8»19,  electron  energy  loss2^,  Penning  ionization^l-^^,  Raman^-^, 
CARS28»29?  an(j  Multiphoton  ionization-^'-^,  among  others.  Caledonia  et  al. 
have  reviewed  much  of  this  work  in  a  recent  report^. 

Our  diagnostic  for  vibrationally  excited  nitrogen  in  discharge  afterglows 
extends  the  pioneering  work  of  Schmeltekopf  et  al.21  and  Young  and  Horn22>23> 
They  relied  on  the  vertical  nature  of  Penning-ionization  transitions  from 
ground-electronic-state,  neutral  nitrogen  to  produce  N2+(B2Eu+,v)  distribu¬ 
tions  characteristic  of  the  ground-state  distributions.  Mixing  metastable 
helium  atoms  with  a  flow  of  molecular  nitrogen  results  in  strong  emission  of 
the  nitrogen  first-negative  system,  N2+(B  2EU+  -X2Eg+).  Since  the  Penning- 
ionization  process  follows  a  Franck-Condon  excitation  pathway,  the  vibrational 
distribution  in  the  neutral,  ground  state  will  determine  the  distribution 
observed  in  the  upper,  ionic  state.  One  problem  with  this  approach  is  that 
care  must  be  taken  not  to  have  any  He+  or  He2+  in  the  flow  of  metastable 
helium.  Both  of  those  species  also  excite  N2+(B)  quite  strongly  in  charge- 
transfer  reactions,  but  with  an  N2+(B)  vibrational  distribution  that  is 
decidedly  non-Franck-Condon33.  This  diagnostic  is  most  sensitive  to  vibra¬ 
tional  excitation  of  the  first  six  vibrational  levels  in  the  ground  electronic 
state.  Our  goal  is  to  be  able  to  estimate  the  overall  vibrational  distribu¬ 
tion  of  the  ground-state  nitrogen  by  combining  our  experimental  determination 
of  the  vibrational  distribution  in  the  lower  levels  with  modeling  calculations 
on  the  temporal  development  of  the  fully  coupled  vibrational-state  manifold^’3. 
Jolly  et  al.2^  have  also  re-investigated  the  Penning-ionization  technique. 
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II .  Experimental 


The  discharge-flow  apparatus  used  in  these  experiments  has  been  described 
in  most  respects  previously34-39.  Ve  summarize  briefly  the  modifications 
required  for  the  measurements  related  to  N2(v).  Figure  1  shows  the  configura¬ 
tion  of  the  flow-  reactor.  A  microwave  discharge  at  the  upstream  end  of  the 
flow  reactor,  through  a  flow  of  helium  and  nitrogen,  creates  the  active  nitro¬ 
gen.  It  enters  a  section  of  2  in. -diameter  Pyrex  ®  containing  three  side 
arms.  Small  amounts  of  SFg  flow  through  the  first  side  arm  to  attach  elec¬ 
trons  created  in  the  Penning-ionization  reactions  and  the  metastable  helium 
atoms  enter  through  the  third  side  arm  to  produce  the  Penning-ionization 
spectrum.  A  0.5m  monochromator,  mounted  on  rails,  views  the  region  of  the 
flow  reactor  just  downstream  from  the  He*  injector  to  observe  the  fluorescence 
created  in  the  Penning-ionization  reaction. 

A  hollow-cathode,  d.c.  discharge  through  a  flow  of  purified  helium 
creates  the  metastable  helium  atoms.  Flowing  the  helium  through  a  molecular- 
sieve  trap  at  liquid  nitrogen  temperature,  upstream  from  the  discharge, 
removes  most  of  the  impurities  in  the  helium,  including  nitrogen  and  oxygen. 
Small  amounts  of  residual  neon  in  the  helium  ( <10  ppm)  have  no  effect  on  our 
observations.  The  discharge  operates  at  350V  with  a  current,  limited  by  a 
200  kQ  resistor,  below  a  milliamp. 

We  showed  previously-^  that  operating  the  hollow-cathode  discharge  at 
high  voltage  and  high  current  tended  to  produce  significant  number  densities 
of  atomic  and  molecular  helium  ions.  These  ions  are  an  anathema  to  the  Penning- 
ionization  measurements  because  they  produce  highly  non-Franck-Condon  distribu¬ 
tions  in  the  N2+(B).  Thus  one  could  be  led  to  false  conclusions  regarding  the 
extent  of  ground-state  vibrational  excitation.  The  absence  of  significant  num¬ 
ber  densities  of  helium  ions  was  confirmed  by  failing  to  observe  N2+(B,  v'  >  2) 
when  the  active-nitrogen  discharge  was  off. 
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Typical  conditions  comprised  flow  rates  of  helium  through  the  microwave 
discharge  and  the  metastable-producing  d.c.  discharge  of  3600  and  1200  pmol  s-^ 
respectively,  nitrogen  flow  rates  between  25  and  250  pmol  s-*,  SFg  flow  rate 
of  4  pmol  s-1 ,  total  pressure  of  1.5  torr,  and  transit  time  between  the  micro- 
wave  discharge  and  observation  region  of  ~11  ms.  Generally,  the  microwave 
power,  which  was  coupled  through  a  McCarroll  cavity,  was  75  watts. 

A  least-squares  fitting  procedure^  analyzes  the  Penning-ionization 
spectra.  We  generate  basis  functions  consisting  of  a  synthetic  electronic 
spectrum  for  a  unit  population  in  each  vibrational  level  of  each  electronic 
state  appearing  in  the  spectral  region  of  interest.  A  linear  least-squares 
routine  then  finds  the  populations  of  each  vibronic  band  which,  when  multi¬ 
plied  by  the  appropriate  basis  function  and  summed  with  overlapping  bands, 
gives  a  composite  spectrum  most  nearly  matching  the  experimental  spectrum. 

III.  Theory  Behind  the  Penning-ionization  Measurements 

The  Penning  ionization  between  metastable  helium  atoms  and  molecular 
nitrogen  is  a  vertical  process.  One  can  calculate  the  vibrational  distribu¬ 
tion  in  the  final  state,  therefore,  knowing  only  the  vibrational  distribution 
in  the  lower  state  and  the  Franck-Condon  factors  that  couple  the  two  states 
together.  Thus 


N 


«  E 


Nyll  Qy/yll 
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where  v'  and  v"  represent  the  vibrational  levels  of  the  upper  and  lower 
states,  respectively,  and  qvrvii  is  the  Franck-Condon  factor  coupling  them. 

We  calculated  .a  set  of  Franck-Condon  factors  over  the  range  of  ground- 
electronic-state  vibrational  levels  of  0-18,  and  N2+(B)-state  vibrational 
levels  of  0-9  using  a  calculat ional  procedure  previously  described. Table  1 
contains  the  results  of  these  calculations. 
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In  theory,  if  one  measures  the  vibrational  distribution  in  the  upper 
state,  that  in  the  ground  state  can  be  determined  simply  by  inverting  the 
Franck-Condon  matrix,  and  multiplying  this  inverse  matrix  on  both  sides  of 
Eq.  (1).  This  procedure  did  not  work  for  us  because  the  measurements  of  the 
upper-state  populations  have  some  uncertainty  associated  with  them,  and  the 
uncertainties  become  greatly  magnified  by  the  matrix  multiplication. 

Jolly  et  al.24  apparently  succeeded  in  using  this  approach  to  analyze  ^(v) 
distributions  created  in  a  glow  discharge.  Ve  do  not  understand  their  success 
in  light  of  our  lack  of  it.  Our  spectral  fitting  approach  for  determining 
N2+(B)  number  densities  ought  to  be  more  accurate  than  their  graphical  inte¬ 
gration  of  one  branch  from  each  band.  Ve  analyzed  our  data,  therefore,  using 
a  model  to  describe  the  N2<X)  vibrational  distributions  which  we  then  used  to 
predict  N2+(B)  vibrational  distributions.  We  then  compared  these  predictions 
with  observation. 

Using  the  Franck-Condon  factors  in  Table  1,  we  computed  the  N2+(B)-state 
vibrational  populations  relative  to  that  in  N2+(B,  v'  =  0)  expected  for  a 
number  of  values  of  vibrational  temperature  of  ground-state  nitrogen  between 
300  and  30,000  K,  assuming  a  Boltzmann  distribution  among  the  levels. 

Figure  2  displays  the  results  of  these  calculations.  If  the  ground-state 
vibrational  level  populations  follow  a  Boltzmann  distribution,  one  can  use 
experimentally  derived  N2+(B)  population  ratios  and  Figure  2  to  find  the 
vibrational  temperatures  of  ground-state  nitrogen  which  corresponds  to  those 
excited-state  population  ratios. 

If  the  ground-electronic-state  vibrational  levels  do  not  follow  a 
Boltzmann  distribution,  one  can  still  compute  expected  excited-state 
populations  from  Eq.  (1),  if  the  ground-state  distribution  is  known. 

Effluents  from  nitrogen  discharges  generally  have  nonequilibrium  vibrational 
distributions,  with  the  higher  vibrational  levels  more  strongly  populated  than 
would  be  predicted  on  the  basis  of  a  Boltzmann  distribution.  Caledonia  and 
Center^  and  Capitelli  and  Dilonardo^  have  shown  that  for  low  vibrational 
levels  the  analytical  distribution  given  by  Treanor  et  al.^  holds: 


N 


v"=0 


=  exp 


1 .4388  -  2w_ x 


e  e, 


e. 


(v"-l) 


1.4388  w  x 
e  e 


(2) 


where  9^  is  the  Boltzmann  vibrational  temperature  referenced  to  v"  =  1,  T  is 
the  ambient  gas  temperature,  and  we  and  ooexe  are  spectroscopic  constants  in 
units  of  cm~l.  The  Boltzmann  vibrational  temperature  is  given  by 
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This  distribution  goes  through  a  minimum,  generally  referred  to  as  the  Treanor 
minimum,  at  vibrational  level  v*f  given  by 
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For  vibrational  levels  greater  than  the  Treanor  minimum,  the  product  v"Nv"  is 
essentially  constant.  The  resulting  distribution  for  v"  >  v*,  therefore,  is 
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where  the  various  parameters  in  Eq.  5  are  determined  by  requiring  that  the  two 
distributions  be  equal  at  v"  =  v*-l. 

The  Penning-ionization  spectrum  calculated  from  a  non-Boltzmann,  ground- 
state  distribution  based  upon  Eqs .  (2),  (4)  and  (5)  is  significantly  different 
from  what  one  would  calculate  from  a  Boltzmann  ground-state  distribution.  The 
Penning-ionization  technique,  therefore  ought  to  differentiate  between  the  two 
ground-state  distributions  rather  easily. 
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IV.  Results  of  the  Penning-Ionization  Measurements 


Figures  3  and  4  show  a  portion  of  the  nitrogen  first-negative  spectrum 
with  the  active-nitrogen  discharge  off  and  on,  respectively.  Clearly,  the 
vibrational  development  of  the  emission  is  greatly  enhanced  by  discharging  the 
gas.  We  determined  vibrational  populations  in  the  upper  state  by  fitting  the 
spectrum  in  the  manner  described  previously-^ > 38 ,  While  the  fits  included 
upper-state  vibrational  levels  0-8,  for  the  most  part,  only  vibrational  levels 
0-7  contributed  any  significant  intensity  to  the  spectra. 

A  number  of  different  and  conflicting  sets  of  Einstein  coefficients  for 
the  N2+(b2eu+  -  X^Eg+)  fill  the  literature^-^.  In  particular,  the  constancy 
of  the  electronic  transition  moment  with  r-centroid  is  unsettled^^S.  To 
analyze  our  data,  we  calculated  a  set  of  Einstein  coefficients  based  upon  the 
relative  electronic  transition  moment  variation  given  by  Brcwn  and  Landshoff^, 
which  has  received  some  support  by  Comes  and  Speier^  and  by  Chang  et  al.^®, 
the  Franck-Condon  factors  of  Albritton  et  al.^9,  an(j  a  radiative  lifetime  of 
N2+(B,v'=0)  of  62.5  ns^1.  They  are  tabulated  elsewhere. 24  We  are  currently 
investigating  this  issue  of  transition-moment  variation.  The  molecular 
potential  constants  given  in  Lofthus  and  Krupenie^*  were  found  to  be  inade¬ 
quate  for  predicting  band  positions  accurately.  Our  fitting  program  uses  the 
recent  potential  constants  determined  by  Gottscho  et  al.^0 

Early  in  this  measurement  program,  we  discovered  that  moving  the 
monochromator  slightly  downstream  from  the  injector,  through  which  the  He* 
entered  the  reactor,  resulted  in  a  much  different  vibrational  distribution  in 
the  upper  state.  There  appeared  to  be  no  reason  why  the  vibrational  distribu¬ 
tion  in  the  active  nitrogen  should  change  over  short  distances.  A  diffuse 
emission  could  be  seen  extending  somewhat  downstream  from  the  well-defined 
flame  created  by  the  Penning-ionization.  Several  centimeters  downstream  from 
the  injector,  the  Penning-ionization  flame  disappeared,  and  only  the  diffuse 
emission  remained.  Figure  5  shows  the  diffuse-emission  spectrum. 
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On  the  assumption  that  the  emission  was  caused  by  energetic  electrons 
exciting  N2+(X),  a  trace  of  SFg  was  added  to  act  as  a  scavenger  for  the 
electrons.  This  addition  eliminated  the  diffuse  flame,  and  also  resulted  in 
N2+(B)  vibrational  distributions  which  did  not  change  with  the  location  of  the 
monochromator  relative  to  the  He*  injector.  Presumably,  the  electrons  created 
in  the  Penning-ionization  reaction  pick  up  some  extra  kinetic  energy  either 
from  stray  microwave  fields  which  have  penetrated  downstream,  or  else  from 
collisions  with  energetic  species  in  the  active  nitrogen.  The  amount  of 
energy  must  be  fairly  considerable,  because  exciting  N2+(B)  from  N2+(X) 
requires  more  than  3  eV.  All  the  results  given  below  were  obtained  in  the 
presence  of  small  traces  of  SFg. 

Analyzing  the  N2+(B,v')  vibrational  distributions  with  the  help  of 
Figure  2  gave  different  ground-state  vibrational  temperatures  depending  upon 
which  upper-state  vibrational  level  was  considered.  The  data  fell  into  three 
groups.  The  temperature  determined  from  v'  =  1  was  only  63  percent  of  that 
determined  from  v'  =  2  and  3,  which  were  themselves  reasonably  consistent. 
Vibrational  temperatures  determined  from  v'  =  4-7  also  were  reasonably  consis¬ 
tent  with  each  other,  but  were  30  percent  larger  than  temperatures  determined 
from  vibrational  levels  2  and  3.  Clearly,  a  Boltzmann  model  for  the  ground- 
state  vibrational  distribution  is  inadequate. 

Fitting  our  distributions  to  the  analytical  model  described  by  Eqs.  (2), 
(4),  and  (5)  was  more  successful.  Figures  6  and  7  show  two  examples.  Clearly, 
the  nonequilibrium  model  does  a  reasonable  job  of  fitting  observations  out 
through  v'  =  6.  We  do  not  understand  the  sudden  discrepancy  at  v'  =  7.  The 
dashed  lines  in  Figures  6  and  7  show  the  distribution  which  would  be  predicted 
were  the  ground-state  to  be  determined  by  a  Boltzmann  distribution.  This 
shows  graphically  the  inadequacy  of  that  model  for  the  ground-state  levels. 

In  the  future  we  plan  to  use  an  improved  model  to  predict  the  population 
distribution  up  to  higher  vibrational  levels. 
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Our  observations  show  that  effective  vibrational  temperatures,  ,  as 
determined  from  fits  of  the  data  to  Eqs.  (2)  and  (5)  rather  than  what  one 
would  calculate  from  Eq.  (3),  tended  to  be  smaller  with  larger  nitrogen  mole 
fractions  flowing  through  the  discharge  (Figure  8)  and  also  at  higher  total 
pressures.  Placing  a  nickel  screen  in  the  reactor,  downstream  ftom  the  dis¬ 
charge,  only  slightly  reduced  the  effective  vibrational  temperature  of  the 
nitrogen  (about  5  percent).  The  nickel  screen,  however,  reduced  atomic 
nitrogen  number  densities  by  more  than  an  order  of  magnitude.  Placing  a 
screen  in  the  afttrglow,  therefore,  is  a  means  to  produce  a  somewhat  cleaner 
flow  of  N2(v)  in  a  discharge-flow  system.  The  Penning-ionization  technique 
appears  to  provide  a  reasonably  accurate  monitor  of  the  vibrational  distribu¬ 
tion  of  ground-electronic-state  nitrogen,  at  least  for  the  lower  vibrational 
levels.  Extending  the  model  to  include  higher  vibrational  levels  would  need 
experimental  confirmation. 

To  extend  the  Penning-ionization  technique  as  a  diagnostic  for  N2(v) 
containing  energies  of  several  eV,  we  investigated  excitation  of  the  nitrogen 
second-negative  system,  N2+(C2EU+  -  X2Eg+).  This  system  can  be  excited  by 
collisions  between  metastable  helium  and  molecular  nitrogen  only  if  the  nitro¬ 
gen  has  at  least  3.8  eV  of  internal  energy.  Because  the  N2+(C)  potential 
curve  is  somewhat  offset  from  that  of  N2(X),  vertical  transitions  connecting 
the  two  states  would  arise  from  highly-excited  ground-state  vibrational  levels. 
Observing  the  second-negative  system,  therefore,  would  provide  unequivocal 
evidence  of  highly  excited  ground-state  nitrogen  in  the  reactor. 

Using  a  quartz  flow  reactor  we  observed  the  second-negative  system  in  the 
ultraviolet  between  185  and  210  nm.  Figure  9  shows  the  spectral  region 
between  188  and  208  nm  with  the  metastable  helium  off  and  on.  In  both  cases, 
the  active-nitrogen  discharge  is  on.  The  two  features  appearing  at  191  and 
199  nm  are  the  Av  =  -6  and  -7  sequences,  respectively,  of  the  nitrogen  second¬ 
negative  system.  The  Av  =  -8  sequence  at  206  nm  is  masked  by  the  2,0  band  of 
the  N0(A2E+  -  X2n)  system.  This  system  is  excited  by  interactions  between 
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N2(A^Eu+),  created  in  the  active-nitrogen  discharge,  and  some  impurity  NO  also 
created  in  the  discharge^. 

Scans  with  the  active  nitrogen  discharge  off,  but  the  metastable-helium 
discharge  on,  resulted  in  no  observed  emissions.  This  latter  observation  is 
important  because  it  establishes  the  absence  of  He+  from  the  metastable  helium 
flow.  Charge- transfer  between  He+  and  N2(X)  is  a  well-known  source  of  the 
second-negative  system^3*51- 

Observing  the  nitrogen  second-negative  system  from  the  interaction 
between  metastable  helium  atoms  and  active  nitrogen  demonstrates  the  presence 
of  metastable  nitrogen  in  the  afterglow  containing  at  least  3.8  eV.  In  addi¬ 
tion  to  N2(X,v),  the  electronic  metastables  N2(A^EU+)  or  N2(a'^Eu“)  are  sig¬ 
nificantly  populated  in  active  nitrogen  afterglows.  Using  alternative  sources 
to  generate  the  electronic  metastables,  we  observed  the  resulting  spectra  when 
they  were  mixed  with  helium  metastables. 

The  N2(A^Eu+)  is  produced  cleanly  by  adding  N2  downstream  from  a  hollow- 
cathode,  d.c.  discharge  through  a  mixture  of  helium  and  argon-2'53.  Scans  of 
the  Vegard-Kaplan  bands  between  220  and  400  nm  established  that  the  N2(A)  num¬ 
ber  densities  from  this  source  were  of  comparable  magnitude  to  those  obtained 
from  the  active-nitrogen  discharge.  Turning  the  metastable  nitrogen  on  and 
off  in  the  presence  of  metastable  helium  showed  no  detectable  change  in  the 
emissions  in  the  region  between  185  and  210  nm. 

N2<a')  is  made  along  with  slightly  enhanced  number  densities  of  N2(A)  - 
about  a  factor  of  3  -  when  the  molecular  nitrogen  flows  through  the  d.c.  dis 
charge  with  the  He/Ar  mixture^ > 54.  Again,  turning  the  metastable  nitrogen  on 
and  off  indicated  no  emissions  between  185  and  210  nm.  Ue  conclude  therefore 
that  interactions  between  metastable  helium  and  N2<A)  and  N2(a')  do  not  pro¬ 
duce  N2+(C). 
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The  appearance  of  N2+(C)  from  the  interaction  between  metastable  helium 
and  active  nitrogen,  therefore,  appears  to  be  the  result  of  Penning  ionization 
of  ground-electronic-state  nitrogen  containing  at  least  3.8  eV  of  vibrational 
energy.  This  corresponds  to  ground-electronic-state  vibrational  levels  of 
v"  >  15. 


V.  Summary 

In  summary,  Penning  ionization  of  N2(v)  by  metastable  helium  atoms 
excites  N2+(B)  up  to  at  least  v'  =  8.  Analysis  indicates  the  ground-state 
nitrogen  vibrational  distribution  follows  a  Treanor  rather  than  Boltzmann 
distribution,  with  effective  vibrational  temperatures  of  up  to  6,000  K  at 
times  up  to  30  ms  downstream  from  the  active-nitrogen  discharge. 


The  diagnostic  is  complicated  by  the  presence  of  free  electrons  created 
in  the  Penning  ionization.  These  free  electrons  absorb  energy  from  the 
active-nitrogen  medium  and  produce  further  excitation  of  N2+(B).  Adding 
traces  of  SF^  removes  the  free  electrons  from  the  reactor,  and  thereby 
eliminates  their  interfering  effects. 


Observing  emission  from  N2+(C),  excited  by  Penning-ionizat ion 
in  the  afterglow,  demonstrates  the  presence  of  N2(v)  containing  at 
3.8  eV  internal  energy.  Future  efforts  will  be  directed  towards  a 
quantitative  assessment  of  the  N2+(C)  observations  to  identify  the 
N2<v)  levels  responsible  for  this  excitation. 
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TABLE  1.  Franck-Condon  Factors  of  N2+(B 


FIGURE  CAPTIONS 


Figure 

1.  Flow  reactor  for  studies  on  the  vibrational  energy  content  of  active 
nitrogen. 

2.  Ratio  of  the  populations  of  N2+(B,  v')  to  N2+(B,  v'  =  0)  created  in  He* 
Penning-ionization  of  N2(X,  v" )  as  a  function  of  the  vibrational 
temperature  of  the  ground-state  nitrogen. 

3.  Spectrum  of  the  Av=- 2  sequence  of  the  N2+(B^Iu+-x2£g+)  system  excited 
in  the  Penning-ionization  of  vibrationally  cold  nitrogen  by  He*(2 
The  light  line  is  the  experimental  spectrum;  the  heavy  line  shows  the 
best-fit  synthetic  spectrum. 

4.  Spectrum  of  the  Av=-2  sequence  of  the  N2+(B  2 Eu+-x2Eg+)  system  excited 
in  the  Penning-ionization  of  active  nitrogen  by  He*(2  ^S).  The  light 
line  is  the  experimental  spectrum;  the  heavy  line  shows  the  best-fit 
synthetic  spectrum. 

5.  Spectra  of  the  Av=-2  sequence  of  the  N2+(b2 Eu+-x2£g+)  system  excited  in 
the  Penning-ionization  of  active  nitrogen  by  He*(2  ^s)  in  the  presence 
and  absence  of  SFg. 

6.  Comparison  between  experimental  and  calculated  vibrational 
distributions  of  N2+(B)  created  in  the  Penning-ionization  of  active 
nitrogen  by  metastable  helium  atoms  for  a  nitrogen  mole  fraction  of 
0.011  (p  =  1.5  torr,  transit  time  from  discharge  =  11  ms). 

7.  Comparison  between  experimental  and  calculated  vibrational 
distributions  of  N2+(B)  created  in  the  Penning-ionization  of  active 
nitrogen  by  metastable  helium  atoms  for  a  nitrogen  mole  fraction  of 
0.046  (p  =  1.5  torr,  transit  time  from  discharge  =  11  ms). 

8.  Effective  vibrational  temperature  from  the  non-equilibrium  model  versus 
nitrogen  mole  fraction  through  the  discharge. 

9.  Spectra  from  active  nitrogen  between  184  and  208  nm  in  the  absence  and 
presence  of  metastable  helium  atoms. 
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Figure  1.  Flow  reactor  for  studies  on  the  vibrational  energy  content 
of  active  nitrogen. 
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Ratio  of  the  populations  of  N2+(B,  v')  to  N2+(B,  v'  =  0)  created 
in  He*  Penning-ionizat ion  of  N2<X,  v")  as  a  function  of  the 
vibrational  temperature  of  the  ground-state  nitrogen. 
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Figure  3.  Spectrum  of  the  Av=-2  sequence  of  the  N2+(B2Eu+-X^Zg+)  system 

excited  in  the  Penning-ionization  of  vibrationally  cold  nitrogen 
by  He*(2  ^S).  The  light  line  is  the  experimental  spectrum;  the 
heavy  line  shows  the  best-fit  synthetic  spectrum. 
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Figure  4.  Spectrum  of  the  Av=-2  sequence  of  the  N2  +  (B  ^ Eu  +  -X^ )  system 
excited  in  the  Penning-ionizat ion  of  active  nitrogen  by 
He*(2  ^S).  The  light  line  is  the  experimental  spectrum;  the 
heavy  line  shows  the  best-fit  synthetic  spectrum. 
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Figure  6.  Comparison  between  experimental  and  calculated  vibrational 
distributions  of  NpMB)  created  in  the  Penning- ionizat ion 
of  active  nitrogen  by  metastable  helium  atoms  for  a  nitrogen 
mole  fraction  of  0.011  (p  -  1.5  torr,  transit  time  from 
discharge  -  11  ms). 
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Figure  7.  Comparison  between  experimental  and  calculated  vibrational 
distributions  of  N2+(B)  created  in  the  Penning-ionization 
of  active  nitrogen  by  metastable  helium  atoms  for  a  nitrogen 
mole  fraction  of  0.046  (p  =  1.5  torr,  transit  time  from 
discharge  =  11  ms). 
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Figure  8.  Effective  vibrational  temperature  from  the  non-equilibrium  model 
versus  nitrogen  mole  fraction  through  the  discharge. 
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ABSTRACT 


We  have  studied  the  kinetics  of  the  lowest  excited  sinylet  state  of 
molecular  nitrogen,  ^(a'1^,  v=0),  in  a  discharge  flow  reactor.  The  meta- 
stables  were  generated  in  a  hollow  cathode  dc  discharge  through  molecular 
nitrogen  highly  diluted  in  argon,  and  detected  by  vuv  fluorescence  of  the  for¬ 
bidden  N2(a,1Eu  -  X^Eg)  band  system.  Observations  of  the  equilibrium  between 
the  N2<a^r[g)  and  Ngta'1^)  states  in  our  reactor  indicate  that  the  metastable 

has  a  radiative  lifetime  >  23  ' '  ms,  assuming  a  radiative  lifetime  of 

~  -  6 

80  US  for  N2(a1ITg)-  Rate  coefficients  for  quenching  by  NO,  CH4 ,  CO,  and 

N20  are  approximately  gas  kinetic,  those  by  H2  •  °2 '  aru3  co2  rouyhly  tenth  gas 
kinetic,  and  that  by  molecular  nitrogen  (1.9  r  0.5)  x  10-13  cm3 
molecule-1  s-1.  The  interaction  between  N2(a')  and  CO  produced  COfA1:]  -  X1^) 

emission  in  about  21  +  ^  percent  of  the  quenching  events. 

—  8 


1 .  INTRODUCTION 


Although  the  lowest  excited  singlet  state  of  molecular  nitrogen, 
N2(a,1Zu_)»  is  potentially  an  important  energy  carrier  in  systems  involving 
excited  nitrogen,  such  as  nitrogen  discharges  and  aurorae,  it  has  been  studied 
only  sparingly.  This  is  probably  because  its  long  radiative  lifetime,  esti¬ 
mated  to  be  between  13  and  500  ms,  makes  it  difficult  to  detect  optically. 
Several  authors  have  shown  that  it  probably  is  the  state  to  which  the  much 
more  studied  ^(a1^)  state  is  quenched. 

Golde^  developed  a  discharge  technique  for  studying  N2(a',  v'=0).  He 
showed  that  this  state  was  coupled  kinetically  to  N2(a,  v*=0)  which  is 
1212  cm-'  higher  in  energy.  The  fact  that  he  could  see  none  of  the  higher 
vibrational  levels  of  N2(a')  he  took  as  evidence  that  the  a  and  a'  states  were 
efficiently  coupled  throughout  the  nested  manifold,  and  that  during  the  few 
millisecond  transit  time  between  his  discharge  and  his  observation  vessel  all 
of  the  higher  lying  a' -state  vibrational  levels  had  relaxed  to  a-state  vibra¬ 
tional  levels.  Radiation  or  further  quenching  deactivated  these  levels  so 
that  some  fraction  of  them  ultimately  ended  in  a',  v'=0.  Later,  Casassa  and 
Golde4  studied  the  spectroscopy  of  ^(a1,  v'=0),  determining  the  transition 
moment  variation  with  r-centroid,  and  reanalyzed  the  absorption  data  of 
Tilford  and  Benesch^  to  estimate  a  radiative  lifetime  of  about  17  ms  for 
N2(a',  v'=0).  In  comparing  the  ratio  of  intensities  of  the  a  and  a'  emis¬ 
sions,  they  estimated  a  lower  limit  to  the  a'  lifetime  of  about  31  ms. 

We  have  studied  the  kinetics  of  reactions  of  N2(a',  v'=0)  with  a  number 
of  molecules  in  a  discharge-flow  reactor.  A  cold-cathode  discharge  source 
similar  to  Golde 1 s  produced  the  metastables.  Our  observations  confirm  those 
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of  Casassa  and  Golde  regarding  the  transition-moment  variation  of  ^(a1,  v'=Ci) 
and  the  equilibrium  ratio  of  a  to  a'  emission.  We  have  further  determined  the 
rate  coefficients  for  quenching  ^{a1,  v'=0)  by  a  variety  of  molecules,  and 
have  determined  the  efficiency  with  which  ^(a1)  excites  the  CO(A  -  X)  bands 
in  an  energy- transfer  reaction  with  CO. 
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2.  EXPERIMENTAL 


The  apparatus  used  in  these  studies  has  been  described  in  some  detail  in 
several  previous  publications.^-1^  Basically  it  is  a  2  in.  diameter 
discharge-flow  reactor  which  incorporates  a  moveable  injector  for  varying  the 
reaction  time  between  the  metastables  and  the  quenchers.  A  cold-cathode 
discharge  operating  at  about  390V  through  a  mixture  of  molecular  nitrogen 
highly  dilute  in  argon  produces  the  metastables.  Observation  of  fluorescence 
from  the  forbidden  N2(a',  v'=0  -  X,  v"=4)  transition  at  171  run  was  our  primary 
diagnostic  for  the  metastable  number  density.  Because  molecular  nitrogen  is  a 
reasonably  efficient  quencher  of  ^(a1),  the  mole  fraction  of  nitrogen  in  the 
discharge  and  the  total  system  pressure  must  be  kept  low,  and  transit  times 
between  the  discharge  and  the  observation  region  need  to  be  minimized. 
Typically  we  operated  at  a  pressure  of  about  0.95  torr  with  a  nitrogen  mole 
fraction  of  about  0.007  and  a  total  transit  time  between  the  discharge  and  the 
detection  region  of  about  13  ms. 

Although  the  discharge  source  produces  some  atomic  nitrogen,  both  N(4S) 
and  N(^D),  and  N2(A^j;u+),  none  of  these  species  should  affect  our  kinetic 
results.  In  principle,  ^(A)  could  be  collisionally  excited  to  ^(a1)  in 
energy  pooling  reactions  involving  other  metastables  which  supply  at  least 
2.23  eV.  This  means  N(^D),  N(^P),  ^(A),  and  ^(X,  v>8)  .  The  absence  of 
emission  at  346.6  nm,  N(^P-^S),  shows  that  N(^p)  number  densities  are  less 
than  106  atoms  cm-^.  At  higher  pressures  {>  5  torr)  and  longer  flow  times 
from  the  discharge,  we  do  see  some  (N^p)  emissions  from  this  source. 
Presumably,  the  N(2p)  j.s  excited  by  the  efficient  energy  transfer  between 
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N2(A)  and  N(4S). 14,18  Metastable  helium  Penning  ionization  measurements  on 
the  effluents  of  this  discharge  source  showed  that  N2(X,v)  also  is  absent  from 
the  gas  stream.16  Our  measurements  on  N2(A)  energy  pooling  indicated  very 
little,  if  any,  singlet  metastable  formation.17  Resonance-fluorescence 
measurements  in  the  vacuum  ultra-violet,  however,  do  indicate  the  presence  of 
N(2D)  metastables  at  number  densities  which  could  be  as  high  as 
1018  atoms  cm-2.  Thus  we  cannot  so  easily  rule  out  the  possibility  of  energy 
pooling  between  N2(A)  and  N(2D)  as  a  source  of  N2(a‘).  Our  icinetic  obser¬ 
vations,  however,  seem  to  indicate  that  this  is  not  the  case.  For  example, 

C02  quenches  N2(a')  efficiently  but  quenches  N2(A)  poorly18  and  reacts  fairly 
slowly  with  N(2d).12  Another  possible  downstream  source  of  N2(a')  could  be 
the  three-body  recorabina tion  of  atomic  nitrogen.8  However,  N-atom  number 
densities  from  this  discharge  source  are  below  1012  atoms  cm"8,  a  number 
density  too  small  to  produce  significant  N2(a')  emission  via  three-body 
recombination  at  pressures  below  a  torr .  We  verified  this  point  by  failing  to 
find  N2(a')  in  the  afterglow  of  a  microwave-discharge  source  of  atomic 
nitrogen  with  atom  number  densities  on  the  order  of  10 12  atoms  cm-2  and  total 
pressures  below  2  torr. 
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3 .  RESULTS 


A .  Spectral  Characterization  and  radiative  lifetime  of  N?(a'^En~»  v‘=0) 

Figure  1  shows  the  spectrum  observed  in  our  reactor.  The  only  signifi¬ 
cant  emissions  between  140  and  180  nm  are  transitions  from  the  lowest  vibra¬ 
tional  level  of  the  a  and  a'  states.  We  estimated  a  relative  monochromator 
response  function  from  the  intensities  of  the  various  Lyman-Birge-Hopf ield 
bands  (^(a  -  X)  ;  LBH)  assuming  a  constant  transition  moment  for  the  a  -  X 
emission. ^  with  this  response  function,  we  could  then  estimate  the  relative 
transition-moment  variation  for  the  Ogawa ,  Tanaka,  Wilkinson,  Mulliken  bands 
(N^la1  -  X);  OTWM )  using  the  Franck  Condon  Factors  given  in  Lofthus  and 
Krupenie . Not  surprisingly,  our  results  agree  quite  well  with  those  of 
Casassa  and  Golde  who  used  a  similar  procedure.  The  relative  response  func¬ 
tion,  and  the  transition  probabilities  calculated  with  the  transition  moment 
variation  we  had  determined  allow  us  to  fit  the  two  band  systems  with  our 
synthetic  spectral  fitting  code . ^ 1  The  heavy  line  in  Figure  1  shows  the 
results  of  the  fit  compared  to  the  experimental  observations  which  are  shown 
with  the  light  line.  Spectral  fits  to  a  number  of  spectra  showed  that  the 
ratio  of  the  intensities  of  the  LBH  to  the  OTWM  bands  is  1.71  ±  0.05  at  300K. 

Golde^  has  already  shown  that  this  ratio  is  independent  of  argon  pressure 
over  the  range  0.5  to  4  torr.  Our  observations  show  that  quenching  the  bands 
significantly  by  adding  N2O  and  CO2  also  does  not  change  this  ratio.  The  two 
emissions  are  therefore  in  a  steady-state  ratio  in  our  reactor.  We  can  use 
this  observation  to  estimate  a  lower  limit  to  the  lifetime  of  the  N2(a') 
state.  We  need  consider  only  the  following  three  reactions: 


( 1) 


N2(a)  is  in  steady  state  within  our  field  of  view  because  its  radiative  life¬ 
time  is  only  80  ys.21~27  <phe  kinetic  rate  expressions  derived  from  reac¬ 
tions  (1)  through  (3)  can  be  combined,  therefore,  to  give  the  ratio  of  number 
densities  of  the  two  singlet  metastables: 

[ N2 (a )  ]  k 1 [ Ar ] 

[N2(a')]  =  (k_1  +  k2)[ArJ  +  k3  *  U) 

Since  the  ratio  of  the  a  to  a'  emission  is  independent  of  pressure, ^  quenching 
must  dominate  radiation  and  the  radiative-decay  rate,  k3 ,  can  be  dropped  from 
Eq.  (4) .  The  number  density  is  the  observed  intensity  multiplied  by  the 
radiative  lifetime.  The  ratio  of  the  rate  coefficients  k]/k_1  is  the 
Boltzmann  equilibrium  ratio, 


k1  „  -1212  cm” 1/kT 

- -  =  2  e  , 

k-1 

where  the  factor  of  two  accounts  for  the  lambda  doubling  in  the  pi  state. 
With  these  last  three  points  in  mind,  we  can  rearrange  Eq.  (4)  to  get 
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la  the  limit  that  deactivation  to  the  a'  state  is  the  only  significant  pathway 
for  a-state  quenching,  we  can  set  a  limit  on  the  radiative  lifetime  of  the  a' 
state  in  terras  of  that  of  the  a  state: 


r 


a  1 


(7) 


The  radiative  lifetime  of  ^(a)  is  still  somewhat  in  doubt,21-27  but  a  value 
+  40 

of  30  Us  apparently  encompasses  the  correct  value.  Thus  we  find  that 


t  a  ' 


■+* 

>  23 


(8) 


This  analysis  is  similar  to  that  given  by  Casassa  and  Golde.4  They  used  a 
different  lifetime  for  ^(a1^),  however,  and  reported  ta '  >  32  ms.  Their 
reanalysis  of  the  experiments  of  Tilford  and  Benesch3  yielded  a  lifetime  of 
17  ms.  Thus  they  concluded  that  the  true  lifetime  probably  lay  between 
17  and  40  ras .  We  feel  that  the  analysis  of  the  ratio  of  a  to  a'  intensities 
should  give  a  more  accurate  assessment  of  the  radiative  lifetime  than  does  the 
vuv  absorption  measurement  of  Tilford  and  Benesch  which  involves  a  comparison 
of  densitometer  tracings  o  f  a'l£u-  and  A3Eu+  absorptions  rather  than  a  direct 
curve-of-growth  measurement.  The  major  uncertainty  in  the  present  approach  is 
the  radiative  lifetime  of  r^Ca  ilg)  .  Clearly  what  is  needed  is  a  good  deter-  \ 

mination  of  that  value. 

B .  Quenching  Rate  Measurements 

The  reactions  important  in  our  kinetic  analysis  involve  quenching  of  the 
^(a1)  by  the  species  of  interest,  Q,  diffusion  to  the  walls  of  the  reactor 
with  subsequent  deactivation,  and  quenching  by  the  N2  and  Ar  bath  gas: 


* 

N 

2 


/.AH 


*-  Q  +  N  (X)  +  Q 

2 


(9) 


( 10) 


t  k 

+  wall  ->  N^fX)  +  wall 

N2*  +  Xj.  -*•  N2(X)  +  Xj.*  ,  (11) 

where  X^  represents  other  quenching  species  such  as  N2,  Ar,  etc.  The 
differential  equation  describing  the  rate  of  change  in  the  number  density  of 
the  metastables  with  time  is 


dt 


(  1.2) 


Because  the  number  density  of  the  quenchers  is  much  larger  than  that  of  the 
metastables,  it  remains  constant  (it  would  in  any  event  if  the  reaction 
involves  energy  transfer  rather  than  reaction).  Thus  we  can  solve  Eq.  (12): 


Keeping  the  reaction  time  fixed  and  varying  the  number  density  of  the  quencher 
gives  the  decay  coefficient  ,r,  defined  by 


d  Zn[N*] 

d[Q]  =  '  k9  Z/*  =  ’  (14) 

where  the  reaction  time  is  given  by  the  ratio  of  the  distance  from  the  injec¬ 
tor  to  the  observation  region,  z,  to  the  bulk  flow  velocity,  v.  Finally,  the 
rate  coefficient  results  from  the  slope  of  a  plot  of  f  versus  z  : 


k 


9 


(  15) 
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This  procedure  corrects  for  imperfect  mixing  at  the  injector.  Figure  2  shows 


semilog  plots  according  to  Eq.  (14)  for  the  quenching  of  N^Ca')  by  NO  at  three 
different  reaction  times.  Figure  3  shows  plots  of  decay  coefficient  versus  z 
for  three  different  quenchers. 

The  above  analysis  assumes  plug  flow.  The  rate  coefficients  derived  from 
plots  according  to  Eqs.  (14)  and  (15)  must  be  divided  by  0.62  to  correct  for 
the  coupling  of  a  radial  density  gradient  in  the  raetastable  number  density 
with  a  parabolic  velocity  profile.  A  number  of  papers  in  the  literature 
discuss  this  effect  in  detail. 28-31 

Wall  deactivation  and  relatively  efficient  quenching  by  molecular 
nitrogen  constrained  the  experimental  conditions  to  fairly  low  pressures  and 
short  reaction  times.  Several  of  the  most  inefficient  quenchers,  therefore, 
required  fairly  large  additions  of  gas  to  the  reactor  in  order  to  obtain 
adequate  decays.  In  order  to  correct  for  any  possible  perturbation  to  the 
metastable  number  densities  caused  by  injecting  large  mass  flows  of  gas  into 
the  main  flow,  we  measured  apparent  decays  when  comparable  amounts  of  argun 
were  added  through  the  injector.  We  then  corrected  the  measured  decay 
coefficients  of  the  quenchers  of  interest  by  subtracting  the  effective  argon 
decay  coefficients  from  them.  This  correction  procedure  assumes  that  argon  is 
an  inefficient  quencher  of  N2(a')«  The  decrease  in  the  effective  argon  decay 
coefficient  with  increasing  reaction  time  indicates  this  assumption  is  valid. 
The  molecules  N2,  CF4 ,  and  SFg  were  corrected  in  this  manner.  This  large 
addition  of  gas  was  still  <  7  percent  of  the  total  gas  flow,  so  it  did  not 
perturb  the  flow  field  sufficiently  to  change  the  flow  correction  factor. 

Table  1  summarizes  our  results.  The  error  bars  on  the  rate  coefficients 
include  the  contributions  from  the  standard  deviations  in  the  least  squares 


fits  according  to  Eqs.  (14)  and  (15),  our  estimates  of  the  uncertainties 
involved  in  the  calibrations  of  our  flow  meters,  pressure  gauge,  thermometer, 
etc.,  and  our  estimate  of  the  basic  reproducibility  of  our  measurements  as 
determined  in  other  kinetic  investigations  which  involved  multiple  Determin¬ 
ations  of  rate  coefficients  under  a  number  of  different  conditions  of  pres¬ 
sure,  total  gas  flow  rate  and  flow  velocity . ® ^ > 1 ^ 

C .  The  excitation  of  CO(A)  by  energy  transfer  from  N?(a') 

When  carbon  monoxide  is  added  to  the  metastables,  emission  from  the  CO 
fourth-positive  system,  CO(A'n  -  X^E),  is  observed.  Golde  and  Thrush^2>33 
have  shown  that  CO  fourth-positive  emission  is  excited  in  the  energy  transfer 
Detween  N2(a)  and  CO.  In  our  system,  however,  the  low  steady-state  number 
density  of  the  N2(a)  precludes  that  state  as  being  the  source  of  the  observed 
C0( A)  emission  (see  below).  We  conclude  therefore,  that  it  arises  from  energy 
transfer  from  N2(a'): 

CO{xh)  --  N2(a,1Eu")  *  C0(A]n)  +  N2(X1Eg+)  (16) 

C0(A1n)  +  COU1!)  +  hv  •  (17) 

The  extremely  rapid  radiative  lifetime  of  the  CO(A)34  prevents  its  being 
deactivated  by  electronic  quenching  under  our  conditions.  Thus  the  CO(A)  is 
in  steady  state  in  our  observation  region,  and  we  can  write 

ICOIAH  T  .  (18) 

IN2(a')  k17  a' 

We  scanned  the  spectrum  between  140  and  185  nm  for  varying  amounts  of 
added  CO.  We  determined  the  number  density  of  the  emitting  CO(A)  molecules 
and  the  N2(a')  and  N2(a)  emission  intensities  by  spectral  fitting.  Figure  4 


shows  one  such  spectrum  with  its  associated  fit.  The  da ta  derived  from  the 


fits  are  plotted  according  to  Eq.  (18)  in  Figure  5  for  vibrational  levels  0-2 


of  CG( A) .  Excitation-rate  coefficients  derived  from  the  slopes  of  the  lines  in 

Figure  5,  assuming  a  radiative  lifetime  for  N2(a')  of  23  +  '  ms ,  are  given 

—  8 

+  10 

in  Table  2.  The  results  indicate  that  only  21  percent  of  the  quenching 

—  8 

events  between  N2(a')  and  CO  result  in  CO  electronic  excitation. 

If  the  a-state  rather  than  the  a'  state  were  responsible  for  exciting 
C0(A,n) ,  then  our  data  would  indicate  an  excitation  rate  coefficient  of 
5  x  10-9  cm3  molecule-1  s-1,  more  than  ten  times  gas  kinetic.  Golde  and 
Thrush-1-1  previously  have  established  that  the  rate  coefficient  for  exciting 
CO( A)  by  N2(a)  is  only  1.5  x  10-1(1  cm1  molecule-1  s-1.  Thus  we  can  confi¬ 
dently  assign  N2(a')  as  the  precursor  for  C0(A)  excitation  in  our  system. 
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4. 


DISCUSSION 


Dreyer  and  Perner^  excited  N2(a')  with  a  pulsed  relativistic  electron 
beam  in  nitrogen  cas  and  observed  its  temporal  behavior  via  resonance 
absorption  on  the  fifth-positive  system  of  nitrogen.  They  analyzed  their  data 
to  determine  a  rate  coefficient  for  N2(a')  quenching  by  molecular  nitrogen  of 
(2.21  ±  0.14)  x  10-1®  cm3  molecule-1  s-1  in  excellent  agreement  with  our 
value.  Van  Veen  et  al.^  used  a  laser  to  excite  N^a1^,  v'=0,1)  via  two- 
photon  pumping.  The  variation  in  the  fluorescence  decay  lifetimes  with 
nitrogen  pressure  over  the  range  of  0.05  to  10  torr  indicated  complete  cou¬ 
pling  between  N2(a)  and  N2(a')  at  pressures  above  1  torr,  so  that  observed 
decays  in  that  pressure  region  were  diagnostic  of  N2(a')  behavior.  They 
reported  a  rate  coefficient  for  quenching  N2(a')  by  nitrogen  of  (2.3  ±  0.2)  x 
10-1®  cm3  molecule-1  s-1 ,  again  in  excellent  agreement  with  the  present 
results . 

Golde  and  Thrush®  observed  the  behavior  of  N2(a)  and  N2(a')  under  various 
conditions  in  a  flow  of  recombining  nitrogen  atoms  and  concluded  that  ratio  of 
the  rate  coefficients  for  N2(a')  quenching  by  Ar  to  that  by  N2  was  unity  and 
that  the  similar  ratio  for  C02  quenching  to  that  by  N2  was  117  ±  20.  Our 
results  clearly  indicate  that  quenching  of  N2(a')  by  argon  must  be  much  less 
efficient  than  quenching  by  nitrogen.  If  the  argon  quenching  rate  coefficient 
were  the  same  as  our  nitrogen  value,  we  calculate  that  the  raetastables  would 
have  a  lifetime  against  quenching  of  0.28  ms  at  the  argon  number  density  in 
our  reactor.  Thus  it  would  not  survive  the  13  ms  traversal  time  from  the 
discharge  to  the  observation  region.  The  fact  that  N2(a')  is  observed  in  the 
observation  region  shows  that  N2(a')  quenching  by  Ar  is  more  than  an  order  of 
magnitude  slower  than  is  quenching  by  N2 .  Our  results  give  a  ratio  of  the  C02 
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rate  coefficient  to  that  for  N2  of  129  t  21  in  good  agreement  with  Golde  and 
Thrush's  observation. 

We  are  unaware  of  any  other  quenching  data  on  N2(a').  Our  rate 
coefficient  for  quenching  by  molecular  oxygen  indicates  that  quenching  would 
be  comparable  to  radiative  decay  at  an  altitude  just  below  105  km  (assuming  a 
radiative  lifetime  of  23  ms).  Thus  one  cannot  neglect  the  importance  of 
N2(a')  as  a  metastable  energy  carrier  in  modeling  energy  degradation  in  the 
disturbed  upper  atmosphere. 
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Table  1.  Results  of  ^(a'^Eu-)  Quenching  Studies 


MOLECULE 

QUENCHING  RATE  COEFFICIENT 
(cm3  molecule"^  s~^) 

N2 

(1.9  ±  0.5)  x  10~13 

°2 

(2.8  ±  0.6)  x  10“n 

co2 

(2.5  ±  0.6)  x  10-1  1 

NO 

(3.6  ±  0.8)  x  10-10 

n2° 

(1.7  ±  0.4)  x  10" 10 

ch4 

(3.0  ±  0.8)  x  10-10 

CO 

(1.1  i  0.3)  x  10~ 10 

H2 

(2.6  ±  0.6)  x  10-11 

cf4 

(2.4  i  1.0)  x  10- 13 

sf6 

(7.8  ±  1.8)  x  10~14 

* 
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Table  2.  Results  of  CO(a'h)  Excitation  by  N2(a 


(cm2  molecule-'  s~')* 


0 

(11.7  +  x  10“ 12 

-  4.4 

1 

(6.8  +  3-l)  x  10-12 

“ “  4  •  J 

2 

(3.8  +  '•])  x  10-12 
-  1.4 

3 

(0.8  +  °'5)  x  1 0“ 1 2 

-  0.4 

TOTAL 

(23.1  +  ic  10-U 

—  o  .  / 

‘Assumes  lifetime  of 

23  +  ^  ms  for  N2(a,1Eu-) 

5 


Figure  1 


Figure  2 

Figure  3 


Figure  4 


Spectrum  of  N2(a',Eu”)  and  ^(a^rig)  13  ms  downstream  from  the 
discharge.  The  light  line  shows  the  experimental  spectrum, 
the  heavy  line  the  synthetic  best  fit  to  the  data. 

The  decay  of  ^(a'1^-)  with  added  [NO]  for  three  different 
reaction  times. 

Decay  coefficients  plotted  against  reaction  time  for  N2(a‘1Eu“) 
quenching . 

Spectrum  of  COCA1  II)  excited  by  energy  transfer  from  N2(a  '  1  Eu_)  . 
The  experimental  spectrum  is  the  light  line,  the  synthetic  best 
fit  to  the  experimental  data  is  the  heavy  line. 


452 


0.8 


0.6 


t  0.4 
co 


0.2 


0.0 


140  150  160  170  180 

WAVELENGTH  (nm) 

A-3281 

FIGURE  1 


453 


INTENSITY  (arbitrary  units) 


1.0 


140  150  160  170  180 

WAVELENGTH (nm) 

I  I  Cl -RE  4 


4  r)6 


A-3282 


APPENDIX  M 


(SR-353  reproduced  in  its  entirety) 


J.  Chem.  Phys.  90,  XXXX  (1989) 


457/458 


PSI— 9057 
SR-353 


EXPERIMENTAL  DETERMINATION  OF  THE  EINSTEIN 
COEFFICIENTS  FOR  THE  N2(B-A)  TRANSITION 


Lawrence  G.  Piper,  Karl  V.  Holtzclaw,  and  B.  David  Green 
Physical  Sciences  Inc. 

Research  Park,  P.O.  Box  3100 
Andover,  MA  01810 


and 


Villiam  A.M.  Blumberg 
Air  Force  Geophysics  Laboratory/LSI 
Hanscom  Air  Force  Base,  MA  01731 


January  1989 


ABSTRACT 


We  have  used  a  branching-ratio  technique  to  measure  the  relative 
variation  in  the  transition-dipole  moment  with  internuclear  separation  for  the 
N2<B-A)  transition.  Our  spectral  observations  cover  the  range  from  500  to 
1800  nm,  and  use  several  different  detectors  and  excitation  sources.  The  data 
from  different  sets  are  consistent  in  the  regions  of  spectral  overlap.  Using 
well  established  values  for  the  radiative  lifetimes  of  N2(B,v'>5)  allows  the 
relative  dipole-moment  function  to  be  placed  on  an  absolute  basis.  From  the 
dipole-moment  function  and  a  set  of  RKR-based  Franck-Condon  factors  which  we 
have  computed,  we  derive  Einstein  coefficients  covering  the  range  v'=0-l2  and 
v"=0-20.  Our  results  indicate  that  currently  accepted  lifetimes  for 
N2(B,v'=0-2)  should  be  revised  upwards  by  20  to  40  percent. 
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I.  INTRODUCTION 


The  currently  accepted  experimental  values  for  the  N2(B-A)  Einstein 
coefficients  (i.e.,  as  tabulated  in  Lofthus  and  Krupenie^-)  rest  upon  the 
variation  in  the  electronic  transition  moment  as  a  function  of  internuclear 
separation  determined  by  Shemansky  and  Broadfoot.^  Using  a  branching-ratio 
technique  and  the  r-centroid  approximation, 3  they  determined  the  relative 
variation  in  the  transition  moment  with  internuclear  separation  between  1.20 
and  1.55A.  They  placed  this  relative  curve  on  an  absolute  basis  by  adjusting 
it  to  give  a  lifetime  for  v'=3  which  agreed  with  one  they  determined  directly 
by  measuring  real-time  fluorescence  decays  following  excitation  of  N2  by  a 
pulsed,  electron  beam. 

This  procedure  gave  lifetimes  for  ^(B^llg)  in  modest  accord  (±15  percent) 
with  those  derived  from  most,^-6  but  not  all,7  recent  theoretical  transition- 
moment  functions  for  v'>5.  The  theoretical  lifetimes  all  become  increasingly 
longer  than  those  derived  from  Shemansky  and  Broadfoot's  work,  however,  for  the 
lowest  vibrational  levels.  This  discrepancy  reaches  50  percent  for  v'=0.  Even 
for  the  higher  vibrational  levels,  where  the  agreement  in  the  total  radiative 
lifetimes  between  theory  and  experiment  is  adequate,  the  branching  ratios  for 
radiation  from  a  given  upper  vibrational  level  to  the  various  lower  levels 
differ  greatly  between  theory  and  experiment. 

Shemansky  and  Broadfoot's  Einstein  coefficients  for  the  lowest  vibra¬ 
tional  levels  of  N2(B)  result  from  extrapolations  of  their  transition-moment 
curve  well  outside  the  region  of  the  measurements  used  to  establish  it.  This 
procedure  is  prone  to  systematic  error.  To  reconcile  experiment  and  theory  we 
have  redetermined  the  relative  transition-moment  variation  as  a  function  of 
internuclear  separation,  using  a  procedure  similar  to  that  used  by  Shemansky 
and  Broadfoot,  but  with  measurements  extending  over  a  much  greater  range,  1.15 
to  1 . 70A. 


The  first-positive  bands  arise  from  the  ^(B^rig  -  A^EU+)  transition,  and 
are  among  the  most  prominent  in  nitrogen  and  air  discharges  of  various  types, 
in  the  disturbed  upper  atmosphere, ® ^  and  in  a  number  of  chemiluminescent 
reactions . 10-14  jn  addition,  laser  oscillation  has  been  observed  from  these 
bands. 16  In  order  to  understand  excitation  mechanisms  of  the  various  first¬ 
positive  sources  or  to  quantify  first-positive  band  lasing,  one  must  have 
accurate  Einstein  coefficients.  Furthermore,  a  complete  set  of  Einstein 
coefficients  allows  one  to  predict  emissions  accurately  in  one  region  of  the 
electromagnetic  spectrum  based  upon  observations  in  another  region, 
e.g.,  infrared  radiances  can  be  estimated  from  visible  intensities. 
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II.  THEORETICAL  BASIS  FOR  MEASUREMENTS 


The  r-centroid  approximation  was  first  put  forth  by  Fraser^  and  has  been 
used  extensively  over  the  last  three  decades  to  evaluate  the  transition-moment 
function  of  a  large  number  of  molecular  systems  experimentally.  Although  this 
method  has  been  criticized,  several  authors  have  demonstrated  that  it  is 
generally  accurate  for  most  molecules  having  only  slowly  varying  transition- 
moment  functions-*- ^-19  (hydrides  and  deuterides  appear  to  be  an  exception). 

The  intensity  of  a  molecular  emission,  in  units  of  photons  cm“3  s~-*,  is 
given  by 


v'  V' 


=  N  A 


v'  v' v' 


(1) 


where  Nv/  is  the  number  density  in  the  upper-state  vibrational  level  and  Av>vm 
is  the  Einstein  coefficient, 
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v '  v 


Vv" 


Th~ 


[  <v 7 |Re(r) !v"> 


(2) 


Here  is  the  transition  frequency  in  cm'*,  the  integral  is  of  the  wave- 

functions  over  the  electric  dipole-moment  operator,  and  the  constants  have 
the’  usual  meanings.  This  matrix  element  usually  can  be  separated  into  two 
parts.  One  part  represents  the  overlap  of  the  wavef unc t ions  of  the  two  levels 
with  each  other,  and  is  called  the  Franck-Condon  factor.  The  other  is  the 
electronic  transition  moment  which  is  a  function  of  internuclear  separation. 

Fraser  showed  that  the  electronic  transition  moment  could  generally  be 
represented  adequately  at  discrete  points  by  a  function  (usually  polynomial  or 
exponential)  of  the  r-centroid.  This  quantity  is  an  average  r  tor  each  (v'-*v" ) 
transition  and  is  defined  as 


A  f,  ) 


<v '  I r I v"> 
< v  | v"> 


(3) 


'v' v" 


By  measuring  the  relative  intensities  of  a  number  of  bands  with  a  common 
upper  vibrational  level,  one  can  map  out  relative  values  for  the  electronic- 
transition  moment  at  various  r-centroids,  and  thereby  determine  the  variation 
in  the  transition  moment  with  internuclear  separation.  Combining  Eqs.  (1)  and 
(2)  and  separating  the  dipole-moment  matrix  into  its  component  parts  gives 
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V' v" 


V'  V' 


IF 


N 


v  f  V  " 


lRe  (Vv") 


(4) 


where  the  Franck-Condon  factor  is 

Vv"  =  l<V>  |v">|2 


(5) 


The  ratio  of  measured  band  intensities  to  the  product  of  the  Franck-Condon 
factor  times  the  cube  of  the  transition  frequency  for  each  transition  then 
gives  a  set  of  reduced  intensities  which  should  vary  one  from  the  other  in  the 
same  way  that  the  electronic  transition  moment  varies  with  internuclear 
separation: 


I  ,  /q  ,  \>, 

V'  V  ^V  V  v' V 

: — ; - 3 — 

v'  v"/c*v'  v"  V  v" 


R  (r 
e 


(Vv 


R  (r  ,  , 
el  v' v1 


(6) 


This  technique  gives  only  the  relative  transition-moment  variation.  The 
transition-moment  function  is  made  absolute  by  normalizing  to  experimental 
lifetime  or  oscillator  strength  data. 


',66 


The  sum  of  the  Einstein  coefficients  from  a  given  upper  level  to  all 
accessible  lower  levels  equals  the  reciprocal  of  the  radiative  lifetime  of  that 
upper  vibronic  leve1 


E  A 


v'  v' 


=  1/T 


V' 


(7) 


Knowing  the  radiative  lifetimes  of  N2(B,v)  allows  one  to  put  the  relative 
transition-moment  function  on  an  absolute  basis  (vide  infra). 


Oscillator-strength  data  can  establish  one  Einstein  coefficient  from  the 
set  of  relative  Einstein  coefficients  absolutely: 


_  2  2 
8n  e 


v' v" 


m  c 
e 


"72  3”  v'v,f 

X  ,  „  u 
v'  vM 


(8) 


In  Eq.  (8)  me  is  the  electron  rest  mass,  e  its  charge  in  esu,  c  the  velocity 
of  light,  Xy'v"  the  transition  wavelength  in  cm,  du  and  d  the  electronic 
degeneracies  of  the  upper  and  lower  states,  respectively,  and  fVfVn  the 
oscillator  strength.  We  are  unaware  of  reliable  oscillator-strength  measure¬ 
ments  for  N2  first-positive  transitions. 

We  have  measured  the  intensities  of  a  number  of  first-positive  bands  over 
the  wavelength  region  between  500  and  1700  nm.  This  set  of  sequences  covers 
the  Av=6  through  Av=-2  sequences  and  samples  a  range  of  internuclear  separa¬ 
tions  between  about  1.15  and  1.7  A.  We  have  placed  our  relative  measurements 
on  an  absolute  scale  using  the  experimental  lifetimes  of  Eyler  and  Pipkin. 20 
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III.  EXPERIMENTAL 


The  experimental  procedure  involved  determining  first-positive  spectra 
over  the  appropriate  spectral  regions  under  several  different  sets  of 
conditions  using  monochromators  with  photoelectric  detectors.  Individual  band 
intensities,  Iv>vm,  for  each  spectrum  were  determined  using  a  computerized 
spectral-fitting  procedure  (see  below).  This  procedure  greatly  reduces 
uncertainties  introduced  by  spectral  overlap  with  other  first-positive  bands 
or  with  other  nitrogen  band  systems. 

Three  different  monochromator/detector  combinations  and  three  different 
excitation  sources  were  used  in  these  studies.  This  allowed  us  to  make 
multiple  determinations  of  the  data  in  the  regions  of  spectral  overlap.  In 
addition,  different  excitation  sources  emphasize  the  formation  of  different 
vibrational  levels.  This  allowed  matching  the  source  to  the  most  appropriate 
wavelength  region. 

One  approach  consisted  of  looking  at  the  excitation  of  nitrogen,  at  pres¬ 
sures  on  the  order  of  a  few  millitorr,  by  a  high  energy  electron  beam  of 
4.5  kV  and  about  15  mA.  The  apparatus  for  these  observations,  the  LABCEDE 
facility  at  the  Air  Force  Geophysics  Laboratory  (AFGL),  has  been  described  in 
detail  previously . ^1  Basically  it  consists  of  a  cylindrical  vacuum  chamber,  lm 
diameter  by  3.4m  long.  The  electron  beam  enters  perpendicular  to  the  chamber 
axis  about  1.2m  from  the  viewing  port  in  one  end.  Gas  flows  through  the  region 
irradiated  by  the  electron  beam  in  a  few  milliseconds.  This  prohibits  the 
buildup  of  beam-created  species  within  the  irradiated  volume.  A  0.3m 
monochromator  with  an  S-l  photomultiplier  monitored  fluorescence  in  the  chamber 
between  570  and  1050  nm.  This  spectral  region  encompasses  the  Av=4  through 
Av=0  sequences  of  the  nitrogen  first-positive  system. 

The  second  approach,  which  emphasized  the  short-wavelength  region  of  the 
spectrum,  involved  observing  atomic  nitrogen  recombination  in  a  He/N2  discharge 
afterglow.  This  source  ‘ends  to  emphasize  excitation  of  the  highest  levels  of 
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N2 ( B ) ,  v'=9-12.  A  0.5m  monochromator  coupled  to  a  thermoelectrically  cooled 
photomultiplier  with  a  GaAs  photocathode  was  used  for  these  measurements.  They 
encompassed  the  spectral  region  between  500  and  850  nm,  covering  principally 
the  Av=6  through  Av=2  sequences  of  the  first-positive  system. 

The  third  approach  consisted  of  infrared  observations  of  a  low  pressure 
(0.1-20  Torr),  low  power  (10-20  Watts)  microwave-discharge  lamp  in  nitrogen. 

In  these  studies  a  0.5m  monochromator  with  an  intrinsic  Ge  detector  scanned  the 
spectral  region  between  700  and  1700  nm.  This  covered  the  Av=2  through  Av=-2 
sequences  of  the  nitrogen  first-positive  system. 

All  systems  were  calibrated  for  relative  spectral  response  using  a 
standard  quartz-halogen  lamp.  The  lamp  and  its  power  supply  were  obtained  from 
Optronic  Laboratories.  When  set  up  as  prescribed  on  the  calibration  sheet,  the 
lamp  irradiance  is  certified  to  have  an  uncertainty  of  ±2  percent  between  500 
and  1600  nm,  and  ±4  percent  between  1600  and  2500  nm.  In  the  latter  two 
systems,  the  light  from  the  lamp  was  reflected  off  a  BaS04  screen  before 
passing  through  the  entrance  slit  of  the  monochromator.  The  BaS04  screen  is  a 
uniform,  diffuse  reflector.  Its  reflectivity  is  99  +  1  percent  between  300  and 
1000  nm,  and  declines  only  slightly  at  longer  wavelengths  (94  percent  at 
1700  nm).22,23  ye  did  correct  for  ..e  decrease  in  reflectivity  at  the  longer 
wavelengths . 

The  reason  for  using  the  BaS04  screen  is  to  ensure  the  monochromator 
optics  are  filled.  When  the  monochromator  views  the  lamp  directly,  only  a 
small  portion  of  the  grating  is  illuminated,  and  the  location  of  the  lamp  image 
on  the  grating  becomes  a  sensitive  function  of  lamp  placement.  Irradiaing 
small  grating  imperfections  could  distort  the  response  of  the  instrument  to  the 
radiation.  Under  filling  the  monochromator  optics  also  can  affect  the  location 
of  the  lamp  image  on  the  detector  element.  Any  non-uniformities  in  detector 
responsivity  with  position  across  the  detector  surface,  therefore,  could  lead 
to  a  distorted  response-function  calibration.  The  BaS04  screen  eliminates 
these  potential  difficulties. 
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We  confirmed  the  response  function  of  the  latter  two  systems  with  a  second 
quartz-halogen  lamp  (alro  from  Optronic  Laboratories)  and  with  a  1250  K  black 
body  (from  Infrared  Associates).  Relative  response  functions  derived  from  all 
three  standard  sources  agreed  to  within  +5  percent.  As  a  final  crosscheck  of 
the  relative  response  of  the  latter  two  systems,  we  compared  observed 
intensities  of  the  0/N0  air  afterglow  continuum  with  published  values.  24-29 
Good  agreement  was  found  in  the  region  between  400  and  650  nm, 24-27  but  not  in 
the  infrared. 28-29  Apparently  most  published  values  for  the  relative  variation 
of  air  afterglow  intensity  as  a  function  of  wavelength  are  erroneous  for 
X>900  nm.  Recent  observations  by  Bradburn  and  Lilenfeld30  over  this  same 
spectral  region  support  our  findings. 
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IV.  FITTING  PROCEDURE 


Individual  band  intensities  were  calculated  from  vibrational  level  popula¬ 
tions  determined  by  a  spectral  fitting  technique. 31  We  fit  regions  containing 
only  one  prominent  transition  from  the  level  of  interest,  i.e.,  fitting  single 
sequences.  The  spectral  fitting  code  determines  vibronic-level  number  densi¬ 
ties  within  an  electronic  state  and  corrects  for  spectral  overlap  between 
different  bands.  The  general  procedure  involves  calculating  synthetic  spectra 
whose  magnitudes  correspond  to  unit  population  for  the  bands  of  interest.  The 
populations  of  these  emitters  are  then  adjusted  in  a  linear  least-squares 
fitting  routine  to  reproduce  the  observed  spectrum. 

Spectral  simulation  begins  by  calculating  line  positions  and  their  rela¬ 
tive  intensities  following  the  methods  of  ttovacs.32  This  procedure  generates  a 
stick  spectrum  which  is  then  convoluted  with  a  triangular  slit  function  whose 
full  width  at  half  maximum  equals  the  monochromator  band  pass.  We  confirmed 
the  appropriateness  of  a  triangular  slit  function  by  scanning  the  monochromator 
over  the  atomic-line  spectrum  emanating  from  a  pen-ray  lamp,  and  then  measuring 
the  shapes  of  the  spectral  lines. 

We  used  the  spectroscopic  constants  of  Roux  et  al.33  to  calculate  band 
positions  and  found  them  to  agree  with  our  observations.  Synthetic  spectra 
generated  from  the  constants  tabulated  in  Lofthus  and  Krupenie^  could  not 
adequately  reproduce  the  wavelengths  of  the  higher  v'  levels  (even  after 
correcting  for  the  sign  error  on  (%ye  for  the  A  state  in  Lofthus  and  Krupenie). 

The  distribution  of  rotational  populations  was  taken  to  be  Boltzmann. 
Rotational  temperatures  for  the  microwave  discharge  spectra  were  chosen  to  give 
the  best  match  to  the  observed  rotational  contours  of  the  bands.  These 
rotational  temperatures  varied  between  450  and  1000  K  in  accord  with  our 
observations  on  microwave-discharge  lamps  using  other  techniques . 3^  The  other 
two  excitation  sources  produced  spectra  with  a  room  temperature,  300  K, 
rotational  distribution. 
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Spectra]  overlap  is  a  problem  in  all  band  sequence^.  Generally  bands  from 
higher  vibrational  levels  of  the  Av=n-1  sequence  overlap  with  the  lower 
vibrational  levels  in  the  Av=n  sequence.  This  problem  is  particularly  severe 
for  the  Av=+3,  +2,  +1,  and  0  sequences.  The  infrared  afterglow  system 
-  B-^IIg)  and  several  sets  of  atomic-nitrogen  lines  contaminate  the  microwave- 
discharge  spectra.  The  infrared  afterglow  system  also  affects  the  atom- 
recombination  data.  In  addition,  overlap  with  Av=n-1  first-positive  bands  is 
particularly  severe  in  these  data  due  to  the  larger  relative  populations  of  the 
high  v'  levels.  The  N2+  Meinel  bands  (A^ nu-X^ £g+)  are  prominent  in  the 
electron-impact  excitation  spectra.  Spectral  constants  of  Roux  et  al.-^  and 
Lofthus  and  Krupenie^  were  used  to  calculate  the  line  positions  of  the  B'-B  and 
A-X  systems,  respectively. 

The  extent  of  the  overlap  determined  the  fitting  strategy.  In  instances 
of  partial  overlap,  the  code  determined  the  populations  of  all  emitting  species 
directly.  With  complete  overlap  or  overlap  from  more  than  one  band  system, 
populations  of  the  contaminating  species  were  estimated  by  fitting  adjacent 
regions  of  the  spectrum  which  afforded  more  reliable  estimates.  These  popula¬ 
tions  were  then  used  to  generate  a  synthetic  spectrum  of  interfering  emissions 
in  the  region  of  interest  that  subsequently  was  subtracted  from  the  data. 
Fitting  the  residual  spectrum  then  determined  the  intensities  of  interest. 

Bands  for  which  the  relevant  feature  coincided  with  another  band  of  comparable 
or  greater  intensity,  or  which  were  severely  overlapped  with  a  significant 
band  for  which  no  independent  population  estimate  could  be  made,  were  not  used 
in  the  transition-moment  determination. 

Figures  1  and  2  illustrate  the  capabilities  of  the  spectral  fitting 
procedure.  They  show  successive  approximations  in  the  fitting  of  the  Av=0 
sequence  in  the  data  taken  on  LABCEDE.  The  data  are  displayed  as  the  light 
line,  the  fit  as  the  smoother,  darker  line.  In  Figure  1  only  the  v'-O  5 
levels  of  NyfB)  and  some  atomic  lines  were  included  in  the  fit.  Large 
discrepancies  are  evident  in  the  fit,  particularly  around  920  and  950  nm. 

These  regions  contain  the  1,0  and  2,1  Meinel  bands.  Including  the  Meinel  bands 
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in  the  fit  and  estimating  contributions  from  the  7,8  and  6,7  first-positive 
bands  gives  the  result  displayed  in  Figure  2.  The  entire  spectral  region  is 
fit  excellently,  lending  confidence  to  the  accuracy  of  the  intensity 
determinations . 

Situations  involving  subtraction  of  interfering  emissions  could  result  in 
«  systematic  errors.  The  particular  Einstein  coefficients  used  in  the  estimates 
determined  the  magnitudes  of  the  intensities  to  be  subtracted . 36  jn  regions 
of  the  spectrum  where  such  subtraction  has  the  potential  to  produce  signifi¬ 
cant  errors,  there  is  overlap  of  data  using  at  least  two  excitation  methods 
that  produce  spectra  with  markedly  different  characteristics.  Systematic 
errors  would  be  manifest,  therefore,  as  differing  trends  in  the  Re(r)  curves 
from  the  different  data.  All  systems  showed  similar  trends,  however,  indicat¬ 
ing  the  relatively  minor  role  of  such  systematic  errors. 

The  r-centroids  and  Franck-Condon  factors  necessary  for  the  evaluation  of 
the  reduced  intensities  and  the  construction  of  the  Re(?)  curve  were  calculated 
by  W.J.  Marinelli  and  T.  Quagliaroli  of  PSI  using  a  procedure  outlined  pre¬ 
viously.  37  Rydberg-Klein-Rees  potentials  were  first  calculated  for  both 
electronic  states  using  the  spectroscopic  constants  of  Roux  et  al.35  an(j  the 
approach  of  Tellinghuisen. 38  The  numerical  eigenfunctions  were  then  evaluated 
using  the  Numerov-Cooley  procedure-^  to  solve  the  radial  Schrodinger  equation. 
The  overlap  integrals  in  the  r-centroid  and  Franck-Condon-factor  calculations 
were  evaluated  using  Simpson's  Rule.  The  results,  which  cover  the  range  of 
v'=0-12  and  v"=0-24,  are  more  extensive  than  other  calculations , 3 > 35  but  agree 
»  quite  well  with  them  in  regions  of  overlap.  We  have  tabulated  them 
el sewhere . 
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V.  RESULTS 


The  reduced  intensities  were  averaged  for  all  spectra  taken  with  a  given 
excitation  method.  A  total  of  four  spectra  were  analyzed  for  the  discharge 
data,  three  for  the  atom-recombination  data,  and  two  for  the  data  taken  with 
the  LABCEDE  facility.  The  data  were  combined  to  form  the  relative  transition- 
moment  curve  by  minimizing  the  sum  of  the  squares  of  the  differences  of  the 
normalized  reduced  intensities  for  each  progression  from  the  analytical 
expression  which  best  represented  all  the  data.  In  practice,  we  determined 
the  oest  fit  of  a  given  analytical  form  to  the  reduced  intensities  of  one 
progression,  e.g.,  v'=2  or  v'=10,  and  then  varied  the  intensity  normalization 
factor  of  e<-  '.h  of  the  other  progressions  to  obtain  the  minimum  least-squares 
deviation  from  that  analytical  expression.  We  then  recomputed  the  best-fit 
analytical  expression  to  the  adjusted  progressions.  This  defined  a  new  line, 
and  the  individual  progressions  were  again  adjusted  to  match  this  line.  After 
several  iterations,  further  adjustment  of  the  progressions  did  not  signifi¬ 
cantly  improve  the  fit  between  all  of  the  data  and  the  analytical  form  which 
best  represented  them.  The  final  result  was  independent  of  the  initial 
progression  chosen  to  begin  the  fitting  procedure. 

We  tried  both  linear  and  quadratic  expressions  to  represent  Re  as  a  func¬ 
tion  of  rv,vn.  The  fitting  procedure  showed  the  quadratic  term  not  to  be 
statistically  significant.  A  few  trials  with  a  cubic  function  also  did  not 
improve  the  fit.  The  best  linear  fit  was 

Re(rV'V")  =  ±  0.04)  -  (0.654  +  0.030)  rvrv»  .  (9) 

Figure  3  shows  the  relationship  between  this  line  and  the  experimental  data. 

The  relative  transition-moment  curve  shown  in  Figure  3  was  placed  upon  ar. 
absolute  basis  by  multiplying  Eq.  (9)  by  a  normalization  factor,  C,  that  would 
bring  radiative  lifetimes  calculated  from  our  set  of  Einstein  coefficients 
into  congruence  with  the  experimental  values  of  Eyler  and  Pipkin. 20  That  is 
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We  found  a  value  of  C  =  (1.89  +  0.03)  x  10-2  provided  agreement  between  our 
calculations  and  all  of  Eyler  and  Pipkin's  experimental  lifetimes.  These 
.  lifetimes,  which  cover  the  range  v'=5-12,  should  be  the  most  accurate 

available  because  they  were  determined  via  laser-induced  fluorescence  in  a 
molecular  beam  under  collision-free  conditions.  We  note,  however,  that  this 

v 

set  of  lifetimes  appear  to  be  10  to  15  percent  shorter  than  those  determined 
in  the  most  reliable  of  the  other  experimental2 »41-44  or  theoretical 
studies. 3-5  our  selection  of  a  different  set  of  lifetimes  fcould  alter  our 
transition-moment  curve  by  five  to  seven  percent.  We  discuss  lifetimes 
further  in  the  next  section. 

The  normalization  factor,  £,  and  the  constants  in  Eq.  (4)  allow  us  to 
re-normalize  the  relative  transition-moment  variation  given  by  Eq.  (8)  so  that 
Re  can  be  expressed  in  Debye: 


Re(D)  =  (1.30+0.03)  -  (0.508+0.02)  r v.v»  .  (11) 

The  normalization  procedure  requires  summing  over  all  significant  transi¬ 
tions.  The  best  test  for  completeness  is  to  sum  the  Franck-Condon  factors 

from  a  given  v'  over  all  v".  The  resultant  sum  should  be  close  to  unity.  For 
20 

v'=0-10,  Z  q  ,  i,  >  0.9,  so  our  set  of  relative  Einstein  coefficients  is 
>  v"=0 

sufficiently  complete  for  accurate  normalization  in  Eq.  (10).  While  the  sum 
for  v'=ll,12  is  somewhat  smaller,  our  results  still  are  sufficiently  complete. 
Even  if  the  Franck-Condon  factor  for  transitions  to  v"=21  made  up  the  balance 
of  the  difference  between  the  Franck-Condon  factor  sum  up  to  v"=20  and  unity, 
the  term  in  Eq .  (4)  is  so  small  that  for  these  additional  transitions  our 
twenty-one  term  sum  over  the  relative  Av»vm  would  be  augmented  by  less  than 
ohe  percent. 
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In  practice  a  small  correction  to  the  sum  in  Eq.  (10)  should  be  made  to 
include  radiative  transitions  to  other  electronic  states,  in  this  case  W^Ay 
and  Werner  et  al.^  indicate  that  the  largest  of  these  corrections 

range  from  about  five  percent  for  v'=12  to  less  than  three  percent  for  v'=5. 
Their  incorporation  into  our  analysis  will  yield  a  normalization  factor,  C, 
which  is  the  same  within  experimental  error  as  the  value  we  have  used.  We 
have  chosen,  therefore,  not  to  make  these  corrections. 


Table  1  contains  a  complete  set  of  N2(B-A)  Einstein  coefficients. 


VI.  DISCUSSION 


A.  Transition-Moment  Functions 

The  present  results  agree  much  better  with  theoretical  calculations^-^  of 
the  transition-moment  function  than  with  the  results  of  previous  empirical 
investigations. 2,45-49  jn  an  cases  agreement  on  the  shapes  is  quite  good  over 
the  range  1.1  to  1.7  A.  Agreement  with  absolute  magnitudes  of  the  transition- 
moment  function  is  also  quite  good,  generally  within  10  to  15  percent.  Our 
measurements  appear  to  have  resolved  the  major  discrepancies  between  theory  and 
experiment . 

Figure  4  compares  our  experimentally  derived  variation  in  electronic 
transition  moment  with  internuclear  separation  with  that  determined  by 
Shemansky  and  Broadfoot  and  with  the  theoretical  calculation  of 
Werner  et  al.^  The  shape  of  our  curve  matches  that  of  Werner  et  al.  quite 
well.  Consequently,  our  results  agree  quite  well  with  theirs  on  the  relative 
variation  of  lifetimes  as  a  function  of  vibrational  level  and  on  the  radiative 
branching  ratios  from  a  common  upper  vibrational  level  to  the  various  lower 
levels.  In  contrast,  Shemansky  and  Broadfoot 's  transition-moment  function  has 
a  much  different  slope  which  results  in  a  much  different  variation  in  lifetimes 
with  vibrational  level  and  set  of  branching  ratios. 

Werner  et  al.  have  reviewed  the  various  experimental  and  theoretical 
transition-moment  functions.  All  of  the  previous  empirical  values^-49  diverge 
even  more  widely  from  our  results  than  those  of  Shemansky  and  Broadfoot. 2 
These  discrepancies  can  result  from  a  number  of  factors  depending  upon  the 
approach  taken.  For  example,  measurements  of  band  intensities  require  extreme 
care  in  determining  instrument  response  function,  in  separating  out  contribu¬ 
tions  to  the  apparent  band  intensity  from  overlapping  radiators,  and  in  avoid¬ 
ing  experimental  pitfalls  such  as  self  reversal  which  apparently  can  be  a 
problem  using  high  intensity  sources . 49-51 


477 


Accurate  response-function  determinations  require  uniformly  filling  the 
optical  system  with  the  calibration  lamp  output  as  well  as  with  the  radiation 
source  being  measured.  Monochromator  optics  usually  cannot  be  filled  by  plac¬ 
ing  a  tungsten  strip  lamp  in  front  of  the  entrance  slit.  Rather,  illumination 
from  the  strip  lamp  must  be  reflected  into  the  monochromator  from  a  diffuse 
white  source  such  as  a  BaSO/,  screen22»23  or  a  spectralon®  target. 

The  first-positive  system  is  subject  to  significant  overlap  from 
interfering  radiators.  This  interference  is  so  pervasive  that  accounting  for 
it  requires  spectral-fitting  capability.  Shemansky  and  Broadfoot  did  use 
spectral  fitting  in  their  analysis.  Other  groups  measuring  band  intensities 
did  not. 

Jeunehomme^  tried  to  infer  a  transition-moment  function  from  lifetime 
measurements  of  various  vibrational  levels.  Generally,  radiative  lifetimes 
cannot  be  determined  accurately  enough  for  this  procedure  to  prove  successful. 
In  even  the  best  circumstances,  radiative  lifetime  variations  with  v'  are 
fairly  insensitive  to  the  form  of  the  transition-moment  function.  This  is 
because  the  total  radiative  decay  rate  (the  inverse  of  the  radiative  lifetime) 
is  the  sum  of  a  number  of  components  which,  taken  together,  sample  a  wide  range 
of  internuclear  separations.  Comparing  to  experimental  variations  in  radiative 
lifetime,  however,  provides  a  good  cross  check  to  a  transition-moment  function 
derived  from  other  measurements  (vide  infrared). 

Several  groups  also  have  attempted  to  measure  transition  moments  by 
observing  the  radiation  emitted  by  high  temperature  gases  in  thermal  equilib¬ 
rium.  49,53-56  jn  this  approach,  the  gas  generally  is  heated  by  a  shock  wave 
or  high  current  arc.  These  measurements  often  are  plagued  by  interfering 
impurity  emissions,  and  generally  lead  to  overestimates  of  transition  probabil¬ 
ities.  Such  measurements  have  proven  to  have  uncertainties  of  factors  of  two 
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B.  Radiative  Lifetimes 

Figure  5  compares  lifetimes  calculated  from  the  present  results  with  those 
of  Shemansky  and  Broadfoot^  as  well  as  the  recent  measurements  of  Eyler  and 
Pipkin^O  and  the  theoretical  calculations  of  Werner  et  al.^  Our  transition- 
moment  function  reproduces  quite  well  the  change  in  lifetime  with  vibrational 
level  observed  by  Eyler  and  Pipkin.  The  discrepancy  between  our  results  and 
those  of  Shemansky,  on  the  other  hand,  is  considerable.  Although  somewhat 
different  in  magnitude,  our  relative  changes  in  lifetime  with  vibrational  level 
match  those  calculated  by  Werner  et  al. 

We  scaled  Werner  et  al.'s  results  to  Eyler  and  Pipkin's  lifetimes  by  first 
fitting  their  transition-moment  variation  to  a  quadratic  function.  This 
quadratic  function  was  used  as  a  relative  Re  function  to  calculate  relative 
transition  probabilities.  The  relative  probabilities  were  then  made  absolute 
by  applying  the  procedure  given  by  Eq.  (10).  This  resulted  in  a  nine  percent 
increase  in  Werner  et  al.'s  transition-moment  curve  at  every  point.  The  result 
of  this  procedure  is  a  set  of  lifetimes  which  agree  with  ours  to  better  than 
five  percent. 

Similar  treatment  of  the  theoretical  transition-moment  functions  of  Yeager 
and  McKoy^  and  Weiner  and  Ohrn^  result  in  even  smaller  increases  in  the 
magnitude  of  their  functions.  Rizzo  et  al.'s^  function,  on  the  other  hand, 
requires  a  12  percent  decrease  in  magnitude  for  it  to  come  into  congruence  with 
ours.  Their  transition-moment  function  results  in  a  set  of  radiative  lifetimes 
more  than  20  percent  shorter  than  those  measured  by  Eyler  and  Pipkin. 20 

A  number  of  groups  have  reported  experimental  values  for  radiative 
lifetimes  of  various  vibrational  levels  of  N2(B) . 20 , 41-44 , 57 , 58  Table  2 
summarizes  some  of  them.  In  most  instances  the  measurements  involved  fitting 
multiexponential  decays  at  a  number  of  pressures  and  extrapolating  one  of  the 
decay  components  to  zero  pressure.  This  approach  generally  is  hazardous  unless 
the  lifetimes  of  the  various  components  in  the  decay  differ  greatly.  The  N2(B) 


state  is  collisionally  coupled  to  various  levels  of  the  and  A^  Eu+ 

states. ^2, 44, 58-63  Coupling  effects  are  manifest  even  at  pressures  as  low  as  a 
mTorr.^2,44,63  This  coupling  typically  leads  one  to  overestimate  the  radiative 
lifetime.  Eyler  and  Pipkin's  measurements,  having  been  made  by  laser-induced- 
fluorescence  techniques  in  the  collision-free  environment  of  a  molecular  beam, 
should  be  free  from  these  effects  of  interstate  coupling. 

Our  results  indicate  lifetimes  for  vibrational  levels  0  and  1  N2(B),  that 
are  36  percent  and  24  percent  longer,  respectively,  than  those  given  by 
Shemansky  and  Broadfoot . *•  The  discrepancy  could  be  even  greater  since 
Shemansky  and  Broadfoot's  values  for  the  highest  vibrational  levels  are  10  to 
15  percent  longer  than  ours.  Scaling  their  results  down  accordingly  would 
lead  to  discrepancies  greater  than  50  percent  for  v'=0,l.  Ortiz  et  al.^7  and 
Heidner  et  al.'s5®  experimental  values  of  13  ±  1  us  and  10+2  us,  respec¬ 
tively,  for  the  radiative  lifetime  of  v'=0  appear  to  confirm  our  result  of 
12.1  us. 
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VII.  SUMMARY  AND  CONCLUSIONS 


Ve  have  determined  the  variation  in  the  transition  moment  with  inter- 
nuclear  separation  for  the  N2(B-A)  transition.  Our  approach  relies  on  the 
redundancy  of  intensity  measurements  from  several  different  sources  using 
several  different  detection  systems  to  reduce  possible  systematic  errors.  We 
have  made  extensive  use  of  spectral  fitting  to  reduce  uncertainties  in  band- 
intensity  determinations  caused  by  overlap  with  other  bands  or  band  systems. 

Our  results  agree  well  with  several  recent  theoretical  calculations,  but  show 
that  substantial  revision  of  currently  accepted  experimental  values  is 
necessary. 

Number  densities  of  N2(B,v)  determined  in  experiments  on  N2(B)  excitation 
by  various  sources  result  from  dividing  measured  band  intensities  by  the  appro¬ 
priate  Einstein  coefficient.  Our  results  indicate  that  using  the  Einstein 
coefficients  in  Lofthus  and  Krupenie  will  result  in  over  estimates  of  30  to 
60  percent  for  the  number  densities  in  v'=10-12  but  under  estimates  of  20  to 
30  percent  in  the  lowest  levels.  This  means  that  relative  number  densities  of 
vibrational  levels  0  and  12  could  differ  by  100  percent  or  more  using  our 
Einstein  coefficients  compared  to  previous  values. 
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Table  1.  Einstein  Coefficients  and  Wavelengths  for  ^(B-A)  Transitions 


Table  2.  ^(B^IIg)  Lifetimes  in  Microseconds 
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FIGURE  CAPTIONS 


Figure  1.  Av  =  0  Sequence  Data  (light  line)  and  Fit  (dark  line)  for  N2(B-A) 
and  Atomic  Lines 

Figure  2.  Av  =  0  Sequence  Data  and  Fit  for  N2(B-A),  N2+(A-X),  and  Atomic  Lines 
with  Estimated  Contributions  from  Other  N2(B-A)  Sequences  Subtracted 

Figure  3.  Experimentally  Determined  Variation  in  the  Electronic-Transi t ion 
Moment  with  r-centroid  for  the  N2  First-Positive  System 

Figure  4.  Variation  in  Electronic  Transition  Moment  with  Internuclear  Distance 
for  the  N2(B-A)  Transition 

"igure  5.  Variation  in  ^(B)  Radiative  Lifetimes  as  a  Function  of  Vibration 
Level 
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Figure 
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Figure  4 
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1 .  INTRODUCTION 


1 . 1  An  Overview  of  the  COCHISE  Apparatus 

The  COCHISE  facility,  situated  at  the  Air  Force  Geophysics  Laboratory, 
is,  in  the  simplest  terms,  a  relatively  large,  cryogenically  cooled  test 
chamber  which  may  be  used  to  study  the  state  of  the  reaction  products  created 
as  the  result  of  a  variety  of  possible  gas  phase  reactions  such  as  chemi- 
excitation  processes, 

A  +  BC  -*■  AB*  +  C  (  1  ) 

'  A  +  B  +  AB*  (2) 

or  vibrational  exchange  or  electronic-state  quenching  reactions 

A*  +  B  -*■  A  +  B*  ,  (3) 

where  the  quantities  A,  B,  C  are  meant  to  represent  atoms  or  molecules  and  the 
asterisk  implies  electronically  or  vibrationally  excited  states. 

The  predominant  IR  active  mode  of  excitation  of  AB*  is  vibrational;  how¬ 
ever,  some  atmospheric  species  (most  notably  N2  and  O2)  possess  electronic 
states  which  can  radiate  in  the  IR.  It  is  often  difficult  to  deduce  detailed 
information  about  these  processes  from  in  situ  atmospheric  measurements  of  IR 
emission;  such  data  often  have  low  spectral  resolution  and  contain  contribu¬ 
tions  from  competing  thermal  and  scattering  mechanisms.  However,  laboratory 
investigations  of  individual  chemical  processes  under  carefully  controlled 
conditions  can  provide  detailed  mechanistic  and  spectroscopic  information; 
this  information  can,  in  turn,  be  used  in  conjunction  with  atmospheric  models 
to  both  interpret  in  situ  data  and  define  further  measurements. 

Infrared  emission  spectroscopic  experiments  designed  to  investigate  gas 
phase  chemiexcitation  reactions  require:  1)  high  sensitivity  to  the  rela¬ 
tively  weak  molecular  radiation  from  the  products  of  reactions  (1)  and  (2); 
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2)  adequate  resolution  to  allow  conclusive  spectral  analysis;  and  3)  specif¬ 
icity  to  the  chemiexcitation  process  under  investigation,  i.e.,  negligible 
interference  from  radiative,  collisional,  and  surface  relaxation  processes. 
These  conditions  work  against  each  other;  the  achievement  of  condition  (3) 
necessitates  the  use  of  small  pressures  and/cr  product  yields,  which  results 
in  low  radiation  levels  and  increases  the  difficulty  of  achieving  condi¬ 
tions  (1)  and  (2).  Furthermore,  the  limitation  of  IR  detection  sensitivity  by 
thermal  background  radiation  in  the  field  of  view  severely  inhibits  spectrally 
resolved  observations  of  long-wavelength  radiation  from  important  atmospheric 
species  such  as  CO2  (15  pm)  and  O3  (10  pm).  Previous  IR  chemiluminescence 
methods,  such  as  arrested  relaxation1-2  and  discharge  flow4'5  coupled  with  the 
use  of  circular-variable-filter  spectrometers, 5  grating  monochroma  tors , 2  and 
Fourier  transform  spectrometers, 1-4  have  been  highly  successful  in  studies  of 
chemiexcitation  at  wavelengths  less  than  -5  pm.  In  general,  however,  these 
methods  suffer  from  thermal  background  limitations;  for  example,  at  10  pm,  the 
spectral  radiance  from  a  typical  AB  species  at  a  column  density  of  ~10°  cm 
(often  representative  of  such  experiments)  would  be  ~9  orders  of  magnitude 
less  than  that  from  a  room  temperature  blackbody. 

Infrared  detectivity  can  be  dramatically  improved  by  resorting  to  a  cryo¬ 
genic  background,  as  shown  in  Figure  1.  Since  the  limiting  noise  mechanism  is 
the  statistical  fluctuation  in  the  photon  flux  associated  with  the  detector's 
environment,  the  detectivity  of  the  system  is  inversely  proportional  to  the 
square  root  of  that  flux.  The  results  plotted  in  Figure  1  are  normalized  to 
the  detectivity  at  293  K.  Although  considerable  enhancement  in  detectivity  is 
obtained  with  an  80  K  background,  the  effects  of  background  radiation  can  be 
virtually  eliminated  by  cooling  the  entire  reaction  chamber/detection  system 
to  temperatures  below  40  K.  As  a  further  benefit  of  cryogenic  temperatures, 
rapid  cryopumping  of  reagent  gases  can  be  attained  near  20  K,  and  radiation 
leakage  from  the  external  vacuum  system  can  then  be  eliminated  by  isolating 
the  reaction/detection  system  during  the  experiments. 

COCHISE  is  designed  for  partial  simulations  and  detailed  experimental 
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COLD  BACKGROUND  DETECTIVITY  IMPROVEMENT 


A-3322 


Figure  1.  Reduction  in  Background  Photon  Flux  by  Cyrogenic  Operation 
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studies  of  high-altitude  IR  excitation  phenomena,  particularly  chemiexcitation 
and  excitation  transfer  processes.  The  primary  design  objective  for  the 
apparatus  was  to  eliminate  the  limitations  to  IR  detectivity  outlined  above. 
The  experimental  direction  taken  in  the  development  of  Cold  Chemiexcitation 
Infrared  Simulation  Experiment  (COCHISE)  has  been  to  reduce  the  temperatures 
of  both  the  reaction  chamber  and  the  detection  system  to  ~20  K  to  optimize  the 
benefits  of  cold  background  and  cryopumping.  A  sta te-of-the-art  IR  detector 
and  a  cryogenic  scanning  monochromator  are  coupled  to  a  long-path  reaction 
cell  to  allow  spectrally  resolved  detection  of  IR  radiation  (2  to  20  ym)  from 
emitting  species  at  concentrations  below  106  molecules  cm--*.  The  techniques 
of  discharge  flow  and  arrested  relaxation  are  combined  to  allow  mixing  of  the 
reacting  species  in  the  field  of  view  at  steady-state  pressures  near  3  mTorr, 
so  that,  in  most  cases,  the  chemiexcited  products  can  be  formed  and  observed 
under  nearly  single-collision  conditions.  The  high-speed  cryopumping  of  the 
reaction  vessel  removes  the  excited  species  from  the  field  of  view  of  the 
detection  system  before  radiative  or  collisional  relaxation  can  occur  and 
eliminates  contamination  by  back -diffusion  from  the  chamber  walls. 

A  schematic  of  the  entire  COCHISE  apparatus  is  given  in  Figure  2.  The 
reaction  cell  and  detection  system  are  enclosed  within  a  cryogenic  thermal 
shroud,  which  is,  in  turn,  enclosed  by  a  main  vacuum  chamber  to  minimize  con¬ 
ductive  losses.  All  surfaces  within  the  shroud,  with  the  exception  of 
selected  optical  components  and  reagant  gas  lines,  are  held  at  temperatures  as 
low  as  -20  K;  this  allows  rapid  cryopumping  of  all  reagent  and  background 
gases  except  He  and  H2*  The  reagent  gases  enter  the  reaction  cell  through 
four  sets  of  opposing  jets  (see  also  Figure  3)  and  interact  along  the  center 
line  of  the  cell;  the  resulting  radiation  is  viewed  along  the  center  line 
through  a  lens  by  a  grating  monochromator  detection  system.  The  data  process¬ 
ing,  temperature  control,  and  system  housekeeping  are  performed  by  an  external 
computer.  The  various  components  of  the  system  are  described  in  more  detail 
in  later  sections. 

The  COCHISE  facility  as  described  above  thus  provides  a  unique  capability 
for  the  study  of  the  fundamental  kinetic  paths  and  product  distributions 
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Figure  2.  Schematic  of  COCHISE  Facility 
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Figure  3.  Scale  Drawing  of  Chemiluminescence  Reaction  Cell.  Ar/N2  or 
Ar/02  mixtures  are  excited  by  microwave  discharges,  D,  and 
mix  with  the  counterflow  of  02  from  nozzles  N  in  an  axisym- 
metric  reaction  zone  at  the  chamber  center.  Baffles,  B, 
restrict  detection  system  f ield-of-view  to  the  on-axis 
reaction  volume.  M  is  the  plane  mirror  which  increases  the 
intensity  of  the  radiation  collected  by  lens  L.  Spectral 
response  calibration  is  performed  using  blackbody  source,  C, 
embedded  behind  a  small  hole  in  the  end  mirror. 
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produced  by  chemical  reactions  between  various  excited  species,  free  radicals, 
and  neutral  molecules.  Such  studies,  although  of  basic  interest  in  their  own 
right,  can  prove  to  be  invaluable  in  the  understanding  of  radiative/energy 
transfer  phenomenology  occurring  at  rarefied  gas  densities  such  as  auroral  or 
nuclear  phenomena  in  the  upper  atmosphere.  To  the  best  of  our  knowledge,  the 
COCHISE  device  is  the  only  facility  in  the  nation  which  can  provide  such 
detailed  information  on  these  fundamental  chemical  processes. 

1 . 2  Applications  and  Accomplishments 

The  COCHISE  facility  has  been  used  to  develop  an  extensive  data  base  on 
the  spectroscopy  of  infrared  (2  to  16  ym)  fluorescence  from  the  effluents  and 
interactions  of  microwave-discharged  N2/Ar  and  02/Ar  mixtures.7  Some  results 
of  these  measurements  include  the  formation  of  NO(v)  from  the  interaction  of 
discharged  N2/Ar  and  02,®  the  formation  of  03(v)  from  discharged  02/Ar,9>10 
and  the  excitation  of  high- lying  Rydberg  states  of  Ar  in  the  discharge.11 
Studies  currently  under  way  include  the  interaction  of  discharged  N2/Ar  with 
CO  and  the  potential  for  formation  of  N20(v)  from  various  processes.12 

Discharged  N2/Ar  mixtures  provide  a  wealth  of  energy-rich  states  of 
nitrogen.  The  most  prominent  effluents  of  such  discharges  are  the  atomic 
species  N(4S,  2D,  2P) ,  the  principle  molecular  metastable  N2(A3EU+),  and 
vibrationally  excited  N2(v).  Other  excited  triplet  and  singlet  forms  of  N2, 
e.g.,  B3rig  and  a^g,  have  relatively  short  radiative  lifetimes  (1  to  100  ps) 
and  thus  may  be  expected  to  disappear  soon  after  exiting  the  active  discharge 
if  there  are  no  other  source  terms  besides  direct  excitation.  These  excited 
species  can  give  rise  to  vibrationally  excited  NO  and  CO  through  exoergic 
reaction  or  energy  transfer  with  atmospheric  species  such  as  02,  CO,  and  C02. 
Our  previous  and  current  experiments  have  examined  the  interaction  of 
discharged  nitrogen,  N2* ,  with  02  and  CO,  which  occurs  primarily  through  the 
N(2D)  metastable. 12  COCHISE  will  further  examine  the  potential  effects  of 
N(2p)  in  these  systems,  as  well  as  to  more  completely  characterize  the  N2*-co 
interactions.  Our  present  studies  of  the  N2*-C0  system  indicate  at  least 
three  different  excitation  processes  responsible  for  emission  in  the  4  to  6  ym 
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region.  One  reaction,  which  results  in  excitation  of  CO(v),  is  probably 
energy  transfer  from  N(2d): 

N(2D)  +  CO  -*•  CO  (v  <  9)  +  N(4S) 

A  second,  unidentified  process,  which  is  most  apparent  at  low  CO  concentra¬ 
tions,  appears  to  excite  CO(v  <  23)  with  extremely  low  (<  10  K)  rotational 
temperatures.  However,  the  COCHISE  data  base  at  this  time  is  still  too  sparse 
to  permit  a  useful  characterization  of  this  process.  This  extraordinary  com¬ 
bination  of  high  vibrational  and  low  rotational  excitation  is  intriguing,  and 
may  be  related  to  gas-surface  interactions  at  the  wall  of  the  cryogenic  reac¬ 
tion  chamber.  We  propose  to  obtain  a  more  complete  data  base  as  a  function  of 

CO  concentration  and  surface  aging.  An  additional  manifestation  of  the  n2*-co 
interaction  is  an  unidentified  band  system  occurring  between  4.3  to  4.5  pm. 

The  spectroscopic  parameters  for  this  band  system  are  remarkably  similar  to 

those  for  ground-state  N2(v  =  1-6),  which  is  not  supposed  to  radiate  in  the 
gas  phase.  Thus,  if  this  assignment  is  correct,  an  excitation  process  on  the 
cryogenic  surface  is  indicated.  We  propose  to  obtain  more  data  for  a  variety 
of  N2  and  CO  conditions,  and  to  perform  a  conclusive  spectroscopic  analysis  of 
these  bands.  In  addition,  we  propose  to  examine  the  effects  of  N(2p)  via  the 
reactions 


n(2p)  +  co  -*•  n(2d)  +  co(v  <;  4) 

+  N(4S)  +  CO(v  <  15) 

While  N(^P)  cannot  be  isolated  from  N(^D),  the  ratio  of  the  two  species  can  be 
varied  to  some  extent  by  adjusting  the  N2  mole  fraction  and  discharge  power, 
so  a  parametric  study  of  C0(v)  fluorescence  will  indicate  whether  N(2P) 
excites  C0(v).  We  are  planning  to  extend  this  type  of  investigation  to  other 
atmospheric  collision  partners,  e.g.  C02  and  N20.  Active  nitrogen  is  known  to 
be  highly  reactive  with  these  species,12  and  we  will  examine  the  IR  fluores¬ 
cence  spectrum  for  emission  from  vibra tionally  excited  C02  (4.3,  15  um),  N20 

(4.5,  7.8,  17  ym),  CO  (4.7  um),  and  NO  (5.3  um).  These  experiments  will 


512 


require  the  development  of  techniques  for  heating  the  reagent  lines  suffi¬ 
ciently  to  prevent  condensation  of  C02  and  N20;  this  capability  was  designed 
into  COCHISE,  but  has  never  been  implemented. 

Recent  COCHISE  measurements  have  also  provided  new  spectral  data  on  N2 
electronic  emission  in  the  short-wavelength  infrared  (SWIR) ,  seemingly  excited 
in  the  microwave  discharge.®  The  W®AU  -*■  B%g  transitions  are  clearly  indi¬ 
cated  between  2  and  5  pm,1'*  some  B®IIg  *  A®EU+  transitions  appear  to  be 
present,  and  several  other  features  are  unidentified  as  yet.  We  will  pursue 
these  studies,  using  a  combination  of  parametric  experiments  with  varying  N2 
flows  and  counterflow  gases  and  a  spectral  fitting  analysis,  to  characterize 
the  excitation  and  quenching  of  this  emission  and  its  spectroscopic  distribu¬ 
tions  and  branching  ratios  in  the  infrared. 

In  addition  to  N2*,  metastable  and  vibrationa llv  excited  forms  of  02 , 

02*,  are  potential  precursors  to  atmospheric  IR  phenomena,  but  are  relatively 
poorly  understood  in  terms  of  kinetics  and  reaction  paths.1®  COCHISE 
experiments  indirectly  point  to  excitation  of  03(v)  by  energy  transfer  from 
02(b  ^  )  and/or  02(v),  and  auroral  enhancements  in  03(v3)  and  C02(v2) 
radiation  as  observed  by  HIRIS  may  be  related  to  similar  processes.16*18  We 
are  planning  to  examine  02*  interactions  with  the  same  reactive  species  used 
for  N2\  namely  CO,  C02  and  N20,  to  survey  the  potential  of  various  forms  of 
02*  as  precursors  to  infrared  radiation.  02*  will  be  generated  by  microwave 
discharge  of  02/Ar.  The  02  states  of  interest  include  02(a1Ag,  b^g*,  c1^-, 
A3EU+,  A'®Au)  and  02(v).  Spectral  surveys  will  cover  the  range  2  to  17  pm, 
with  emphasis  on  potential  IR  radiators  such  as  NO,  CO,  N20,  C02  and  03.  This 
exploratory  work  will  be  used  to  identify  key  02*  states  and  reactions  for 
further  study. 

The  potential  for  chemical  generation  of  N20  in  the  upper  atmosphere  is 
a  subject  of  continuing  debate.  However,  reactions  of  atmospheric  species 
which  can  result  in  formation  of  N20  at  low  pressures  have  yet  to  be 
conclusively  identified.  The  possible  reactions  would  require  at  least  one 
excited-state  reactant  species,  and  thus  may  occur  most  prominently  in 


laboratory  discharges  or  strong  aurorae.  Furthermore,  these  reactions  are 
sufficiently  exoergic  to  produce  vibrationally  excited  N20,  resulting  in 
infrared  chemiluminescence  near  4.5  {V3),  7.8  (v^),  and  17  (V2)  ym.  We  are 
using  the  COCHISE  facility  to  investigate  these  phenomena,  and  have  observed 
emission  from  N20(v3)  formed  in  microwave  discharges  of  flowing,  dilute 
mixtures  of  N2/02/Ar  at  ~1  Torr.^  Other  measurements  of  the  interactions  of 
discharged  02/Ar  with  N2  and  of  discharged  N2/Ar  with  02  give  no  measureable 
formation  of  N20(v).  These  results  contraindicate  the  formation  of  N20(v) 
from  the  reactions  of  02(v)  with  N2,  metastable  02  with  N2,  and  metastable  N2 
with  02.  We  suggest  that  excited  and/or  metastable  forms  of  both  N2  and  02 
are  required  in  the  critical  reaction(s).  This  could  entail  reactions  of  one 
of  the  forms  of  N2*  with  02*  or  perhaps  C>3(v).  We  propose  to  pursue  this 
issue  in  further  COCHISE  experiments.  The  first  step  in  this  investigation  is 
to  examine  the  IR  spectra  of  discharged  N2/02/Ar  mixtures  over  a  wide  range  of 
N2  and  02.  partial  pressures,  temperature  and  total  pressure,  to  ascertain 
under  what  conditions  N20(v)  formation  in  the  discharge  is  optimized.  These 
results  should  give  some  clues  as  to  the  identities  of  the  key  N20  precursors, 
since  we  know  from  separate  experiments  how  the  production  of  O2* ,  N2* ,  and 
O3 ( v )  varies  with  these  conditions.  Such  measurements  will,  of  course, 
provide  similar  information  on  NO,  O3,  and  N02  formation  which  would  be 
applicable  to  discharged  air. 

A  further  and  necessary  development  of  this  line  of  research  will  require 
the  modification  of  COCHISE  to  permit  discharge  excitation  of  the  counterflow 
gas.  Owing  to  the  difficulty  of  dissipating  heat  from  high-power  (~50W) 
microwave  units,  this  modification  would  probably  take  the  form  of  a  lower 
powered,  DC  discharge.  Such  a  modification  permits  the  study  of  _in  situ 
fluorescence  resulting  from  the  direct  interactions  of  two  or  more  different 
excited  species.  Example  reactions  which  could  be  studied  in  this  configu¬ 
ration  include  N2(A)  +  0  (to  look  for  NO(v)  formation),  N2(A)  +  C>3(v), 

N2(A)  +  02(1A,  ) ,  etc.  Such  reactions  are  little-known  but  could  play  an 

important  role  in  strong  electron  bombardments  where  metastable  production 
rates  are  large.  Thus  this  proposed  modification  to  COCHISE,  while  rather 
extensive,  would  significantly  increase  the  power  of  the  facility  to  study 
high-energy  atmospheric  phenomena. 
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The  observation  in  COCHISE  of  IR  emission  from  high-lying  Ar  Rydberg 
states  excited  in  the  microwave  discharges  has  opened  the  possibility  for 
performing  similar  studies  of  0(1)  and  N(I)  radiation.  Based  on  the  results 
of  the  Ar(I)  experiments,  these  measurements  would  use  discharged  02/Ar  and 
N2/Ar  mixtures  at  ~0. 1  Torr,  with  the  highest  possible  diatomic  partial 
pressure.  The  Ar  is  required  for  the  discharge  to  persist,  and  its  spectral 
interference  will  have  to  be  removed  by  subtracting  background  Ar  spectra. 
Also,  it  will  probably  prove  advantageous  to  reconfigure  the  field  of  view 
such  that  the  discharge  tubes  are  more  directly  viewed  by  the  detection 
system.  This  can  be  done  by  a  combination  of  re-orienting  the  reaction  cell 
and  introducing  reflective  surfaces  into  the  chamber.  These  measurements 
would  provide  information  on  the  spectroscopy  and  excitation  mechanisms  for 
these  states  in  microwave  discharges,  and  the  resulting  data  base  could  be 
linked  to  other  on-going  plasma  experiments  at  AFGL. 


2.  THE  COCHISE  APPARATUS 


2.1  The  Chamber 


A  schematic  of  the  COCHISE  apparatus  is  shown  in  Figure  2.  The  cryogenic 
portion  of  the  apparatus  is  thermally  isolated  inside  a  cylindrical  main 
vacuum  chamber  ~3m  long  and  ~1.5m  in  diameter  (Figure  2).  The  thermally 
insulating  vacuum  enclosure  is  maintained  at  ~10~10  atm  by  a  50-cm  (20-in.) 
cold-baffled  diffusion  pump  (Consolidated  Vacuum  Corp. )  backed  by  a  large  two- 
stage  mechanical  pump  (Leybold-Heraeus  DK-200) .  The  light-tight  aluminum 
shroud,  ~2.5m  long  and  1.25m  in  diameter,  is  supported  on  stainless  steel  and 
phenolic  standoffs  and  is  completely  surrounded  by  a  multilayer  thermal  shield 
to  provide  insulation  from  the  300  K  thermal  radiation  of  the  outer  wall. 

This  shield  consists  of  19  alternating  layers  of  aluminized  Mylar,  Dacron 
bridal  veil,  and  0.76  mm  (0.030  in.)  Scott  Industrial  Foam  and  permits  a  total 
radiative  load  of  <  50W  from  the  300  K  outer  wall  to  the  20  K  shroud.  This 
corresponds  to  a  net  reduction  of  over  2  orders  of  magnitude  in  the  radiative 
heat  flux,  allowing  a  hugh  increase  in  the  efficiency  and  load-carrying  capac¬ 
ity  of  the  cooling  system.  The  shroud,  reaction  cell,  and  detection  system 
are  cooled  to  ~20  K  by  a  6Q0W  closed-cycle  helium  refrigeration  system  (Cryo¬ 
genic  Technology,  Inc.).  Extended-stem  valves  permit  control  of  coolant  flow 
in  three  independent  loops  supplying  the  shroud,  reaction  cell,  and  microwave 
discharge  heat  sink  (see  below). 

Component  temperatures  are  monitored  at  roughly  25  locations  by  cali¬ 
brated  platinum  resistance  thermometers.  Direct  current  voltage  is  supplied 
from  an  external  source;  the  resulting  voltage  drop  across  each  sensor  is 
amplified,  processed  by  analog-digital  conversion  in  a  PDP-8E  (Digital  Equip¬ 
ment  Corp.)  computer,  converted  to  absolute  temperature,  and  displayed  on 
various  storage  oscilloscopes  or  hard  copy  devices.  Temperatures  of  gas  flow 
lines  and  optical  components  are  controlled  in  13  different  areas  by  resistive 
heaters  mounted  on  the  appropriate  hardware.  Temperature-controlled  portions 
are  isolated  from  surrounding  areas  by  stainless  steel  or  phenolic  standoffs. 
Power  to  the  heaters  is  percentage-controlled  by  the  computer  via  solid-state 
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switches  so  that  the  temperature  of  each  heat  station  is  independently  sampled' 
and  regulated  to  within  ±0.5  K. 

Because  the  internal  cryogenic  chambers  are  rapidly  cryopumped,  they  may 
be  vacuum-isolated  from  each  other  and  from  the  external  vacuum  system  during 
the  measurements  to  seal  off  external  light  leaks.  This  is  accomplished  by 
means  of  specially  modified,  cryogenic,  latching  solenoid  valves  of  a  type  not 
commercially  available.  These  valves,  which  are  operated  from  an  external 
control  panel,  can  be  opened  and  closed  by  separate  solenoids  and  will  remain 
in  either  position  in  the  absence  of  applied  power;  microswitches  mounted  on 
the  valve  shafts  provide  positive  readouts  of  the  valve  positions. 

The  physical  arrangement  shown  in  Figure  2  emphasizes  flexibility  in  use 
of  the  apparatus.  The  main  vacuum  chamber  and  cryogenic  shroud  have  full¬ 
opening  end  doors  and  a  number  of  access  ports  through  which  all  the  services 
for  the  experiments,  e.g.,  cryogenic  coolant,  microwave  power,  reagent  gases, 
and  instrumentation,  are  introduced.  All  diagnostic  and  control  devices  are 
mounted  externally.  The  ready  availability  of  extra  access  ports  and  vacuum 
feedthroughs  allows  a  high  degree  of  adaptability  to  future-diagnostic  con¬ 
cepts,  such  as  use  of  lasers  or  atomic  resonance  line  sources.  Details  of  the 
chamber  subsystems,  i.e.,  the  reaction  chamber,  refrigerator  system,  etc., 
will  be  presented  in  later  sections. 

2 . 2  Reaction  Cell  and  Reagent  Gas  Lines 

A  schematic  of  the  chemiluminescence  reaction  chamber  is  given  in 
Figure  3.  The  copper  reaction  cell  is  cylindrically  symmetric,  60  cm  in 
length,  40  cm  in  diameter.  To  control  the  cryopumping  speed  ^and  thus  the 
reagent  gas  residence  time  -  see  below),  the  reaction  cell  temperature  can  be 
independently  varied;  temperature  uniformity  to  ±0.1  K  over  the  entire  surface 
is  achieved  through  proper  positioning  of  coolant  flow  lines  and  heater 
elements. 
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The  reagent  gases  are  metered  into  the  system  by  motorized  Granville- 
Phillips  leak  valves;  flow  rates  are  monitored  by  calibrated  mass  flowmeters 
(Brooks).  One  gas  mixture  flows  through  four  sapphire  microwave  discharge 
sidearms  on  one  side  of  the  cell;  another  gas  mixture  flows  through  four 
counterflow  sidearms  on  the  opposite  side  of  the  cell.  Up  to  three  gases  can 
be  mixed  prior  to  introduction  into  the  discharge  cavities;  two  gases  can  be 
separately  mixed  for  the  counterflow.  Upon  penetration  of  the  outer  vacuum 
chamber,  the  room- temperature  gases  undergo  a  two-stage  heat  exchange  process. 
The  first  stage,  outside  the  shroud,  is  a  copper  plate  cooled  with  liquid 
nitrogen;  the  second  stage,  inside  the  cryogenic  shroud,  permits  variable 
temperature  selection  over  the  30  to  500  K  range.  The  gases  then  pass  through 
temperature-controlled  copper  lines  to  the  appropriate  inlet  sidearms  and  into 
the  reaction  volume.  Just  prior  to  entering  the  reaction  cell,  gases  in  the 
sapphire  discharge  tubes  are  ionized  and  excited  by  four  modified  McCarroll- 
Evenson7  microwave  discharge  cavities  (Opthos  Instruments)  powered  by  Raytheon 
PGM  10  power  supplies  (2450  MHz,  100W)  with  a  variable  duty  cycle.  Typical 
operation  is  at  50  to  90  percent  of  maximum  power  with  duty  cycles  of  either 
100  or  50  percent  at  23  Hz,  creating  steady-state  conditions  while  the  dis¬ 
charges  are  on;  all  four  discharges  are  driven  by  a  single  pulse  generator  to 
insure  synchronous  operation  to  within  a  few  microseconds.  Excess  heat  gener¬ 
ated  by  the  discharge  plasmas  is  transmitted  to  a  copper  heat  sink  plate  where 
it  is  removed  by  the  helium  coolant;  the  heat  transfer  rate  is  controlled  via 
an  extended-stem  valve  on  the  coolant  as  described  above.  In  normal  oper¬ 
ation,  the  discharges  are  held  at  the  desired  gas  temperature  (typically 
80  to  100  K  in  experiments  to  date)  to  within  1  or  2  K. 

As  shown  in  Figure  3,  the  reagent  gases  are  expanded  into  the  interaction 
volume,  with  equal  mass  flow  rates,  through  four  sets  of  diametrically  opposed 
inlet  jets  which  are  equally  spaced  along  the  cylindrical  cell.  Flow  condi¬ 
tions  in  the  reaction  cell  have  been  modeled®  for  a  range  of  experimental 
conditions;  free  expansion  occurs  in  the  reaction  cell,  and  the  gas  reaches  a 
limiting  velocity  in  which  nearly  all  the  thermal  energy  in  the  inlet  tube  is 
converted  into  kinetic  energy.  Near  the  cylinder  axis,  an  axisymme tr ic  stag¬ 
nation  point  occurs,  and  the  thermal  energy  of  the  gas  molecules  returns  to 
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that  of  the  gas  in  the  inlet  tubes.  Gas  density  variations  along  the  cylinder 
axis  (reaction  volume)  are  predicted  to  be  <  10  percent.  The  pressures  in  the 
reaction  zone  and  in  the  sidearms  are  monitored  by  MKS  Baratron  capacitance 
manome ters. 

At  typical  flow  rates  used  in  the  past,  ~4  standard  liter/min  of  various 
Ar/N2/02  combinations  on  each  side,  the  resulting  pressures  are  ~1  Torr  in 
each  sidearm  and  ~3  mTorr  in  the  reaction  zone;  the  corresponding  gas  resi¬ 
dence  times  are  ~2  msec  from  the  discharge  to  the  expansion  region,  ~0.5  msec 
from  the  expansion  region  to  the  field  of  view,  and  ~0.3  msec  in  the  field  of 

view. ^ 


The  chemiluminescent  reactions  occur  primarily  in  the  stagnation  region 
where  the  opposing  flows  meet.  The  IR  radiation  is  viewed  by  the  detection 
system  through  a  7.6  cm  (3  in.)  f/5  antireflection  (AR)-coated  germanium  lens; 
the  resulting  observed  volume  is  roughly  cylindrical  along  the  optical  axis  so 
that  only  radiation  occurring  within  ~4  cm  of  the  axis  is  observed.  A 
polished  gold  mirror  at  the  far  end  of  the  ce1  enhances  the  collection 
efficiency  of  the  system.  Internal  baffles  on  the  mirror  and  lens  discrimi¬ 
nate  against  scattered  light  from  the  discharges.  The  walls  of  the  reaction 
cell  are  not  in  the  field  of  view.  The  spectral  response  of  the  optical 
system  is  routinely  calibrated  with  a  variable- temperature  blackbody  source 
imbedded  in  the  center  of  the  mirror. 

2.2.1  Gas  Manifold  and  Flow  Calibration 


A  diagram  of  the  gas  manifold  is  shown  in  Figure  4.  The  solenoid  valves 
are  located  on  the  high  pressure  side  of  the  lines  and  the  controlling 
switches  are  on  the  COCHISE  control  panel.  There  are  two  solenoid  valves,  one 
for  discharge  gases  and  another  for  counterflow.  Similarly,  there  are  two 
pump  out  valves;  also  one  each  for  discharge  and  counterflow.  Figure  4  shows 
the  gas  manifold  for  each  gas  line.  There  are  currently  three  lines  on  the 
discharge  side  (Ar,  N2,  and/or  02 )  and  two  on  the  counterflow  side  (Ar,  N2  or 
C>2  or  CO) . 
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Figure  4.  Gas  Manifold 


In  order  to  pump  out  the  lines,  the  pump-out  valve  must  be  engaged  and 
the  Whitey  and  board  valves  open  as  well.  This  is  the  mode  in  which  the 
discharge  and  counterflow  lines  are  leak-checked.  A  rate  of  rise  in  the  TC 
readout  is  monitored  with  the  pump-out  valve  closed.  Typical  values  are 
<  0.1  micron/min.  The  lines  may  be  pumped  out  up  to  the  pressure  gauges  on 
the  gas  cylinders  by  opening  the  solenoid  valves.  This  is  typically  done  each 
morning  prior  to  beginning  experiments. 

To  operate  the  gas  manifold  system,  the  following  should  be  done  in 
order:  1)  pump  out  the  lines  up  to  the  cylinder  regulator  ,  2)  close  the 
bypass  valve,  3)  set  up  the  desired  flows  by  adjusting  the  leak  control 
switches  on  the  COCHISE  control  panel  (be  sure  to  turn  on  the  flow  meters) , 
and  4)  close  the  pump-out  valves  and  open  the  reaction  cell  valves  (there  are 
a  total  of  eight,  four  each,  for  counterflow  and  discharge). 

Typical  total  mass  flow  rates  are  0.11  g/s  from  each  of  the  counterflow 
and  discharge.  This  corresponds  to  flow  rates  of  ~4  standard  liters/min  from 
each.  It  is  important  for  mass  balance  to  be  maintained  so  whenever  one 
component  of  the  discharge  counterflow  is  changed  the  other  component  must 
also  be  adjusted.  Figures  5  through  14  show  the  calibrations  for  the  gases 
Ar,  C>2,  N2(CO)  for  each  of  the  flow  meters.  The  units  in  Figures  5  through  8 
are  g/s  while  those  of  the  other  curves  are  in  umoles/s  or  cm3/min. 

2.2.2  Field  of  View 


In  order  to  discriminate  against  surface  emissions,  it  is  critical  for 
the  field  of  view  to  underfill  the  mirror.  The  lens  and  position  of  the 
monochromator  were  chosen  to  accomplish  this  goal  but  until  recently  there  was 
no  experimental  verification.  Two  different  configurations  were  used;  one  for 
simplicity,  the  other  for  accuracy.  The  configurations  are  shown  in 
Figures  15a  and  b,  respectively. 
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Figure  5.  Discharge  Ar  Meter,  g/s  units 
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Figure  6.  Discharge  O2  Meter,  g/s  units 
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Figure  7.  Counterflow  Nt  Meter,  g/s  units 
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Figure  9.  Counterflow  Meter,  cm^/min  units 
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Figure  12.  Discharge  0?  Meter,  cmJ/min  units 
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Figure  13.  Discharge  Ar  Meter,  praol/s  uni 
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Figure  15.  Field-of-View  Appratus  Configurations 


The  apparatus  configuration  shown  in  Figure  15a  was  the  easiest  to  con¬ 
struct  and  so  was  tried  first.  With  the  detector  placed  along  the  horizontal 
diameter  of  the  reaction  cell  at  the  position  the  mirror  would  occupy,  the 
curve  in  Figure  16  was  obtained.  This  figure  shows  the  field  of  view  to  be 
~6  to  8  cm  wide.  Since  the  mirror  is  10  cm  wide,  it  is  clearly  underfilled. 
With  the  source  moved  2  in.  farther  away  from  the  lens,  Figure  17  was 
obtained.  The  intensity  has  increased  in  the  field  of  view  but  the  width  of 
the  field  has  decreased.  Figure  18  illustrates  the  results  for  a  position  of 
the  source  displaced  out  of  the  line  of  sight  by  1/2  in.  The  field  of  view  is 
offset  and  the  intensity  reduced. 

The  studies  with  this  conf iguration,  although  informative,  do  not  accu¬ 
rately  represent  the  true  field  view,  what  is  required  is  a  configuration 
similar  to  that  in  Figure  15b.  The  monochromator  was  put  in  zeroth  order  with 
the  slits  wide  open  (3  mm)  and  a  source  placed  at  the  position  of  the  mirror. 
Using  a  1/8  in.  square  source,  Figure  19  was  obtained.  This  figure  shows  a 
narrower  field  of  view  than  Figures  16  through  18.  Signal- to-noise  limita¬ 
tions  prevented  more  than  one  full  decade  of  data.  Nevertheless,  the  results 
indicate  that  the  mirror  is  underfilled  and  that  COCHISE  is  unlikely  to 
directly  observe  surface  emissions. 

2. 3  Detection  System  and  Signal  Processing 

2.3.1  Spectral  Detection  System 

The  detection  system  is  isolated  from  the  reaction  system  by  an  aluminum 
bulkhead  which  blocks  scattered  radiation  from  the  cell  and  prevents  the 
passage  of  reagent  gases  into  the  detection  system.  Infrared  radiation 
emanating  from  a  cylinder  ~8  cm  in  diameter  centered  along  the  axis  of  the 
reaction  cell  (or  from  the  blackbody  source  in  the  center  of  the  end  mirror, 
in  the  case  of  calibration  measurements)  passes  through  the  reaction  cell  lens 
into  the  detection  system  via  a  7.6  cm  AR-coated  germanium  window  in  the 
center  of  the  bulkhead.  The  radiation  then  passes  through  an  order-sorting 
filter  and  an  optional  tuning  fork  chopper  into  a  cryogenic  grating  mono¬ 
chromator,  where  it  is  dispersed  and  focussed  onto  a  liquid-helium  cooled 
detector. 
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2.3.2  Foreoptics 


A  series  of  long-wavele.  gth-pass  and/or  bandpass  filters  is  deployed  by  a 
specially  constructed,  solenoid-driven,  rotating  filter  wheel.  For  most 
applications  to  date,  long-wavelength-pass  filters  with  sharp  short-wavelength 
cutoffs  at  nominally  2,  4,  and  8  ym  have  been  used  to  isolate  the  correspond¬ 
ing  first-order  spectral  regions  of  2-4,  4-8,  and  8-16  ym. 

A  23  Hz  tuning  fork  chopper  (American  Time  Products  type  L40  or  Philamon 
type  TFLC  34CH-7178)  positioned  in  front  of  the  monochromator  entrance  slit 
can  be  used  to  modulate  the  incoming  radiation.  This  chopper  is  used  in  all 
blackbody  measurements  and  is  optional  in  the  chemiluminescence  measurements, 
as  will  be  described  below. 

2.3.3  Cryogenic  Monochromator 

In  the  development  of  a  cryogenic  spectrometer  system  for  the  COCHISE 
facility,  a  scanning  grating  monochromator  proved  to  be  the  most  straight¬ 
forward,  reliable,  and  economical  approach.  The  instrument  chosen  was  a  0.5m 
asymmetric  Czerny-Turner  system  (Minuteman  305  CM),  which  is  mounted  on  a 
Invar  frame  to  minimize  distortion  at  low  temperatures.  The  spherical 
collecting  and  focusing  mirrors  and  the  grating  are  mounted  in  Invar  supports, 
which  are  in  turn  mounted  directly  to  the  basic  frame.  The  entire  assembly  is 
covered  with  a  light-tight  aluminum  housing  and  is  held  at  ~20  K  during  normal 
operation. 

Two  interchangeable  replica  gratings  are  employed,  their  selection 
depending  on  the  wavelength  region  of  interest.  Both  gratings  are  6.4  cm 
square,  which  gives  the  instrument  an  aperture  ratio  of  f/6.9.  For  wave¬ 
lengths  below  8  ym,  a  grating  blazed  at  3  ym  with  150  lines/mm  (at  room  tem¬ 
perature)  allows  a  nominal  first-order  reciprocal  dispersion  of  0.013  ym/mm  in 
the  exit  slit  plane;  for  longer  wavelengths,  a  grating  blazed  at  10  ym  with 
75  lines/mm  gives  0.027  ym/mm.  Thus,  for  a  slit  width  of  0.5  mm,  the 


resolution  element  is  ~2.6  cm“^  at  5  ym  and  ~1.4  cm“^  at  10  ym.  Such  resolu¬ 
tion  permits  not  only  determination  of  detailed  vibrational  structure  for  most 
small  atmospheric  molecules  (typical  vibrational  spacings  of  ~25  cm-1  )  but 
also  the  observation  of  some  rotational  structure  in  many  cases. 

The  entrance  and  exit  slit  assemblies  consist  of  rotatable  disks  contain¬ 
ing  a  series  of  fixed  2  cm  high  slits  ranging  in  width  from  0.1  to  3.0  mm 
(2  mm  exit  slit  has  sapphire  insert);  the  slits  can  be  independently  indexed 
by  manual  controls  outside  the  vacuum  chamber.  The  grating  is  rotated  by  a 
conventional  sine-drive  mechanism  using  a  precision  ball  screw  designed 
specifically  for  cryogenic  service;  the  scan  drive  controls  are  located 
outside  the  vacuum  system.  Three  shafts  penetrate  the  outer  and  inner  walls 
of  the  system  and  are  coupled  to  the  slit  assemblies  and  grating  drive  screw 
via  stainless  steel  bellows. 

2.3.4  Cryogenic  Infrared  Detector 

The  IR  preamplifier-detector  module  was  designed  and  fabricated  at  the 
Electrodynamics  Laboratory  of  Utah  State  University.  The  photovoltaic  detec¬ 
tor  element  is  an  arsenic-doped  silicon  cube,  3.0  mm  on  an  edge  (Santa  Barbara 
Research  Center);  an  f/0.25  parabolic  reflector  images  the  monochromator  exit 
slit  onto  the  detector  and  thus  permits  collection  of  virtually  all  radiation 
passing  through  the  slit.  A  contiguous  preamplifier  unit  utilizes  an  inte¬ 
grated  JFET  operational  amplifier  (Burr-Brown  3521R)  in  a  direct-coupled 
negative-feedback  operational  scheme  that  previously  has  been  applied  success¬ 
fully  to  ultrahigh  impedance  cryogenic  detectors.1^1  The  module  also  contains 
an  IR  emitting  diode,  temperature  monitors,  and  associated  electronic  elements 
and  devices.  The  preamplifier  and  detector  are  maintained  at  temperatures  of 
~220  and  9  K,  respectively,  by  proportional  control  systems.  An  external 
console  provides  system  control  functions,  detector  bias,  signal  conditioning 
circuitry,  and  output  indicators  and  terminal.  This  detector,  when  operated 
at  ~9  K,  is  sensitive  to  1.2  to  22  ym  radiation  and  has  a  peak  noise  equiva¬ 
lent  power  of  ~10-^  W/Hz^^. 
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Liquid  helium  cooling  of  the  detector  package  is  provided  by  a  20-liter 
liquid  helium  vessel  suspended  from  a  port  on  the  top  of  the  vacuum  chamber. 
The  neck  of  the  vessel,  which  serves  as  a  f ill-and-vent  line,  is  flexibly  con¬ 
nected  by  a  bellows  to  a  port  in  the  20  K  shroud  to  effect  a  light-tight  pene¬ 
tration.  Thermal  coupling  between  the  detector  and  the  liquid  helium  vessel 
is  accomplished  by  a  copper  bar  and  thermal  strap,  with  sapphire  electrical 
standoffs  to  isolate  the  detector-preamplifier  case. 

2.3.5  Signal  Processing 

The  detector  output  is  processed  by  phase-sensitive  detection  to  discrim¬ 
inate  against  random-phase  background  noise.  The  radiation  reaching  the 
detector  is  modulated  in  one  of  two  ways.  In  most  chemiluminescence  observa¬ 
tions,  the  microwave  discharges  are  pulsed  at  ~23  Hz  with  a  square  wave, 

50  percent  duty  cycle;  since  chemical  reaction  does  not  occur  in  the  absence 
of  the  reactive  species  formed  in  the  discharge  plasma,  this  results  in  a 
square  wave  oscillation  in  the  chemiluminescence  signal  (provided  the  chopping 
frequency  is  slow  enough  relative  to  the  flow  speed  so  that  diffusional  blur¬ 
ring  of  the  waveform  does  not  occur).  An  alternate  method  of  data  collection 
is  to  operate  the  discharge  continuously  and  modulate  the  radiation  with  the 
chopper.  However,  the  pulsed  discharge  method  has  the  advantage  of  discrimin¬ 
ation  against  stray  radiation,  which  is  not  directly  linked  to  species  pro¬ 
duced  in  the  discharge  plasma.  The  chopper  is  also  used  to  modulate  the 
radiation  from  the  internal  blackbody  source  during  calibration  measurements. 

The  ac  detector  output  passes  through  an  optional  bandpass  amplifier 
(PARC  Model  113)  and  into  a  lock-in  amplifier  (PARC  Model  124)  which  is 
synchronized  either  to  the  frequency  of  the  chopper  pickup  coil  or  to  the 
external  pulse  generator  used  to  pulse  the  microwave  discharges.  The  final 
demodulated  signal  is  digitally  averaged,  displayed,  and  stored  in  real  time 
by  the  PDP-8/E  computer.  The  combination  of  grating  scan  rate,  electronic 
time  constant,  and  digital  averaging  for  a  given  spectral  scan  is  carefully 
chosen  so  that  there  are  at  least  seven  averaged  data  points  and  three  time 
constants  per  spectral  resolution  element;  these  criteria  insure  that  the 
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apparent  spectral  detail  is  indeed  limited  by  the  instrument  resolution.  The 
data  are  initially  stored  on  disk  by  the  PDP-8/E  computer  and  are  subsequently 
transferred  via  magnetic  tape  to  other  computers  for  detailed  analysis. 


2.3.6  Spectral  Responsivity 

The  spectral  responsivity  of  the  COCHISE  optical/detection  system,  as 
determined  from  calibrations  with  the  internal  blackbody  source,  is  plotted  as 
a  function  of  wavelength  in  Figure  20.  The  2  to  16  pm  spectral  region  was 
traversed  by  means  of  two  gratings  and  three  long-wavelength-pass  filters. 

The  solid  curves  were  obtained  using  a  10  pm  blazed  grating.  First-order 
measurements  employ  a  4  pm  long-pass  filter  (5  to  9  pm)  and  an  8  pm  long-pass 
filter  (8  to  16  pm),  while  second-order  measurements  can  be  made  using  a  4  pm 
long-pass  filter  together  with  a  sapphire  window  (4  to  8  pm).  The  dashed 
curves  were  obtained  in  first  order  using  a  3  pm  blazed  grating  with  2  and 
4  pm  long-pass  filters.  Since  the  noise  level  of  the  instrumentation  is 
~10”6V,  the  maximum  responsivity,  at  ~10  pm,  corresponds  to  an  optimum  NESR 
near  10“12W  cm-2  sr'1  pm-1  or  ~10^  photons  cm”2  s”1  pm-1  for  the  effective 
optical  path  length  of  100  cm.  For  typical  species  of  interest  (e.g.,  03 ) , 
the  transition  probability  is  ~10  s”^,  and  the  bandwidth  is  ~1  pm;  thus  the 
minimum  detectable  concentration  of  such  vibrationally  excited  species  is  of 
the  order  of  105  molecules  cm”2.  This  unsually  high  sensitivity  at  such  long 
wavelengths  is  due  to  the  absence  of  noise  from  thermal  background  radiation 
and  to  the  long  optical  path  which  is  available  for  the  measurements. 

2.4  Refrigeration  System 

COCHISE  requires  a  rugged,  high  capacity  refrigeration  system  in  order  to 
cool  down  and  maintain  a  chamber  temperature  of  ~16  K.  The  refrigeration  sys¬ 
tem  chosen  was  a  Cryogenic  Technology  Model  1400  which  is  capable  of  providing 
gaseous  He  at  the  required  temperature  and  mass  flow.  Pictorial  diagrams  of 
the  refrigerator  can  be  found  in  the  system  operation  manuals  which  are  kept 
in  the  COCHISE  work  space.  A  diagram  of  the  plumbing  of  the  refrigeration 
system  illustrating  those  components  relevant  to  COCHISE  operation  is  shown  in 
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Figure  20. 


Approximate  Spectral  Responsivity  of  COCHISE  Optical  System 
for  Two  Grating  Blaze  Wavelengths  and  Three  Order-Sorting 
Long-Pass  Filters  as  Indicated.  Curves  A,B,D,E  were 
obtained  in  first  order;  curve  C  was  obtained  in  second 
order  using  a  sapphire  long-wavelength  cutoff  filter. 


Figure  21.  Figure  22  shows,  in  greater  detail,  the  plumbing  of  the  refriger¬ 
ator  lines  entering  the  chamber.  Each  of  the  indicated  valves  (shroud,  heat 
sink  plate,  and  reaction  cell)  allow  individual  control  of  refrigerant  flow  to 
these  subsystems.  These  valves  are  physically  located  underneath  the  COCHISE 
chamber  and  are  evident  in  Figure  2.  Labelled  photographs  of  the  refrigerator 
system  may  be  found  in  the  manuals  stored  in  the  COCHISE  work  space. 

Typical  operation  of  the  refrigerator  system  is  decribed  below.  The 
system  operates  in  two  basic  modes:  normal  and  cooldown.  Normal  operation 
will  be  described  first: 

Normal:  (low  speed,  1  or  2  compressors) 

•  He  from  supply  manifold  and  low-P  return  gas  from  RF  (5  to  9  psig) 
go  into  compressor(s) . 

•  Compressed  He  is  regulated  at  230  to  240  psig  by  "bypass  regula¬ 
tor";  excess  goes  through  recovery/bypass  valve  (in  recovery  posi¬ 
tion)  into  storage  tank. 

•  Storage  tank  should  be  at  25  to  50  psig.  If  pressure  is  too  low  (as 
when  gas  gets  cold  in  cooling  lines,  then  some  bleeds  from  tank  into 
low  P  side),  the  low  side  P  cannot  be  controlled  by  low  side 
regulator.  If  this  happens,  let  some  in  from  the  supply  manifold. 
This  gas  goes  through  the  compressor ( s )  and  the  excess  is  shoved 
through  the  bypass  regulator  into  the  storage  tank. 

•  High  P  return  He  enters  four-stage  heat  exchangers.  Valve  V806 
controls  amount  which  goes  through  LN2  first  stage.  Other  stages 
use  heat  transfer  with  outgoing  He. 

•  Cooled  compressor  gas  undergoes  expansion  in  parallel  engines. 

Inlet  and  outlet  valves  for  each  engine  are  operated  by  big  cam  on 
top  of  RF  cabinet.  These  valves  need  to  be  adjusted  so  that  they 
open  and  shut  properly;  adjust  best  when  cold. 
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Figure  21.  Refrigerator  System 
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Figure  22.  COCHiSE  Refrigera 


P  drop  between  high  side  and  engine  inlet  should  be  ~(235  to  180) 
psig,  the  smaller  the  better.  Too  much  P  drop  is  robbing  you  of 
efficiency,  and  may  be  due  to  freeze-up  of  contaminants  in  heat 
exchangers. 

Low  P,  cold  (~10  K)  He  passes  through  load  and  returns  to  RF  at 
~  1 5  to  20  K. 

Return  coolant  He  (low  P)  passes  through  heat  exchangers  and  cools 
high  P  gas  coming  from  compressor ( s ) .  For  best  efficiency,  the 
temperature  at  the  exit  of  the  heat  exchangers  should  be 
-15°  to  -20®C.  This  temperature  is  controlled  by  valve  V806,  which 
diverts  some  of  the  high  P  gas  through  as  LN2  first  stage.  If  this 
valve  is  open  too  far,  there  is  too  much  cooling  by  LNg  and  not 
enough  warming  of  the  low  P  He.  The  result  is  that  LN2  is  wasted 
and  the  low  P  He  is  too  cold  going  to  the  compressors.  In  this 
event,  the  Cu  pipe  will  frost  up.  The  best  setting  for  V806  is 
~50  percent  open  (fully  open  =  9  turns). 

Warmed,  low  P  (5  to  9  psig)  He  returns  to  the  compressor  and  gets 
compressed;  the  cycle  starts  over  again. 

Reducing  low  side  P  from  9  to  4  psig  decreases  the  efficiency  but 
may  be  easier  on  the  compressors.  This  is  an  attempt  to  prolong  the 
life  of  the  seal  on  the  compressor  inlet,  the  most  failure-prone 
part  of  the  compressor.  Use  when  load  is  small,  i.e.,  no 
discharges. 

Note:  Never  operate  <4  psig;  operation  <0  psig  can  cause  damage  to 

the  compressors. 

If  high  P  side  exceeds  -245  psig,  the  compressor  relays  will  start 
to  make  a  godawful  noise  as  they  cut  in  and  out  in  order  to  not 
exceed  maximum  P.  Keep  the  setting  below  240,  since  this  tends  to 
drift  up  over  several  hours.  235  is  best. 
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Cooldown :  (high  speed,  1  compressor:  an  option)  but  typical  operation 


now  low  speed,  2  compressors 

Recovery  system  set  up  as  in  normal ;  high  P  He  from  compressors 
regulated  at  ~230  to  210  psig  and  storage  tank  stays  at 
25  to  50  psig. 

High  P  He  from  compressors  goes  entirely  through  LN2  first  stage. 
This  is  done  by  closing  the  red-handle  ball  valve  up  inside 
plumbing.  This  is  all  the  precooling  that  the  gas  gets,  since  there 
is  no  return  gas  on  the  other  side  of  the  heat  exchangers. 

High  P  He,  cooled  by  LN2,  undergoes  expansion  in  parallel  engines. 
Optimum  P  drop  is  the  same  as  in  normal . 

Cold,  low  P  He  goes  through  system  load. 

The  gas  returning  from  the  load  bypasses  the  heat  exchangers  and 
goes  directly  to  the  compressors,  as  follows.  Return  valve  on  RF  is 
closed.  All  He  from  load  passes  through  Cu  tubing  out  of  plenum, 
through  the  parallel  heaters  (variable  control),  through  the  bypass 
ball  valve  (which  is  closed  for  normal )  next  to  the  platform  and 
directly  to  the  compressors.  The  return  gas  and  heat  exchanger  exit 
gas  temperature  monitors  are  bypassed  and  therefore  their  readings 
are  meaningless.  The  bypass  is  continued  until  the  gas  returning 
from  the  load  is  cold  enough  to  do  some  good  in  the  heat  exchangers. 

In  bypass  mode,  plug  in  bypass  line  heater.  Th.s  heater  is  thermo- 
statted  for  50  °F  to  keep  this  line  from  freezing. 


Bypass  mode  should  be  raa in ta ined  until  the  shroud  is  -100  K. 


When  shroud  warm  end  ~100  K  or  at  the  end  of  the  day,  begin  tran¬ 
sition  to  normal  (turn  off  heater,  close  bypass).  REMEMBER  TO 
TURN  OFF  HEATER  SEVERAL  MINUTES  BEFORE  CLOSING  THE  BYPASS  VALVE, 
OR  AN  EXPLOSION  WILL  RESULT 


•  Once  bypass  is  closed,  system  will  cool  fast  except  for  spectrom¬ 
eter.  To  speed  this  up,  backfill  with  He.  Watch  P  in  vacuum  space 
as  well  as  in  spectrometer  and  do  not  exceed  <  10“^  Torr  on  ioniza¬ 
tion  gauge.  This  will  correspond  to  ~10Q  mTorr  in  the  spectrometer 
which  will  slowly  decrease.  Repeat  every  time  this  gets  below 

~ 1 0  mTorr. 

•  IMPORTANT:  remember  to  purge  LHe  dewar  with  He  during  the  entire 
cooldown  in  order  to  prevent  icing. 

•  Start  with  LN2  on  automatic  throttle.  After  1  hr,  switch  to 

un throttled  LN2.  Reinstall  LN2  flow  controller  when  switching  over 
to  normal. 

2 . 5  Temperature  Diagnostics  and  Control 

Not  all  of  the  COCHISE  components  can  be  operated  at  cryogenic  temper¬ 
ature.  The  gas  lines  (discharge  heat  sinks,  counterflow  and  discharge  mani¬ 
folds,  and  pressure  port  plate)  are  typically  operated  at  90  K  to  avoid  con¬ 
densation  of  reagent  gases.  Figure  23  shows  the  vapor  pressure  of  all  the 
gases  used  in  COCHISE  experiments  at  various  temperatures.  Operation  at  90  K 
is  adequate  for  most  commonly  used  gases  (N2,  CO,  Ar,  and  02)  but  much  higher 
temperatures  are  required  for  NO,  N20,  and  C02. 

The  optical  components  (bulkhead  window,  reaction  cell  lenses,  reaction 
cell  mirror)  are  typically  operated  40  to  80  K  in  order  to  prevent  condensa¬ 
tion  of  reagent  and  product  gases.  Since  these  components  are  within  the 
field  of  view  operation  at  40  to  60  K  is  preferred  to  minimize  the  thermal 
background. 
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Figure  23.  Vapor  Pressure  of  Reagent  Gases; 
melting  temperature 


indicates 


The  various  channel  assignments  for  the  temperature  monitor/controls  are 
listed  in  Table  1 .  Control  and  readout  is  possible  only  when  the  program 
MOTHER  is  running.  The  teletype  prints  the  refrigerator,  shroud,  reaction 
chamber,  heat  sink  and  spectrometer  temperatures  every  5  min.  Output  of  these 
sensor  readings  are  necessary  during  cooldown  and  warmup  and  cannot  be 
altered.  Readings  of  any  of  the  other  sensors,  including  those  output  to  the 
teletype,  may  be  displayed  on  the  two  console  screens  by  setting  the  A/D 
numbers  into  the  six  locations  beneath  the  left  console.  A  total  of  six  chan¬ 
nels  may  be  monitored  simultaneously.  One  note  of  reference,  the  second  chan¬ 
nel  can  accept  only  an  even  A/D  valve. 

The  temperature  sensors  used  in  COCHISE  have  a  variable  resistance  with 
temperature.  Figure  24  shows  a  typical  plot.  The  resistance  is  linear  over 
the  temperature  range  40  to  400  K.  Below  40  K  the  resistance  is  decidedly 
nonlinear  and  accurate  temperature  readings  are  afforded  only  above  25  K. 

This  is  adequate  for  most  of  the  COCHISE  sensors.  However,  a  much  more  accu¬ 
rate  temperature  reading  in  the  5  to  20  K  region  is  required  for  the  detector. 
This  is  achieved  with  a  sensor  which  increases  in  resistance  with  decreasing 
temperature.  A  plot  of  the  temperature  versus  sensor  resistance  is  shown  in 
Figure  25.  This  sensor  is  mounted  on  the  copper  bar  which  is  thermally 
attached  to  the  detector.  Readout  is  achieved  through  a  port  marked  LHe  on 
the  console  connected  to  a  DVM. 

The  heater  controls  are  located  on  the  COCHISE  console  panel.  A  power 
switch  is  located  to  the  left.  Each  heater  unit  is  composed  of  a  Variac  and 
Ihree-posi tion  switch.  In  the  up  position,  the  heater  is  controlled  by  the 
computer.  With  the  switch  in  the  neutral  position,  the  heater  is  off  and  in 
the  down  position,  manual  control  is  permitted.  The  Variacs  may  be  set  to  any 
desired  position  allowing  the  heating  rate  to  be  adjusted.  The  spectrometer 
heater  is  not  located  on  the  console,  it  is  located  on  the  shelf  above  the  gas 
cylinders.  This  heater  is  also  controlled  by  a  Variac  but  has  only  a  manual 
control.  Computer  control  is  not  necessary  since  this  heater  is  only  used 
during  warmup  to  keep  the  spectrometer  warmer  than  the  shroud  and  reaction 
chamber  thus  preventing  condensation  of  gases. 


Table  1.  Channel  Assignments  for  CCCKISE  Temperature  M^n  i  tor  ./Control 


A/D 

Heater 

Location 

Sensor 

Cabling 

Identity 

Sensor 

I/O 

Port 

37 

Shroud,  warm  end 

HP  2 

BV  (3,4) 

36  ** 

Reaction  cell 

CT  60 

BI  (5,6) 

35 

8 

Discharge  gas  manifold 

HP  8 

BIV  (1,2) 

34 

Discharge  master  heat  sink  plate 

CT  53 

BV  (7,8) 

33 

Refrigerant  return 

CT  61 

BVI  (1,2) 

32 

9 

Bulkhead  window 

HP  3 

BIV  (7,8) 

31 

10 

Reaction  cell  lens 

CT  69 

BVI  (5,6) 

30 

7 

Reaction  cell  mirror 

CT  56 

BVI  (7,8) 

27* 

4 

Discharge  heat  sink  #8 

CT  65 

Bill ( 7 , 8 ) 

26 

(4) 

Discharge  heat  sink  #7 

CT  64 

Bill (5,6) 

25* 

3 

Discharge  heat  sink  #6 

CT  58 

Bill ( 3 , 4 ) 

24 

(3) 

Discharge  heat  sink  #5 

CT  57 

BIII( 1 ,2) 

23* 

2 

Discharge  heat  sink  #4 

CT  55 

BII  (7,8) 

22 

(2) 

Discharge  heat  sink  #3 

CT  54 

BI I  (5,6) 

21* 

1 

Discharge  heat  sink  #2 

CT  52 

BII  (3,4) 

20 

(1) 

Discharge  heat  sink  #1 

CT  51 

BII  (1,2) 

17 

1 1 

Inner  hreat  exchange  plate 

HP  5 

BIV  (5,6) 

16 

Shourd,  cold  end  (spectrometer 

valve  plate) 

HP  7 

AI  (5,6) 

15 

Outer  heat  exchange  plate 

HP  1 

BV  (5,6) 

14 

6 

Pressure  port  plate 

HP  4 

BIV  (3,4) 

13 

5 

Counterflow  manifold 

CT  62 

BV  (1,2) 

12  *** 

Spectrometer  interior 

HP  6 

AI  (3,4) 

1  1 

Spare 

CT  67 

BI  (7,8) 

10 

12 

Blackbody  source 

CT  59 

BVI  (3,4) 

♦Heater  channels  1,2, 3, 4  are  controlled  by  temperature  data  from  heat 
sinks  2, 4, 6, 8:  A/D  channels  21,23,25,27. 

♦♦Reaction  cell  temperature  is  controlled  by  separate  analog  device. 


♦♦♦Spectrometer  heating  is  accomplished  via  CW  Variac  on  side  shelf 
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2.6  The  Computer  System 


The  PDP-8E  computer  system  is  the  control  and  data  acquisition  system  of 
COCHISE.  The  program  MOTHER  affords  all  of  the  control  and  data  acquisition 
tasks  necessary  to  operate  COCHISE.  This  section  discusses  three  aspects 
relevant  to  operation  of  the  computer  system:  1)  booting  the  computer, 

2)  Tektronix  CRT  commands  for  running  MOTHER  and  data  file  manipulation,  and 

3)  teletype  commands  to  open  data  files,  start  data  collection,  and  affect 
temperature  control. 

The  booting  procedure  is  shown  in  Table  2  and  is  relatively  simple  and 
straightforward.  If  this  booting  system  fails,  try  again.  If  the  computer 
will  still  not  operate,  there  may  be  a  problem  in  the  computer  hardware  or  in 
the  disk  drive.  A  service  call  to  Digital  Corp.  (273-6630)  may  be  necessary. 
The  identification  number  of  the  computer  system  is  8E  (system  type),  6631 
(serial  number). 

The  Tektronix  CRT  commands  and  teletype  commands  are  shown  in  Tables  3 
and  4,  respectively.  It  is  important  to  note  that  Tektronix  CRT  commands  can 
only  be  used  when  MOTHER  is  not  running,  only  teletype  commands  can  be  entered 
in  MOTHER. 

The  computer  system  has  a  total  of  four  disks,  the  two  sides  of  the  one 
on  top  and  the  two  sides  of  the  bottom  one.  They  are  numbered  sequentially 
RKAO,  RKA1 ,  RKA2,  and  RKA3 .  Disk  0  is  used  to  store  all  of  the  COCHISE 
programs  including  MOTHER.  Data  storage  should  not  be  attempted  on  this  disk. 
Data  storage  is  automatically  perfomed  on  Disk  1.  Disks  2  and  3  are  used  for 
backup. 


Table  2.  Computer  Startup  Procedure 


1 .  Turn  power  on 

2.  Disk  drive  0  to  run,  wait  until  "RDY-ONCYL"  lights  illuminated  - 
"WT  PROT”  light  off 

3.  Switch  register  to  all  zeros,  press  "Ext  Addr  LD"  switch, 
verify  "EMA"  lights  off 

4.  Switch  register  30  octal,  "7,8"  up  -  press  "Addr  Load"  switch 

5.  Switch  register  6743  octal  -  "0,1,3,4,5,6,10,11"  up  -  raise 
"Dep"  switch 

6.  Switch  register  5031  octal  -  "0,2,7,8,11"  up  -  raise  "Dep"  switch 

7.  Repeat  Step  4 

8.  Press  "Clear"  switch 

9.  Press  "Cont"  switch 

10.  Monitor  dot  should  appear  on  Tektronix  4010 

11.  Type  in  date  in  form  XX-YYY-ZZ  (XX  =  day,  YYY  *  month,  ZZ  =  year 
and  acceptable  years  are  70-79) 
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Table  3.  Tektronix  CRT  Commands 


Running  MOTHER: 

DA  XX-YYY-ZZ  (CP)  -  Set  date,  day  (XX),  month  (YY),  and  year  (ZZ) .  The 

acceptable  entries  for  year  are  70-79. 

DA  (CR)  -  Ask  for  the  current  date 

R  MOTHER  (CR)  -  Run  MOTHER 

Data  File  Mani.  ulation: 

Dir  Tel:  <RKA1 :/F=5  (CR)  List  contents  of  RKA1  disk  to  printer 

Copy  RKA3 : <RKA1 : XXXXXX. YY ,  etc.  (CR).  Copy  file  (up  to  5)  from  Disk  l' 

to  Disk  3 

R  CVRTBA  (CR)  -  Converts  the  binary  DB  files  to  ASCII  DA  files.  Enter 
file  name  including  DB  with  no  dot  or  carriage  return. 
The  files  must  be  on  Disk  1  (RKA1 ) .  This  program  will 
have  to  be  run  for  each  DB  file  to  be  converted. 

Copy  TEL: <RKA1 : * . DB  (CR)  -  Copy  comments  of  all  DB  files  to  terminal. 

Use  only  on  DB  files. 

Delete  RKA 1 : * . DA  (CR)  -  Delete  files  from  a  disk.  Individual  file 

names  may  be  entered  or  *  for  all. 

Verification  will  be  asked  for  each  file. 

May  be  used  on  both  DA  and  DB  file  types. 

R  TVDATA  (CR)  -  To  inspect  an  individual  file  (screen  will  prompt 
for  name) 

Tape  Operations: 

To  load  tape: 

Plug  in  power  cord  to  tape  drive 
Mount  tape 
Hit  load 

Copy  MTAO : </Z  (CR)  Zeros  directory  of  tape 

Copy  MTAO : <RKA1 : * . DA  (CR)  Copies  files  to  tape.  Only  DA  files  are 

desired  on  tape. 

Copy  MTAO:*. DA  (CR)  Lists  contents  of  tape  (*  may  be  used  for  DA) 

Type  MTAO : XXXXXX. DA  (CR)  List  file  from  tape  on  terminal,  including 

comments . 
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Table  4.  Instructions  for  Data  Acquisition  in  COCHISE  MOTHER 


Data  acquisition  is  performed  utilizing  Channel  0,  GAIN  =  0  (i.e., 

10  V.F.S);  Channel  0  is  sampled  at  a  120  Hz  rate  and  averaged  over  NP 
samples.  NP  is  set  at  120  initially,  so  this  leads  to  a  data  rate  of 

1  Hz,  where  each  of  these  stored  data  points  is  an  average  of  NP(120) 
samples.  For  example,  if  the  user  sets  NP  =  60,  then  the  data  rate  is 

2  Hz  where  each  of  those  stored  data  points  is  an  average  over  60  samples. 
Thus  the  "stored  data"  rate  is  (120/NP)  Hz,  and  each  of  the  stored  data 
points  is  the  average  of  NP  samples. 

Teletype  commands. 

DF  (CR)  Which  names  and  opens  a  file  for  data  storage.  MOTHER'S 

response  to  the  ;‘DF"  command  is: 

$$$$$$ 

JJJY#  # 

where  JJJ  is  the  Julian  day  of  the  year,  Y  indicates 
the  year  197Y,  and  ##  is  the  serial  number  of  the  data 
file  for  that  date.  Calendar  information  is  obtained 
from  the  system  date.  The  actual  name  of  the  data  file 
will  be  JJJY##. DA. 

SCOMMENTS  Input  comments  which  will  be  written  into  the  data  file. 

Each  line  of  comments  must  be  preceded  by  the  $  sign. 

NP  XXXX  (CR)  Inputs  new  value  of  NP;  23<NP<4095 

DS  (CR)  Instructs  MOTHER  to  being  data  sampling 

DE  (CR)  Instructs  MOTHER  to  stop  data  sampling,  write  the 

stored  data  into  the  file  opened  by  the  last  DF  command 

and  close  the  file. 

TM  nn  XXX. X  (CR)  Sets  A/D  station  nn  to  desired  temperature 

TI  mmmm  (CR)  Sets  clock  (24  hr)  to  correct  local  time 

.(text)  Enters  comments  on  teletype  output  only  -  no  storage 

by  computer 

Control/C  Stops  running  MOTHER,  goes  back  to  monitor 


3.  TYPICAL  OPERATION  AND  MAINTENANCE 


The  purpose  of  this  section  is  to  provide  detailed  instructions  and 
checklists  such  that  someone  familiarized  with  COCHISE  will  be  able  to  operate 
the  system  and  acquire  data.  This  section  is  not  intended  as  a  cookbook  for  a 
novice;  COCHISE  is  far  too  complicated  to  be  operated  independently  without 
prior  experience  and  the  potential  danger  of  damage  to  the  apparatus  resulting 
from  inexperience  or  an  error  in  judgement  is  substantial.  Nevertheless  this 
section  will  be  sufficiently  detailed  such  that  it  will  serve  as  a  valuable 
learning  tool  for  the  novice  and  a  useful  guide  for  the  experienced  operator. 

This  chapter  is  broken  up  into  several  subsections  including:  1)  sealing 
the  tank,  2)  pumpdown  and  cooldown,  3)  blackbody  calibrations,  4)  taking  dis¬ 
charge  data,  5)  overnight  and  warm-up  procedures,  and  6)  routine  maintenance. 

3 . 1  Sealing  the  Tank 

Closing  and  sealing  the  tank  is  itself  a  rather  simple  task  but  there  are 
several  procedures  to  be  performed  prior  to  sealing  the  tank  which  are  neces¬ 
sary  to  ensure  smooth  operation.  These  include  replacing  the  monochromator, 
verification  of  detector  isolation,  and  checking  all  sensors. 

The  procedures  for  spectrometer  installation  are  shown  in  Checklist  1. 
There  are  several  critical  elements  which  need  to  be  emphasized.  Ensure  the 
spectrometer  is  in  zeroth  order  before  putting  it  back  into  the  chamber.  This 
produces  the  greatest  accuracy  in  the  wavelengths.  Locating  zeroth  order  once 
the  chamber  is  cryogenic  results  in  inferior  wavelength  accuracy,  possibly  due 
to  non-uniform  changes  in  the  sine  bar.  Make  sure  all  cables,  heat  sink 
straps,  and  drives  have  been  reattached.  Use  cable  ties  at  the  connection 
points  to  make  sure  cables  hold  together.  Next,  make  sure  the  detector  is 
electronically  isolated  from  the  rest  of  the  chamber.  This  may  best  be 
checked  by  checking  the  resistance  between  the  detector  power  supply  housing 
(located  on  top  of  the  chamber)  and  the  chamber  itself  or  the  copper  bar. 
Infinite  resistance  should  result  in  both  cases.  If  not,  it  is  possible  that 
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Checklist  1 


Spectrometer  Installation  Checklist 


1.  Heater  cable  and  temperature  sensor  cable  from  interior  attached  to  their 
terminal  block. 

2.  Line  up  both  entrance  and  exit  slits  at  their  3  mm  positions  (widest). 

3.  Shine  flashlight  into  one  slit  and  verify  that  with  fiducial  marks  lined 
up  on  the  grating  drive  shaft  zero  order  image  emerges  center  0  in  the 
other  slit. 

4.  Insert  spectrometer  into  shroud,  centered  on  tracks  and  butted  against 
backstops  on  tracks. 

5.  Verify  that  external  controls  arc  set  to: 

a.  Entrance  slit  3.0 

b.  Exit  slit  3.0 

c.  Grating  drive  00000 

6.  Connect  chopper  drive  thermal  block. 

7.  Attach  connecting  shafts  to  slit  and  grating  drive  shafts,  verifying 
that  no  relative  motion  has  occurred  before  tightening  split  collars. 

8.  Attach  connecting  leads  to  internal  and  external  temperature  sensor 
blocks. 

9.  Connect  two-pin  spectrometer  base  heater  connector  to  mate  on  case 
wiring,  connect  four-pin  connector  to  mate  from  bulkhead. 

10.  Verify  that  nylon  in  sapphire  spacers  remain  in  detector  mounting  pads 
and  attach  detector  to  case  with  nylon  studs.  Tighten  no  more  than  1/8 
turn. 

11.  Attach  detector  cover  can  and  apply  mylar  foil. 

12.  Insert  molecular  sieve  tray  and  clamp  to  frame. 

13.  Install  three  hold-down  bolts;  tighten  1/2  turn. 

14.  Before  closing  door,  check  that  grating  drive  functions,  chopper 
operates,  slits  may  be  moved,  detector  passes  warm  test,  temperature 
sensors  give  correct  room  temperature  response. 

15.  Make  sure  microdot  detector  leads  are  screwed  on  tigh' or  microphonics 
will  result. 

NOTE:  When  removing  unit  from  shroud,  return  external  grating  drive  counter 

to  00000  before  disconnecting  drive  shaft. 
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some  of  the  sapphire  spacers  were  not  replaced  and  the  detector  is  effectively 
grounded.  This  situation  must  be  corrected  prior  to  sealing  the  tank. 

The  last  item  of  importance  is  the  verifiction  of  senscr  operation.  Turn 
on  the  computer,  run  MOTHER  and  display  all  of  the  sensors.  All  should  read 
room  temperature.  The  LHe  sensor  can  be  checked  at  the  marked  port  on  the 
COCHISE  panel  or  directly  in  the  chamber  (it  is  mounted  on  the  Cu  bar).  This 
sensor  should  read  only  a  few  ohms  at  room  temperature.  If  a  sensor  is  not 
operational,  check  the  resistance  through  the  ports  on  the  top  of  the  chamber. 
Frequently,  loose  port  connections  are  responsible  for  anomalous  sensor  read¬ 
ings.  If  the  port  test  fails,  isolate  and  test  the  sensor  directly  and 
replace  if  necessary. 

If  all  is  operational,  the  chamber  may  be  closed  up.  Check  indium  seals 
and  replace  entirely  or  only  broken  sections  (the  indium  seals  fit  in  grooves 
on  the  inner  chamber).  RTV  may  be  used  to  hold  the  indium  firm  while  bolting 
the  caps  in  place.  Bolt  on  the  insulated  covers  and  close  the  outer  door. 
Close  the  door  clamps  to  hold  the  outer  door  firmly  shut.  The  chamber  is  now 
ready  for  pumpdown. 

3 . 2  Pumpdown  and  Cooldown 

This  subsection  provides  detailed  checklists  for  the  pumpdown  and  cool¬ 
down  cycles  of  COCHISE.  Checklist  2  covers  pumpdown,  Checklist  3  covers  the 
cleaning  of  the  refrigerant  lines  (a  necessary  precursor  to  cooldown).  Check¬ 
list  4  lists  the  procedures  for  cooldown.  Checklist  5  illustrates  the  tech¬ 
nique  for  cooling  the  spectrometer  by  backfilling  with  helium,  and  Checklist  b 
shows  how  to  fill  the  detector  with  liquid  helium.  The  cleaning  of  the 
refrigerant  lines  may  be  performed  at  any  time  prior  to  cooldown  and  must  be 
performed  prior  to  every  run  to  remove  condensible  contaminants  (e.g.,  water, 
oil,  etc. ) . 
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Checklist  2 


Pumpdown  Procedures 


1.  Maxe  sure  right  angle,  bypass,  and  foreline  valves  are  closed  on  COCHISE 
central  panel  (see  Figure  26) . 

2.  Turn  on  roughing  pump  until  P  =  50p  on  foreline  gauge. 

3.  Open  bypass  to  rough-pump  the  chamber,  make  sure  the  spectrometer, 

bulkhead,  and  reaction  chamber  valves  are  open. 

4.  Ordinarily  the  diffusion  pump  has  oeen  left  under  vacuum.  If  not,  this 
must  be  roughed  out  as  well: 

a.  Open  foreline  valve  (output  of  diffusion  pump) 

b.  Continue  pumping  diffusion  pump  and  chamber  until  P  ~ 1  x  10~3  Torr 

5.  Turn  on  diffusion  pump  heaters  and  cooling  water  pump.  The  pump  will  take 

-1  hr  to  warm  up  before  it  begins  to  pump.  Sensors  will  shut  down  the 
diffusion  pump  if  it  gets  too  hot.  Open  right  angle  valve;  bypass  closes 
automatically. 


6.  Continue  pumping  for  2  to  3  days. 

7.  When  P  -1  x  10-^,  the  system  is  ready  to  ;e  cooled. 
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Checklisi  3 


Cleaning  the  Refrigerant  Lines 


1  . 

2. 


3. 


4. 

5. 

6 . 

7. 

8. 

9. 

10. 
1  1  . 
12. 

13. 

14. 

15. 
1b. 

17. 

18. 

19. 

20. 


Make  sure  chilled  water  is  running  to  the  compressors  (water  lines  loca¬ 
ted  above  each  compressor) . 

With  the  20  liter  Dewar: 

a.  Using  a  small  pump,  pumpdown  the  20£  Dewar  for  at  least  1/2  hour  to 
remove  H2O 

b.  Verify  with  a  TC  gauge 

c.  Backfill  with  helium  up  to  5  to  10  psi  to  prevent  an  air  leak  into 
the  compressor  lines  when  attaching  hoses. 

On  the  refrigerator: 

a.  Put  lock-in  bar  through  flywheel  to  prevent  turning 

b.  Put  in  jacking  rods  (two  on  each  side,  four  total)  to  make  sure 
valves  stay  open 

Pump  out  heat  exchanger  with  internal  pump  to  10-20y. 

Close  chamber  and  supply  valves  so  COCHISE  is  isolated. 

Open  chamber  bypass  valve. 

Turn  on  one  of  the  compressors  and  allow  circulation  for  several  minutes. 
Fill  20-liter  Dewar  with  LN2 

a.  Typically  20  liters  will  last  only  1  hr 

b.  Keep  a  filled  160-liter  LN2  Dewar  nearby 

Attach  Dewar  to  refrigerator  inlet  and  outlet  (near  floor)  while  the 
system  is  running. 

Close  valve  between  refrigerator  inlet  and  outlet  (near  floor) 

Run  each  compressor  for  -> 2  hr. 

When  one  is  done,  turn  off  and  turn  on  another. 

With  the  last  compressor  done,  it  is  time  to  clean  the  COCHISE  chamber 
1 ines . 

Open  return  gas  valves  (on  side  near  floor) 

Open  chamber  supply,  close  chamber  bypass.  Run  for  -2  hr. 

After  2  hr,  isolate  the  chamber  again  by  opening  the  chamber  bypass  and 
closing  the  supply  and  return  gas  bypass. 

We  are  now  back  to  the  original  scheme.  Open  cryosorb  valve. 

Disconnect  20-liter  Dewar  lines  and  replace  caps.  Must  be  done  guickly 
to  minimize  Me  loss.  Use  spray  lubricant  it  necessary. 

bhut  off  compressors,  remove  all  jacks  and  the  locking  bai. 

With  experience,  this  procedure  taxes  '  tuli  ‘lay. 


Checklist  4 


Cooldown  Procedures 


t.  Turn  on  computer  and  run  program  MOTHER 

2.  On  large  He  storage  tank: 

a.  Ensure  pressure:  25  to  35  psi 

b.  Otherwise  fill  with  He  from  a  cylinder  through  manifold  to  the  left 
of  the  frige. 

2a.  Start  He  purge  of  detector  LHe  reservoir. 

3.  Cool  down  refrigerator  first  (~1  hr): 

a.  Open  chamber  bypass,  turn  on  LN2  (line  pressure  should  be  ~40  psi) 

b.  Remove  LN2  sensor  apparatus  in  refrigerator  and  replace  with  bypass 
tube  (ensures  maximum  LN2  flow).  Leave  bypass  tube  in  until  the 
entire  system,  including  the  chamber,  is  out  of  bypass  mode. 

c.  Turn  on  two  compressors  on  low  speed  (#1,  #2,  or  #3) 

d.  Open  V806  valve  nine  turns  (full  open)  and  close  red-handled  valve 
(see  Figure  21).  This  ensures  maximum  LN2  cooling. 

4.  When  refrigerator  return  temperature  ~25  K,  put  system  on  chamber  bypass 
(begin  cooling  chamber,  ~1  day).  See  also  section  2.4  for  refrigerator 
operation. 

a.  Open  return  gas  bypass 

b.  Open  chamber  supply 

c.  Close  chamber  bypass 

d.  Turn  on  bypass  heaters  (preset  to  50°F) 

e.  Open  red-handled  valve  in  refrigerator  (see  Figure  21) 

f.  V806  remains  a  few  turns  open 

5.  When  the  shroud  temperature  gets  to  150  K  (exit  bypass  mode) 

a.  Shut  heater  off  and  wait  1  min  to  verify  pipes  are  cool  to  the  touch 

b.  Shut  off  return  gas  bypass 

6.  LN2  sensor  apparatus  inside  the  refrigerator  may  be  replaced 

7.  May  begin  to  cool  spectrometer  by  backfilling  with  He  when  shroud  is 
<50  K  (see  Checklist  5) 

8.  When  all  temperatures  are  <30  K,  the  detector  may  be  filled  with  LHe. 

This  procedure  should  take  no  more  than  4  days  and  the  final  chamber 
pressure  should  be  1  x  10“^  Torr.  Two  typical  cooldown  cycles  showing 
the  refrigerator,  shroud,  and  spectrometer  temperatures  are  shown  in 
Figure  27. 


9. 
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re  27.  Cooldown  Cycle 


Checklist  5 


Flushing  Spectrometer  with  Helium 

1.  Close  monochromator  and  bulkhead  valves  (on  central  panel) 

2.  Open  He  tank  valve  (on  central  panel) 

3.  Flush  He  line 

4.  Set  baratron  to  read  spectrometer  pressure 

5.  Turn  on  three-way  valve 

6.  Close  monochromator  and  bulkhead  valves 

7.  Adjust  fill  valve  until  pressure  =  50  to  60  microns 

8.  Close  He  tank  valve  after  each  fill 

9.  Repeat  until  desired  temperature  is  achieved  ( < 30  K) 

10.  Re-open  monochromator  and  bulkhead  valves  when  done 


Checklist  6 


Cooling  the  Detector  with  Liquid  Helium 


1.  Place  the  withdrawal  tube  in  the  LHe  supply  Dewar. 

2.  Open  detector  reservoir  off  gas  valve  to  vent  pressure  and  withdraw 
gauge  assembly  from  detector  LHe  reservoir. 

3.  Attach  the  metal  flexible  hose  to  withdrawal  tube  and  cool. 

4.  Insert  flexible  hose  end  into  detector  reservoir  and  tighten  O-ring  seal. 
Keep  off  gas  valve  open. 

5.  Attach  a  He  line  (<5  psi)  to  pressurize  Dewar  and  fill  reservoir  (main¬ 
tain  at  1  to  1-1/2  psig). 

5a.  Pre-cool  the  transfer  line  by  flowing  LHe  before  inserting  into  the 
reservoir  neck. 

6.  Insert  transfer  line  into  reservoir,  fill  and  withdraw  when  done. 

7.  Insert  LHe  reservoir  depth  detector  into  reservoir  and  determine  volume 
from  gauge  output  and  calibration  shown  in  Figure  28. 

8.  Remove  LHe  depth  detector  apparatus. 

9.  Reinstall  gauge  assembly  and  close  off  gas  valve. 

10.  Ensure  positive  He  pressure  on  reservoir  from  cylinder  to  prevent 
condensation. 

11.  This  procedure  will  take  an  hour  for  the  first  fill  and  30-45  minutes  for 
successive  transfers. 


56  7 


INCHES  READ  ON  LHE  METER 


3. 3  Blackbody  Calibrations 


Before  gas  is  introduced  into  the  chamber,  the  spectral  responsivity  of 
the  spectrometer  must  be  determined.  This  must  be  done  after  every  cooldown 
cycle  prior  to  operating  the  discharges  regardless  of  whether  or  not  the  spec¬ 
trometer  had  been  moved.  The  movement  of  COCHISE  components  during  cooldown 
(absolute  magnitude  is  unknown)  requires  recalibration  of  the  spectral 
responsivity  for  every  run.  The  calibration  is  achieved  by  taking  spectra  of 
a  blackbody  at  a  known  temperature  for  each  filter  to  be  used  (i.e.,  SWIR, 
MWIR,  etc.).  The  blackbody  used  for  this  purpose  is  one  of  the  330(1  resistive 
heaters  mounted  through  the  center  of  the  mirror  with'  its  flat  front  surface 
painted  black  for  maximum  emissivity. 

Taking  blackbody  spectra  is  relatively  routine;  the  chore  is  to  take 
enough  such  that  an  accurate  spectral  calibration  results.  In  this  regard, 
the  SWIR  2  to  2.5  gn  region  poses  the  only  serious  difficulty.  Even  at  the 
highest  operable  blackbody  temperature  (370  K,  operation  above  370  K  has  been 
observed  to  char  the  epoxy),  the  intensity  drops  off  rapidly  in  this  wave¬ 
length  region.  The  current  technique  to  obtain  a  SWIR  blackbody  spectrum  is 
to  increase  the  temperature  to  ~350  K  and  obtain  several  spectra  at  one  of  the 
lower  sensitive  scales  (200  or  500  mV).  A  few  spectra  are  then  taken  at 
higher  sensitivity  (50  mV)  and  these  sections  are  merged  with  the  reduced 
sensitivity  spectra  after  implementation  of  the  calibration  program.  The 
generalized  procedure  to  take  calibration  spectra  is  shown  in  Checklist  7. 
Sample  MWIR  and  SWIR  blackbody  spectra  are  shown  in  Figures  29  through  31. 
Figure  31  shows  a  spectrum  of  the  2  to  3  gm  section  of  the  SWIR.  Figures  32 
and  33  show  the  resulting  responsivity  curves.  Over  the  years,  these  curves 
have  not  deviated  by  greater  than  a  factor  of  four. 


Checklist  7 


Taking  Blackbody  Spectra 


1.  Heat  the  optics  to  ~40  K  and  the  blackbody  to  300  to  370  K. 

The  blackbody  may  require  3/4  hr  to  stabilize 

The  computer  command  to  set  the  temperature  is  shown  in  Table  4 

2.  Select  the  filter  (see  Checklist  8) 

3.  Turn  on  detector 

a.  Press  power  on  * 

b.  Wait  till  temperature  is  ~220  K 

c.  Turn  on  bias 

d.  Verify  bias  >  25V  otherwise  replace  batteries 

4.  Operate  in  dc  mode: 

a.  Turn  chopper  on 

b.  Display  output  on  scope 

c.  If  the  chopper  is  "chattering'1  turn  amplitude  down  (located  on  top 
of  the  chamber) 

d.  Output  of  detector  feeds  into  PAR  channel  A 

e.  PAR  reference  to  Kron-Hite  output 

5.  Select  slit  width  and  scan  rate,  scale  (mV)  and  RC  filtering  constant  and 
data  collection  at  (DP) 

a.  Table  5  shows  recommended  scan  rates,  DP  value,  and  filtering 
constants  for  various  slit  widths 

b.  Although  we  operate  in  first  order,  the  second  order  row  is  appro¬ 
priate  to  interpreting  Table  5, 

* 

6.  Prepare  to  scan: 

a.  Set  up  monochromator  near  X  start  -  typical  MWIR  value  is  7.2 

(3.6  ym) ,  SWIR  is  3.6  (1.8  ym)  .  \ 

b.  Set  up  data  file;  include  comments,  scan  rate,  DP  value,  X  start, 
mV  full  scale,  time  constant,  and  slit  width. 

c.  Start  the  monochromator  and  begin  data  collection  when  the 
monochromator  reaches  X  start. 

7.  Stop  data  collection  when  desired  by  entering  the  DE  command  followed  by 
(CR) . 

8.  Stop  monochromator  drive. 
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Figure  31.  2  to  3  pm  Section  of  SWIR  Blackbody.  Data  taken  at  10  m 

full-scale,  3s  time  constant,  3  mm  slits,  at  a  scan  rate 
of  320/90  with  the  blackbody  =  369.5  K.  Data  from 
November  1985 


Figure  32.  MWIR  Responsivity  Curv« 
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Figure  33.  SWIR  Responsivity  Curve 


Checklist  8 


The  Filter  Wheel 


1.  The  filter  wheel  contains  six  positions,  ostensibly,  2,  4,  6,  8,  10,  12. 

2.  Currently  only  positions  2  and  4  contain  the  appropriate  filters.  The 
other  positions  are  either  open  or  blocked. 

3.  A  DVM  readout  aids  in  determining  the  filter  wheel  position  with  an 

analogous  resistance  for  each.  However,  reliable  DVM  readings  are  rarely 
obtained  so  the  filter  wheel  position  must  be  recorded.  >• 

4.  The  wheel  may  be  rotated  with  the  adjacent  power  supply.  With  this 
supply  on  and  the  voltage  set  to  <5  volts,  10  clicks  of  the  green  panel 
button  will  move  the  filter  wheel  to  the  next  position. 
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3.4  Taking  Discharge  Data 


Discharge  data  are  the  essence  of  COCHISE  research.  The  experiments  are 
necessarily  diverse.  However,  there  are  several  routine  steps  which  may  be 
described  here.  Briefly,  an  experiment  is  performed  as  follows:  1)  the 
experimental  matrix  is  always  laid  out  in  advance,  2)  once  the  calibrations 
have  been  completed,  the  gas  lines  are  heat  soaked  for  1  hr  at  120  K  (all 
discharge  heat  sinks,  both  manifolds,  and  the  pressure  port  plate)  then  set  to 
90  K,  3)  the  desired  gas  flows  (both  discharge  and  counterflow)  are  set  up  and 
a  few  preliminary  spectra  are  taken  scanning  the  experimental  variables, 

4)  discharge  data  is  almost  always  taken  AC  although  DC  operation  has  also 
been  performed,  5)  the  experimental  variables  (e.g.,  percent  discharge  N2, 
percent  counterflow  of  added  gas)  are  examined  in  detail  under  varying 
resolution  conditions,  6)  with  time  permitting  interesting  phenomena  from  tha 
current  or  previous  runs  are  also  examined,  and  7)  when  the  mass  loading  on 
the  chamber  becomes  too  great  (P  >  2  x  10“6),  it  is  time  for  warmup. 

The  routine  aspects  of  taking  discharge  data  which  may  be  useful  to  the 
novice  include:  1)  setting  up  the  gas  flows,  2)  operating  the  discharges, 

3)  some  sample  spectra  from  previous  COCHISE  experiments,  and  4)  AC  versus  DC 
operation.  Each  of  these  will  now  be  discussed. 

The  gas  flows,  both  discharge  and  counterflow,  are  typically  set  up  such 
that  each  mass  flow  rate  is  0.11  g/s.  The  flow  calibrations  in  Figures  5 
through  14  will  aid  in  setting  the  flows  of  each  component.  The  initial  flows 
are  generally  set  up  with  the  reaction  cell  valves  off  and  the  pump-out  valves 
open.  This  minimizes  gas  loading  to  the  chamber.  With  the  flows  established 
aud  flawing  into  the  reaction  chamber,  it  is  time  to  turn  on  the  discharges. 

The  power  supplies  for  all  four  microwave  discharges  are  located  in  a 
movable  panel.  Once  the  warming  cycle  has  been  completed,  the  power  may  be 
turned  on  and  the  cus^ barges  initiated  with  a  Tesla  coil  discharge  to  an 
exposed  wire  located  beneath  the  shelf  on  the  COCHISE  chamber.  Each  of  the 
discharges  has  a  forward  and  reflected  power  meter  so  operation  may  be 


verified  here  as  well  as  by  displaying  the  output  on  a  scope.  After  each  run 
the  discharges  are  turned  off  to  reduce  the  heat  loading  to  the  chamber. 

Each  COCHISE  run  is  remarkably  different  from  all  previous  experiments. 
There  are,  however,  several  common  discharge  phenomena  which  are  frequently 
encountered  and  examined  in  COCHISE.  Sample  spectra  of  these  phenomena  are 
shown  as  a  reference  guide  to  the  experienced  and  novice  COCHISE  user. 

Figures  34  and  35  show  the  NO  fundamentals  and  overtones  in  the  MWIR  and  SWIR, 
respectively.  These  data  were  obtained  by  reacting  Ar/02  in  the  counterflow 
with  discharged  Ar/N2*  The  CO  fundamentals  in  the  MWIR  are  shown  in 
Figure  36.  This  spectrum  was  generated  by  reacting  a  counterflow  of  Ar/CO 
with  discharged  Ar/N2*  A  flow  of  Ar/N2  through  the  discharges  reacting  with 
Ar  counterflow  will  reveal  the  N2*  IR  transitions.  Figure  37  shows  the 
N2(w1A  *  a1!!)  and  +  B3R)  Wu-Benesch  emissions  in  the  SWIR  obtained  in 

this  manner.  One  last  spectrum  will  be  included  here  for  completeness  is  one 
which  is  generated  with  discharged  Ar  only,  producing  Ar(I)  emission.  Figure 
38  shows  the  Ar(I)  emission  observed  in  the  SWIR. 

One  last  point  to  be  covered  in  this  section  is  AC  versus  DC  operation. 
Most  discharge  spectra  are  taken  in  AC  mode  in  order  to  discriminate  against 
emitters  out  of  phase  with  direct  discharge  excitation  but  DC  data  collection 
has  also  been  performed  in  order  to  maximize  signal.  The  differences  between 
AC  and  DC  operation  are  shown  in  Table  6. 

3 . 5  Overnight  and  Warmup  Procedures 

A  large  number  of  systems  have  to  be  shut  off  every  evening  so  a 
checklist  is  provided  here  to  ensure  that  none  are  forgotten.  Checklist  9 
shows  the  overnight  shut-down  procedures.  The  most  critical  of  these  are 
turning  off  all  heaters,  opening  all  chamber  valves,  and  turning  off  the  gas 
flows  and  the  detector.  If  a  compressor  were  to  fail  at  some  point  in  the 
evening  (which  nas  happened)  forgetting  to  turn  off  any  of  the  above  could 
lead  to  a  catastrophe.  This  must  be  kept  in  mind  whenever  the  system  is  left 
unattended  for  any  period  of  time. 
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Figure  34.  NO  Fundamental  Emission 


Figure  35.  NO  Overtone  Emission 


MAXIMUM  INTENSITY  =  59x10  Wcm  sr  Mm 
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WAVELENGTH,  urn 

Imax  *  2.12  x  10-!>  M  cm"2  s r'l  um~l, 

Low  N2,  Ar  Counterflow  On 

A-24B3 

Figure  37a.  Spectrum  of  the  AC  SWIR  Emission  Observed  from  Discharged 
Nitrogen  in  Argon  with  an  Argon  Counterflow.  COCHISE  file 
number  326534.  Resolution  *  0.013  ym.  The  flow  condi¬ 
tions  were  2500  ymole/s  Argon,  20  ymole/s  nitrogen  through 
the  discharge  with  a  flow  of  2700  ymole/s  of  argon  in  the 

counter flow.  ,  , 
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Figure  37b.  Calculated  Best-Fit  Spectrum  to  the  Data  Shown  in  Figure  37a 
Using  N2(w1au  >  a1!^)  Emission  as  a  Basis  Set.  Tr  =  200  K, 
Resolution  **  0.013  ym. 
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Figure  37c.  SWIR  Emission  Spectrum  Observed  from  Discharged  Ar/N2  with  Ar 
Counterflow  Under  Conditions  of  High  N2  Flow.  COCHISE  file 
326531.  Resolution  =  0.013  um.  The  flow  conditi< ns  were 
2500  ymole/s  Ar  and  250  ymole/s  N2  through  the  discharge  and 
2700  umole/s  Ar  through  the  counterflow. 
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Figure  37d.  Calculated  Best-Fit  Spectrum  to  the  Data  Shown  in  Figure  37c 
Using  N2(W3Au  +  B3Ug)  and  N2(w’au  a 1  n g )  Emission  as  Basis 
Sats.  Tr  ( W-*-B )  =  100  K,  Tr(w+a)  =  300  K,  resolution  = 

0 . 0  1  3  u  m . 
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Figure  38.  High  Resolution  Spectrum  of  Ar  Emission  from  2  to  4  pm, 
0.007  pm  Resolution,  0.1  Torr.  The  spectrum  was  taken 
using  a  2  pm  long-pass  interference  filter  and  a  grating 
blazed  at  3  pm,  operating  in  first  order. 


Table  6.  AC  Versus  DC  Data  Collection 


Parame ter 

AC 

DC 

Detector  output 

Channel  A  PAR 

Channel  A  PAR 

PAR  reference 

Discharge  output  through 
light-link  board 

Krohn-Hite  output 

Discharge  control 

Square  wave 

External  drive 

Scope  1  (2-channel 
scope ) 

Signal  output 

Chopper  output,  signal 
output 

Scope  2  (4-channel 
scope ) 

Discharge  output 

Discharge  output 

Chopper 

Off 

On 
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Checklist  9 


Overnight  Shutdown  Procedures 


1.  Gas  flows,  solenoid  valves,  discharges,  scan,  chart  drives  off. 

2.  Disk  up,  computer  locked. 

3.  Bias,  detector  off. 

4.  Console  off:  PAR 

Chopper/preamp 
Recorder 
RC  heat 
Grating  drive 

Tektronix  two-channel  scope 
4-channel  scope 
Flowmeters  and  control  board 
Heater  control  board 

5.  RC,  bulkhead  valves,  open,  SP  valve  always  open,  panel  locked. 

6.  Check  status  of  160  liter  Dewar  (relevant  only  if  LN2  cooling  of  gae 
line  is  used) 

7.  LN2  to  77  K  plate  off,  exit  open  (relevant  only  if  LN2  cooling  of  gas 
line  is  used). 

8.  Microwave  units  off. 

9.  Gas  cylinders  off,  toggle  valves  closed. 

10.  Side  shelf:  Unplug  spectrometer  heater 

Filter  ohm  meter  voltage  source  off 
LHe  level  meter  off 

11.  Refrigerator  to  overnight  r-ondition:  two  compressors,  low  speed 

LN2  valve  three  turns  open 

12.  Message  on  teletype  (.Goodnight). 
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With  the  experiment  completed,  it  is  time  to  warmup  the  chamber.  Great 
care  must  be  exercised  in  implementing  this  procedure  as  severe  explosions 
causing  extensive  chamber  damage  have  occurred  during  warmup  in  the  past. 
During  warmup  the  critical  procedures  are  turning  off  all  heaters  (except  the 
spectrometer  heater  which  must  be  on)  leaving  all  chamber  valves  open,  and 
turning  off  the  diffusion  pump.  Checklists  10  and  11  show  the  warmup  pro¬ 
cedure  in  great  detail. 

3.6  Routine  Maintenance 


A  section  such  as  this  can  by  no  means  be  complete  as  further  experience 
will  continually  add  more  items  to  this  list.  What  follows  is  a  compilation 
of  issues,  which  from  current  experience,  are  deemed  important  to  the 
long-term  smooth  operation  of  COCHISE. 

1.  Grease  all  refrigerator  fittings  (with  a  lithium  grease)  at  least 
once  a  year.  Do  not  overgrease. 

2.  Replace  refrigerator  activated  charcoal  traps  after  3000  hr  of 
operation. 

3.  Always  keep  a  positive  helium  pressure  on  the  detector  reservoir. 

This  will  keep  moisture  out. 

4.  Check  bias  (detector)  batteries  before  each  run  (should  be  >25  volt). 
Change  after  every  second  run. 

5.  Check  and  replace  spent  gas  tanks.  Order  more  if  needed. 

6.  Ordering  liquid  refrigerants: 

LN2  requires  a  2-day  lead,  deliveries  are  needed  twice  a  week 

LHe  requires  a  4  to  5-day  lead,  one  500  liter  Dewer  is  needed  for 
a  complete  run. 
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Checklist  10 


Warmup  Procedures 


1  . 
2. 


3. 


4. 

5. 

6. 


7. 


8. 


9. 
10. 
1 1  . 


12. 


13. 

14. 


Go  through  overnight  shutdown. 

Put  Baratron  on  1000  Torr  head  (switch  DC  output  cable  in  back  to 
1000  Torr  head. 

Attach  Baratron  leads  to  chart  recorder: 

Disconnect  D/A  converter  leads 

Use  both  pens 

Put  on  slowest  speed 

Use  two  different  scales,  100  and  10  Torr 

Set  plumbing  for  Baratron  to  record  Discharge  #1  pressure. 

Raise  disks. 

Leave  MOTHER  operating. 

Open  spectrometer  slits  to  3.0  mm 

Plug  in  spectrometer  heater  and  set  to  40  to  50  volts. 

Monitor  in  MOTHER 

Shroud  -  warm  end,  RC,  discharge  manifold,  RF,  shroud  -  cold  end 
37  36  35  33  16 


Spectrome ter 
12 


Turn  off  diffusion  pump  heaters  on  valve  board. 

Turn  off  LN2  supply  to  refrigerator  (valve  near  window). 
Ensure  spectrometer  is  heating. 


Turn  off  the  second  of  two  compressors  operating  3 
(leave  compressor  #1  on). 


Turn  on  recovery  switch. 


Refrigerator 


Turn  off  last  compressor. 


On  side  shelf,  turn  off  vacuum  space  vacuum  gauge. 
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Checklist  1 1 


Continued  Warmup  Procedure 


1.  All  temperatures  >  220  K  (2  to  3  days). 

2.  Turn  off  spectrometer  heaters. 

3.  Close  right  angle  pump  valve  -  on  board 

Close  foreline  valve 
Turn  off  roughing  pump 

4.  Ensure  that  all  chamber  valves  are  open. 

5.  Use  LN2  offgas  fill  line  to  pressurize  chamber  to  1/2  atmosphere  while 
monitoring  dial  gauges  for  excessive  pressure  differentials. 

6.  Monitor  1000  Torr  pressure  gauge  to  ensure  no  large  pressure 
differentials  are  created. 

7.  When  the  temperatures  are  ~290  K,  use  LN2  offgas  fill  line  to  bring 
chamber  to  1  atmosphere. 

8.  Release  door  lock  so  it  will  crack  when  1  atmosphere  is  reached. 

9.  The  pressure  may  be  read  on  three  pressure  meters  located  beneath  the 
chamber  door. 

10.  The  chamber  may  now  be  opened. 

11.  Never  open  chamber  below  290  K  or  severe  condensation  will  result. 
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7. 


Changing  diffusion  pump  oil: 


The  oil  may  be  changed  by  draining  through  the  bottom  viewing  port 
and  filling  through  a  port  in  the  side  arm  after  removing  the  small 
direct  drive  pump  and  its  Cu  tubing.  Make  sure  the  diffusion  pump  is 
cool  and  vented  (the  whole  chamber  will  have  to  be  vented).  The  pump 
requires  1  gal  of  oil  (we  currently  use  702  silicone  fluid) .  This 
should  be  done  every  2  or  3  years.  The  oil  does  not  seriously 
degrade  but  the  volume  does  decrease.  The  drain  port  seals  with  a 
copper  washer,  be  sure  to  replace  this  when  finished. 

8.  Normal  operation  of  refrigerator: 

High  pressure  meter  should  read  225  psi  (lower  if  cold) 

Return  temperature  -10*C  to  -20°C 
Low  side  pressure  5  to  8  psi 

9.  All  mechanical  pumps  should  have  their  pump  oil  changed  every  year. 
This  includes  the  gas  manifold  pump  located  under  the  floor  boards. 
The  Leybold-Heraeus  direct  drive  forepump  (located  in  the  back  room, 
against  the  far  wall)  should  have  its  oil  changed  after  every  run. 

To  drain,  remove  drain  plug  at  the  base  and  start  the  pump  (briefly) 
two  or  three  times.  Fill  with  direct  drive  pump  oil  in  the  fill  port 
on  top . 
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4.  CONCLUSIONS  AND  PROPOSED  CHANGES 


The  goal  of  this  manual  has  been  to  provide  a  thorough  discussion  of  all 
facets  of  COCHISE  operation.  It  is  by  no  means  exhaustive  but  it  is 
sufficiently  detailed  to  be  a  valuable  learning  tool  for  the  novice  and  a 
useful  guide  for  the  experienced  operator. 

Several  changes  are  currently  planned  for  the  COCHISE  apparatus.  These 
include:  installation  of  a  new  detector  with  greater  sensitivity  and  faster 

response  time,  replacement  of  the  current  grating  drive  with  a  stepper  motor, 
complete  reconstruction  of  the  gas  manifold  to  reduce  system  leaks,  and  the  ' 
replacement  of  the  PDP-8E  computer  system  with  a  Micromint  for  temperature 
control  and  an  HP-PC  for  data  acquisition.  These  changes  will  be  performed  as 
a  cooperative  venture  between  PSI,  AFGL  and  Utah  State  University/Stewart 
Radiance  Lab.  The  involved  PSI  personnel  are  the  authors  of  this  manual,  and 
the  AFGL  and  USU/SRL  scientists  are  Dr.  Steven  Miller  and  Mr.  Dale  Sinclair, 
respectively.  In  view  of  these  proposed  changes,  it  should  be  recognized  that 
several  sections  of  this  manual  will  be  obsolete  within  a  few  years  and 
extensive  rewriting  will  be  required  to  keep  the  manual  up  to  date. 
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ABSTRACT 


Eleven  transitions  in  the  W^Au-*-B^IIg  (W-B)  and  w^Au>a^IIg  (w-a)  systems  of 
nitrogen  have  been  observed  in  the  infrared  including  the  previously 
unobserved  (1,0)  and  (2,1)  W-B  features  at  6.5  and  7.65  pm,  respectively.  The 
fluorescence  spectra  were  observed  in  a  cryogenic  reaction  chamber  at  pres¬ 
sures  of  ~  3  mtorr  (0.4  Pa),  following  expansion  of  flowing  N2/Ar  mixtures 
excited  by  microwave  discharges  at  ~  1  torr.  Einstein  coefficients  for  the 
w-a  system,  calculated  using  a  published  transition  moment  function,  predict 
the  radiative  lifetimes  of  the  lower  vibrational  levels  of  the  w^y  state  to 
be  a  factor  of  three  longer  than  earlier  estimates.  Using  a  spectral  simula¬ 
tion  and  linear  least-squares  fitting  technique,  the  published  W-B  and  calcu¬ 
lated  w-a  branching  ratios  are  verified  for  the  transitions  observed  across 
the  2  to  4  pm  region.  The  observed  vibrational/electronic  state  distributions 
are  not  characteristic  of  those  expected  for  direct  excitation,  but  appear  to 
result  from  extensive  collisional  coupling  among  excited  states  of  nitrogen 
which  occurs  in  the  high  pressure  region  prior  to  expansion. 


* 


f 
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Introduction 


There  are  several  nitrogen  electronic  transitions  which  give  rise  to 
emission  in  the  2  to  8  um  infrared  region.  1-0  These  systems  include 
W^Au~B^ng ,  w^Au+a^IIg,  and  a^g+a'^Ey-.  Emissions  from  the  B'  -3Zu--*-B3rig  and 
B-^IIg->-A^Eu+  systems  also  extend  into  this  spectral  region.  These  systems  are 
commonly  observed  in  nitrogen  discharge  plasmas  and  several  have  been  indi¬ 
cated  as  potential  infrared  auroral  transitions.10  The  detailed  excited  state 
populations  and  vibrational  distributions  responsible  for  the  observed  fluo¬ 
rescence  provide  information  on  mechanisms  of  excitation  and  energy  disposal 
in  excited  nitrogen.  Interore tation  of  auroral  or  plasma  spectra  for  these 
band  systems,  however,  requires  prior  verification  of  the  transition  branching 
ratios.  We  report  here  observation  of  emission  from  several  vibrational 
levels  in  w°Au-*-B3ng  (v'  =  1-5)  and  in  w^u+a^Hg  (v*  =  0-2),  and  verification 
of  their  radiative  branching  ratios  across  the  2  to  4  pm  region. 

Experiments 

These  experiments  were  performed  in  the  COCHISE  (COld  CHemiexcitation 
Infrared  Stimulation  Experiment)  cryogenic  discharge  afterglow  apparatus, 
which  is  described  in  detail  elsewhere.11  Excitation  of  nitrogen  electronic 
states  is  produced  in  four  parallel  microwave  discharges  (2450  MHz,  50  w)  of 
flowing  ^/Ar  mixtures  at  ~  1  torr  total  pressure.  A  diagram  of  the  reaction 
chamber  is  shown  in  Figure  1.  After  exiting  the  discharge  tubes  the  gas 
expands  into  a  low  pressure  (~  3  mtorr),  cryogenically  pumped  chamber  (~  20  K) 
where  the  molecules  enter  the  collimated  field  of  view  of  a  scanning 
monochromator/infrared  detector  assembly.  Residence  times  in  the  discharge 
tubes  are  on  the  order  of  3  to  5  ms;  an  average  time  of  flight  of 


500  ±  100  ps  is  required  for  the  gases  to  exit  the  discharge  tubes  and  enter 


the  field  of  view.  Opposing  flows  of  argon  are  used  to  create  a  quasi-static 


interaction  region  along  the  centerline  of  the  field  of  view,  resulting  in 
partial  rethermalization  of  the  expansion-cooled  rotational  distributions. 

The  gas  residence  time  in  the  field  of  view  is  ~  0.3  ms.  Gaseous  helium 
refrigerant  maintains  all  internal  temperatures  at  20  K,  excepting  the  gas 
lines  and  optics  which  are  held  at  80  and  40  K,  respectively. 

The  infrared  emissions  are  observed  by  a  cryogenic  0.5  m  Czerny-Turner 
monochromator  equipped  with  a  liquid-helium-cooled  arsenic  doped  silicon 
detector  and  a  grating  blazed  at  3  ym.  'T'he  discharges  are  chopped  with  a 
23  Hz,  50  percent  duty  cycle  square  wave,  and  data  collection  is  performec 
with  a  computer-interfaced  lock-in  amplifier.  Data  over  the  2  to  4  ym  region 
were  taken  with  the  phase  shift  between  reference  and  signal  waveforms  locked 
in  at  the  3.3  or  3.6  ym  features.  The  phase  shift  for  the  4  to  8  ym  spectra 
was  locked  at  the  emissions  at  6.5  or  7.5  ym. 

The  absolute  uncertainty  in  the  wavelengths  (due  to  monochromator  drive 
error)  of  the  spectral  scans  is  ±  0.003  ym.  The  data  were  corrected  for  instru¬ 
ment  responsivity  using  blackbody  calibration  spectra  taken  in  the  300  to  370  K 
range.  The  uncertainty  in  the  accuracy  of  the  blackbody  temperature  is  t  3 
which  results  in  a  relative  uncertainty  of  ±  20  percent  in  the  measured 
2. 5/5.0  ym  intensity  ratios.  This  method  of  calibration  has  been  verified  in 
studies  of  the  fundamental/overtone  ratios  for  NO(v)  vibraluminescence . 1 2 

Spectra  were  taken  for  ^/Ar  mixtures  with  N2  mole  fractions  of  0.006  to 
0.12,  with  a  mass-balanced  counterflow  of  argon  added  to  thermalize  the  observed 
emissions.  Spectra  were  also  taken  without  nitrogen  to  identify  the  background 
Arl  Rydberg  emission.'-*  The  data  were  taken  typically  at  a  resolution  of 
0.013  ym  in  the  2  to  4  ym  region  and  0.040  ym  for  the  4  to  8  nn  region. 
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Spectral  Analysis 


The  interpretation  of  the  2  to  4  pm  discharged  nitrogen  data  is  compli¬ 
cated  due  to  the  coincidence  of  several  of  the  W0Au+B°Hg  (W-B)  and  w^Au+a^rig 
(w-a)  features.  The  transitions  that  are  not  distinguishable  from  each  other 
at  a  resolution  of  0.013  pim  include  the  W-E  (4,1)  and  w-a  (2,1)  at  2.4  ym,  the 
W-B  (3,1)  and  w-a  (0,0)  at  3.6  ym(  and  the  W-B  (4,2)  and  w-a  (2,2)  at  4.0  pm. 
These  overlapping  bands  cannot  be  identified  or  quantified  by  inspection, 
therefore  a  spectral  fitting  technique  is  required. 

The  spectral  analysis  technique  used  for  these  studies  has  also  been 
successfully  applied  to  the  analysis  of  NO(A2£),14  electronic  emission,’5 
and  IF(  b3iJq+)  .  1 6  The  synthetic  spectra  are  generated  by  a  computer  code  based 
upon  the  work,  of  Kovacs'7  and  implemented  through  major  modification  of  a 
program  written  by  Whiting,  et  al.1®*19  The  program  considers  singlet- 
singlet,  doublet-doublet,  and  triple t- triple t  dipole  transitions  for  diatomic 
molecules  in  the  optically  thin  limit.  Line-by-line  transition  frequencies 
are  computed  from  energy  eigenvalues  determined  by  exact  solution  of  the 
Schrodinger  equation,  where  the  matrix  elements  of  the  upper  and  lower  state 
Hamiltonians  are  specified  from  tabulated  spectroscopic  constants.  Following 
well-established  procedures,17*20  the  eigenfunctions  of  the  states  are 
described  as  linear  combinations  of  eigenfunctions  of  hypothetical  pure  Hand's 
case  (a)  eigenstates.  The  transition  amplitudes  are  given  by  the  trans¬ 
formation  of  the  dipole  moment  function  between  the  upper  and  lower 
states.17*20  The  computed  infinite  resolution  emission  spectrum  is  con¬ 
volved  with  the  instrument  scan  function  (in  this  case,  a  symmetric  triangle 
with  full  width  at  half  maximum  as  the  spectral  resolution)  to  create  basis 
sets  for  each  vibrational  level  which  are  then  fit  to  the  experimental 
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spectrum  using  a  linear  least  squares  method.  The  fitting  procedure  yields  a 
determination  of  the  product  of  the  upper  state  number  density  and  the  spon¬ 
taneous  emission  coefficient  of  the  transition.  Rotational  distributions  are 
treated  by  simple  Boltzmann  expressions,  so  that  band-integrated  vibrational 
state  number  densities  and  transition  probabilities  may  be  used. 

Lambda  doubling,  which  is  on  the  order  of  1  cm-1  for  the  W3Au-*B3IIg 
system, 21  is  not  treated  in  this  analysis.  The  W3Au-»-B3ng  emission  system 
totals  27  branches  without  lambda  doubling.  Furthermore,  the  coupling  for  the 
W  state  is  close  to  Hund's  case  (b)  while  the  B  state  coupling  exhibits  tran¬ 
sition  from  (a)  to  (b)  at  higher  rotational  levels.  Despite  these  complex¬ 
ities,  comparison  of  the  line  positions  of  the  synthesized  W3Au-*-B3IIg  (2,0) 
feature  to  published  high  resolution  data21  indicates  the  positions  of  prin¬ 
cipal  lines  are  reproduced  to  within  1  cm"1.  Similarly,  the  agreement  between 
the  synthesized  and  high  resolution22  w^Au+a^IIg  (0,0)  transition  is  excellent, 

<;  0.1  cm"1.  Since  most  of  the  data  presented  here  have  been  taken  and  fit  to  a 
resolution  of  0.013  yin  (12  cm"1  at  3.3  ym)  ,  deficiencies  in  the  spectral 
synthesis  program  on  the  order  of  1  cm"1  or  less  are  too  small  to  be  observed. 
The  spectroscopic  data  sources  for  the  considered  N2  states  were  Huber  and 
Herzberg23  and  Roux22  et  al.  for  w^y  and  a1]^,  Cemy24  et  al.  for  w3au  and 
Effantin25  et  al.  for  B3ng. 

Each  vibrational  transition  is  treated  as  an  independent  basis  function 
in  a  linear  least-squares  fit  to  the  total  spectrum  to  determine  the  band- 
integrated  intensities,  Nvi  Ayiyii.  To  determine  the  excited  state  number 
densities,  accurate  values  for  Ayi  vi.  are  needed.  The  transition  probabilities 
used  for  the  W-B  system  are  the  ab  initio  values  calculated  by  Werner  et  al.2& 
and  shown  in  Table  1.  These  values  are  more  than  a  factor  of  two  smaller 
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than  the  empirical  estimates  of  Covey,  Saum,  and  Benesch,27  but  exhibit  essen¬ 
tially  the  same  branching  ratios;  we  prefer  the  results  of  Werner  et  al.28 
owing  to  their  use  of  a  more  accurate  transition  moment  function.  There  are 
no  published  transition  probabilities  for  the  w^Au->-a^IIg  system.  Consequently, 
we  calculated  them  using  the  Franck-Condon  and  r-centroid  data  tabulated  in 
Lofthus  and  Krupenie28  and  the  transition  moment  function  determined  by  Yeager 
and  McKoy.29  The  results  are  shown  in  Table  2.  The  radiative  lifetimes  for 
the  first  five  vibrational  levels  are  predicted  to  be  1500,  830,  470,  430,  and 
360  us,  respectively,  at  variance  with  earlier  estimates28  of  500  to  100  ys. 

These  radiative  lifetimes  are  more  than  an  order  of  magnitude  longer  than 
those  used  by  Cartwright10  for  auroral  modeling  predictions. 

Using  the  spectral  synthesis  linear  least-squares  fitting  technique 
described  above,  eleven  bands  from  the  N2(  W3Au->-B3ng)  and  N2(w1Au>a1  ng)  systems 
have  been  identified.  Inclusion  of  the  B3Hg*A3Eu+ ,  B1  3Eu“->-B3ng ,  a1  ng->-a  '  1  £u- , 
and  B3TIg-»-w3Au  systems  of  nitrogen  in  the  fits  indicated  no  evidence  of 
emission  from  these  systems  in  the  observed  spectra.  This  is  consistent  with 
the  rapid  radiative  decay  expected  for  the  N2(B,B' ,a)  states  in  transit  from 
the  discharge  exits  to  the  field  of  view. 

» 

The  rotational  temperatures,  Tr,  used  in  the  fits  were  determined 
empirically  by  comparison  of  intensity  distributions  and  branch  structures  of 
the  data  with  theoretical  spectra.  The  (2,0)  band  was  used  for  the  W-B 
system  with  the  result  Tr  =  100  K,  indicating  this  emission  system  is  nearly 
rotationally  thermalized.  With  the  rotational  temperatures  of  the  W-B  system 
fixed  at  100  K,  simultaneous  W-B  and  w-a  fits  were  used  to  determine  the 
rotational  temperature  for  the  w-a  system.  The  optimal  Tr  was  chosen  as  that 
which  best  reproduced  the  intensity  distribution  of  the  first  two  branches  of 
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the  3.6  ym  emission  feature  which  corresponds,  in  large  part,  to  the  R  and  Q 


branches  of  the  w-a  (0,0)  emission.  The  result  for  w-a  was  Tr  =  300  K, 
although  200  K  produced  better  fits  in  some  data  taken  at  low  nitrogen  mole 
fraction.  The  deviation  in  the  calculated  populations  incurred  using  w-a 
basis  sets  with  Tr  =  200  K  versus  Tr  =  300  K  is  <10  percent  in  the  w-a  popula¬ 
tions  and  <5  percent  in  the  W-B  populations. 

r 

Results 

Identification  of  the  observed  spectral  features  was  facilitated  by  a 
variation  of  the  relative  W-B  and  w-a  intensities  as  a  function  of  nitrogen 
mole  fraction.  At  low  nitrogen  mole  fraction  the  w-a  features  are  more 
prominent  allowing  unambiguous  identification.  Under  these  conditions  the 
dominant  spectral  feature  at  3.6  ym  is  the  (0,0)  w-a  transition.  Figure  2 
shows  the  data  and  best  fit  to  the  3.5  to  3.9  ye  region  of  data  obtained  at  a 
nitrogen  mole  fraction  of  1.2  x  10~2.  The  figure  also  shows  the  contributions 
to  the  fit  from  each  of  the  w-a  and  W-B  basis  sets.  The  best  fit  was  obtained 
for  a  w-a  basis  set  with  Tr  =  200  K.  At  higher  nitrogen  mole  fractions  the 
W-B  features  predominate,  in  particular  the  (2,0)  feature  at  3.3  ym.  Figure  3 
shows  data  taken  at  high  nitrogen  mole  fraction  and  fit  to  both  W-B  and  w-a 

*  systems  over  the  full  2  to  4  ym  region.  Six  transitions  in  the  W-B  system  are 

identified,  three  from  the  Av  =  3  series  and  three  from  the  Av  =  2  series. 

A 

The  W-B  (5,3)  band  at  4.5  ym  cannot  be  distinguished  from  the  noise  level,  in 
keeping  with  its  small  transition  probability  relative  to  that  for  the  (5,2) 
band  (cf.  Table  1).  No  W-B  emission  from  vibrational  levels  higher  than 
v'  =  5  is  observed  under  any  nitrogen  mole  fraction  condition.  Five  transi¬ 
tions  in  the  w-a  system  are  also  identified,  the  (1,0)  and  (2,1)  bands  in  the 
Av  =  1  sequence  and  the  (0,0),  (1,1),  and  (2,2)  bands  in  the  Av  =  0  sequence. 
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No  clearly  identifiable  emission  from  v'  =  3  in  the  w-a  system  is  observed. 


All  additional  features  in  Figure  3  are  due  to  Arl  Rydberg  emission,  which 
appears  as  scattered  light  from  the  discharges.  We  have  reported  on  these 

emissions  previously.^ 

Spectra  taken  of  discharged  nitrogen  in  the  6  to  8  ym  region  also  show 
evidence  of  W-B  features.  Fits  to  the  data  provide  unambiguous  identification 
of  the  previously  unobserved  (1,0)  and  (2,1)  transitions  at  6.2  to  6.8  urn  and 
7.4  to  8  pm,  respectively.  Figures  4a  and  4b  show  data  and  fit  to  the  W-B 
(1,0)  and  (2,1)  W-B  features  respectively.  The  W-B  (2,1)  feature  is  the  only 
emission  observed  with  this  apparatus  in  the  7  to  8  ym  wavelength  range. 

Figure  5  shows  the  determined  populations  in  the  v'  =  3  level  of  the  w3&u 
state  and  v'  =  0  of  the  w^u  state  as  a  function  of  nitrogen  mole  fraction. 

The  w-B  and  w-a  emission  systems  coexist  in  the  effluent  of  a  microwave  dis¬ 
charge  under  all  nitrogen  mole  fraction  conditions.  Both  systems  exhibit 
similar  kinetic  behavior  except  at  N2  mole  fractions  in  excess  of  10  percent. 
Under  these  conditions  the  w^u  (V  =  0)  population  drops  off  more  rapidly 
than  the  W^Ay  (V  *  3).  The  W-B  (2,0)  feature  exhibits  dissimilar  and 
variable  temporal  behavior  as  measured  by  phasing  relative  to  the  other  W-B, 
w-a  emission  features  as  a  function  of  nitrogen  mole  fraction.  The  W-B  (2,0) 
feature  is  in  phase  only  at  the  higher  nitrogen  mole  fractions.  At  the  lower 
N2  mole  fractions,  the  phase  lag  of  the  (2,0)  band  is  on  the  order  of  u  (22  ms), 
corresponding  to  cessation  of  the  discharge  pulse.  Since  the  radiative  life¬ 
time  of  the  W^au  (v*  »  2)  level  is  short  (1200  ys)  compared  to  the  modulation 
frequency  (23  Hz),  and  is  comparable  to  the  radiative  lifetimes  of  other 
observed  states  which  exhibit  no  phase  lags,  radiative  decay  effects  have  been 
ruled  out  as  the  cause.  Investigations  of  the  (1,0)  and  (2,1)  W-B  bands  also 
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indicate  variable  phasing  between  these  emissions.  These  two  emissions  are 
in  phase  with  each  other  only  at  low  nitrogen  mole  fractions.  No  phase  shift 
effects  are  encountered  vith  the  w-a  emissions  or  with  the  W-B  (v*  >  3)  bands. 

The  relative  population  distributions  for  the  vibrational  levels  of  the 
W3AU  and  w^Au  states  are  shown  in  Figures  6a  and  6b,  respectively.  Also  shown 
in  Figures  6a  and  6b  are  the  distributions  back-corrected  for  the  500  ys  radi¬ 
ative  decay  which  is  incurred  during  transit  from  the  end  of  the  discharge 
tubes  to  the  field  of  view.  In  the  absence  of  collisional  feed  sources  in  the 
reaction  cell,  the  corrected  distributions  are  those  which  would  exist  at 
the  discharge  exits.  Both  W3AU  and  w^Au  states  exhibit  strongly  relaxed  dis¬ 
tributions,  with  the  maximum  population  occurring  at  the  lowest  vibrational 
levels.  Also  shown  in  Figures  6a  and  6b  are  the  determined  upper  limit  rela¬ 
tive  populations  for  W3au  (v*  =  6)  and  w*Au  (v*  -  3).  The  value  for  w3au 
( v '  =  6)  corresponds  to  the  noise  level  of  the  data.  The  w-a  (3,2)  band 
corresponds  to  a  distinguishable  feature  which  lies  at  2.5  ym,  but  the 
observed  intensity  appears  to  be  principally  attributable  to  a  coincident  Arl 
Rydberg  line . 1 3  (The  w-a  (3,3)  band  at  4.2  yra  is  not  observable  owing  to  its 
very  small  relative  transition  probability  (cf.  Table  2)). 

Emission  from  v'  >  1  levels  in  the  ^(a^g+a1  ^y”)  system  also  occurs  in 
the  2  to  4  ym  region.  Owing  to  the  short  radiative  lifetime  of  the  a1 Eg  state 
(~100  ys) ,  observation  of  this  band  system  in  the  field  of  view  due  to  dis¬ 
charge  excitation  is  unlikely.  Fits  to  the  data  of  the  a-a'  (v1  =  1-3) 
features  indicate  none  of  these  features  are  distinguishable  from  the  noise 
level  in  the  data,  with  determined  upper  limit  populations  of  <2  x  106  cm-3. 
This  is  consistent  with  the  anticipated  production  of  N2(a,  v'  >  1)  in  the 
field  of  view  by  the  observed  w-a  radiative  cascade.  Similarly,  the  absence 
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of  identifiable  emission  in  the  B' Eu“  +  B^IIg  system  gives  an  upper  bound  of 
~  lO^  cm-3  for  the  ^(B1)  number  density  in  the  field  of  view. 

Discussion 

In  comparison  of  the  data  and  the  predictions  of  the  spectral  synthesis 
code,  we  observe  only  minor  disagreements  in  spectral  position  and  branch 
intensity  distribution.  Discrepancies  in  the  reproduction  of  the  branch  head 
positions  are  within  the  spectrometer  drive  error  of  0.003  pm.  Failure  to 
exactly  match  the  relative  branch  intensities  within  a  given  band  may  be 
partially  due  to  non-Boltzmann  rotational  temperature  distributions.  Assump¬ 
tion  of  a  single  Boltzmann  distribution  may  also  cause  discrepancies  in  the 
detailed  reproduction  of  the  relative  branch  intensities.  The  magnitude  of 
these  discrepancies  is  small,  and  at  low  resolution  (6X>0.013  ym)  the  band 
shapes  and  integrated  intensities  are  relatively  well  reproduced.  Because  the 
w-a  and  W-B  systems  can  be  discriminated  at  different  nitrogen  mole  fractions, 
the  fits  confirm  the  (3,0)/(3,1),  ( 4 , 1 )/ ( 4 , 2 )  W-B  and  (1,U)/(1,1),  ( 2 , 1 )/( 2 , 2 ) 
w-a  branching  ratios,  obtained  from  Tables  1  and  2,  to  ±  20  percent. 

The  observed  vibrational  and  electronic  state  distributions  are  quite 
different  from  those  one  might  expect  for  discharge-flow  excitation  under  the 
conditions  of  these  experiments.  Typically,  at  low  N2  mole  fractions,  the 
residence  time  in  the  active  discharge  plasma  (~  2800  cm/s  at  1.4  torr  in  each 
tube)  is  3  to  5  ms,  with  the  discharge  region  extending  essentially  to  the  end 
of  the  tube.  For  N2  mole  fractions  near  0.1,  the  discharge  region  "shrinks" 
to  about  2  ms,  allowing  1  to  2  ms  of  flow  time  downstream  of  the  dosed  volume 
prior  to  expansion  into  the  observation  volume.  For  typical  conditions  of 
this  type  of  discharge  ( [ e~ ]  ~  1011  cm-3,  E/N  ~  10~16V  cm2,  characteristic 
electron  energy  ~  6  to  8  eV) ,  we  expect  the  principal  excitation  mechanism  for 
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N2(W,W)  to  be  via  direct  electron  impact  on  ground  state  nitrogen.  This  pro¬ 
duction  term  would  be  balanced  by  losses  due  to  electron  impact  processes 
(dissociation,  superelastic  de-excitation),  collisional  deactivation,  radi¬ 
ative  loss,  and  flow  out  of  the  discharge.30 

For  the  discharge  conditions  used  in  these  studies,  the  dominant  loss 
terms  are  dissociation,  surface  quenching,  and  flow,  which  should  not  affect 

y 

the  vibrational/  electronic  state  distribution  formed  in  the  excitation 
^  process.  Since  electron  impact  excitation  gives  rise  to  Franck-Condon  vibra¬ 

tional  distributions  relative  to  ground  state  nitrogen,  we  would  expect  to 
observe  significant  fluorescence  from  higher  vibrational  levels  (v1  =  4-9)  of 
the  W  and  w  states,  with  a  preponderance  of  the  triplet  over  the  singlet, 
qualitatively  similar  to  the  state  distributions  predicted  by  Cartwright10  for 
an  aurora.  In  contrast,  however,  we  observe  relatively  "cold"  vibrational 
distributions  with  roughly  equal  populations  in  the  triplet  and  singlet 
states,  despite  their  1.5  eV  energy  difference.  These  observations  suggest 
that  there  are  rapid  collisional  "scrambling"  reactions  occurring  in  the  dis¬ 
charge  which  can  significantly  alter  the  excited  state  population  distribu¬ 
tions.  A  possible  example  of  such  a  reaction  would  be  the  enhancement  of  the 
J  lower  levels  of  the  W  state  by  fast  reaction  of  nitrogen  atoms  with  singlet 

metastables : 


N(4S)  +  N2(a'1Zu-)  +  N2(W3Au,  v  4  5)  +  N(4S) 

Similar  rapid  intersystem  crossings  may  occur  through  collisions  of  excited 
N2*  with  N2  and  Ar.  Indeed,  previous  observations  of  fluorescence  from 
N2;s3ng)  excited  by  an  electron  beam  gave  evidence  of  strong  collisional 
coupling  between  the  neighboring  vibrational  levels  of  the  B  and  W  states.31 
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The  excitation  mechanism  for  N2(w)  is  less  apparent  owing  to  its  high 
energy  relative  to  the  neighboring  states  of  N2* .  The  non-thermal  rotational 
distributions  we  observe  for  the  w+a  transitions  suggest  formation  of  N2(w)  by 
a  reactive  collision,  although  direct  energy  transfer  to  N2  cannot  be  ruled 
out.  In  the  latter  case,  however,  the  precursor  states  would  have  to  be 
medium  to  high  vibrational  levels  of  one  or  more  of  the  neighboring  N2*  elec-  „ 

tronic  states.  An  intriguing  possibility  is  excitation  of  N2(w)  via 
metastable-metastable  interactions,  e.g.,  N2*  +  N2*  or  N*  +  N2* ;  precursor  4 

number  densities  on  the  order  of  10^  cm~^  (in  the  active  discharge  at 
1017  cm“3  total  number  density)  and  a  near  gas-kinetic  rate  coefficient  would 
be  sufficient  to  account  for  the  quantity  of  N2(w)  observed  in  the  field  of 
view. 

The  large  phase  lag  between  N2(W,  v'  *  1,2)  and  the  remainder  of  the 
excited  N2  states  is  difficult  to  understand.  Evidently  these  states  are  not 
formed  as  rapidly  (and/or  are  destroyed  more  rapidly)  in  the  active  discharge 
as  W(v'  *  3-5)  and  w(v'  *  0-2),  but  are  rapidly  formed  by  collisional  energy 
transfer  following  cessation  of  the  discharge.  This  could  occur  through  a 
relatively  inefficient  process  such  as  energy  transfer  with  N2,  such  that  at 
lower  N2  mole  fractions  a  substantial  amount  of  N2(W,  v**1,2)  cannot  be  formed 
until  after  the  precursor  and  N?(W)  loss  processes  in  the  active  discharge  are 
terminated.  The  precursor  may  then  be  a  sufficiently  energetic  afterglow 
species  such  as  N2(A,  v>9)  or  N2(X,  v>33). 

The  collissional  processes  for  energy  redistribution  among  the  states  of 
n2*  implied  by  these  data  must  occur  with  relatively  high  efficiency,  with 
collision  efficiencies  ranging  from  near  unity  for  a tom-meta stable  and 
metastable-metastable  reactions  to  perhaps  1  percent  for  collisions  with 
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N2 ( X ,  v=0).  Such  processes  will  clearly  affect  the  energy  flow  in  N2  plasma 
systems  at  pressures  of  a  few  torr  or  greater.  Furthermore,  some  of  the  col- 
lisional  energy  shuffling  reactions  may  be  fast  enough  to  affect  N2*  state 
distributions  in  the  auroral  upper  atmosphere,  particularly  in  the  80  to 
100  km  altitude  regime  where  collisional  quenching  competes  with  radiative 
deactivation  of  excited  species.  This  possibility  could  be  tested  by  observa¬ 
tions  of  altitude-dependent  N2(W,w)  fluorescence  in  a  strong  aurora  and  com¬ 
parison  of  the  observed  vibrational/electronic  distributions  with  those 
expected  for  direct  auroral  excitation. 

Conclusions 

We  have  observed  infrared  (2  to  8  ym)  fluorescence  from  nitrogen  in  the 
Av  =  0  and  1  progressions  of  the  w^u+^IIg  band  system  and  in  the  Av  =  1,  2, 
and  3  progressions  of  the  W^Au+B^IIg  band  system  resulting  from  microwave  dis¬ 
charge  excitation  of  N2/Ar  mixtures.  Using  a  spectral  synthesis  procedure 
with  literature  values  for  the  spectroscopic  constants  provides  adequate 
description  of  the  observed  spectra  for  spectral  resolutions  with  full  width 
at  half  maximum  greater  than  or  equal  to  0.007  ym  (6.4  cm-1  at  3.3  ym) . 
Experimental  observations  of  relative  intensities  of  different  transitions 
arising  from  a  common  upper  state  support  the  theoretical  predictions  for  the 
dipole  moment  functions  of  Werner^G  for  W-B  and  Yeager  and  McKoy^  for  the  w-a 
system.  We  report  w-a  Einstein  coefficients  calculated  from  the  data  given  in 
Ref.  29;  these  values  differ  markedly  from  those  used  by  Cartwright10  in  his 
estimation  of  the  relative  N2(w1Au)  auroral  number  densities.  Recalculation 
with  the  new  w-a  Einstein  coefficients  leads  to  an  increase  in  the  predicted 
high-altitude  N2<w1Au)  relative  number  density  by  an  order  of  magnitude. 
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The  identification  of  the  (1,0)  and  (2,1)  W-B  bands  represents  the  first 


time  these  bands  have  been  experimentally  observed.  The  observed  W-B  and  w-a 
vibrational  distributions  provide  evidence  of  complex  collisional  coupling  of 
the  discharge-excited  states  of  N2,  in  particular  a  collisional  feed  of  W-^AU. 
More  detailed  time-resolved  kinetic  studies  will  be  required  to  resolve  this 
issue . 
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Table 

1.  Einstein  coefficients  for  the 

(in  103  s-1)  as  calculated  by 

W3Au-*-B3ng 

Warner  et 

transitions 

al.26 

v'/v" 

0 

1 

2 

3 

4 

5  6 

0 

0.000 

1 

0.221 

2 

0.735 

0.084 

3 

0.868 

0.764 

0.014 

4 

0.649 

1.525 

0.457 

0.000 

5 

0.377 

1.619 

1.554 

0.183 

0.000 

6 

0.  187 

1.222 

2.311 

1.157 

0.043 

0.000 

1 


A 
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Table  2. 

Computed 

(in  1 03 

Einstein 

s-  1 ) 

coefficients  for 

the  w  ^  A 

u+a 1 Hg  transitions 

> 

\ 

> 

0 

1 

3 

4 

5  6 

0 

0.636 

0.013 

1 

0.979 

0.206 

0.013 

2 

0.660 

1.03 

0.444 

0.010 

3 

0.308 

1.199 

0.767 

0.002 

0.006 

4 

0.117 
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Figure  Captions 

Figure  1.  Diagram  of  the  COCHISE  reaction  chamber.  The  physical  dimensions 
of  the  cell  are  0.6  m  in  length  and  0.4  m  in  diameter. 

Figure  2.  Data  (light  line)  and  best  fit  (dark,  line)  (a)  to  the  (0,0)  and 
1  (1,1)  N2< w1  Au-*-a 1  Hg)  and  (3,1)  N2( w3Au*B3Hg)  emission  systems  at 

rotational  temperatures  of  200  and  100  K,  respectively.  The 

) 

spectrum  was  obtained  at  a  discharge  nitrogen  mole  fraction  of 
1.2  x  10""2  anci  a  resolution  of  0.0067  ym.  The  maximum  intensity  of 
the  spectrum  is  3.45  x  10“  ^  W  cm“2  sr"1  um-1.  The  w-a  and  W-B 
basis  sets  which  comprise  the  best  fit  in  (a)  are  shown  in  (b)  and 
(c),  respectively. 

Figure  3.  Data  (light  line,  shown  with  the  best  fit  and  separately  above)  and 
best  fit  (dark  line)  to  the  indicated  ^(w1  Au+a'  ng)  and 
N2(w3Au+B3iig)  emission  systems  at  rotational  temperatures  of  300 
and  100  K,  respectively,  and  at  a  spectral  resolution  of  0.013  pm. 
The  spectrum  was  obtained  at  a  discharge  nitrogen  mole  fraction  of 
j  9.1  x  10”2.  The  maximum  intensity  of  the  spectrum  is  3.4  x  10“^ 

W  cm-2  sr"'  pm"1.  The  observed  features  not  accounted  by  the  fit 
are  due  to  Arl  Rydberg  emission.12 
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Figure  4.  Data  (light  line)  and  best  fit  (dark,  line)  to  the  N2(W2Au+B2IIg) 

Av  =  1  transitions.  Shown  in  (a)  is  the  (1,0)  band  fit  at  a 
rotational  temperature  of  100  K.  The  spectrum  was  taken  at  a 
nitrogen  mole  fraction  of  7.4  x  10”2  and  has  a  maximum  intensity  of 
2.6  x  10“ 10  W  cm”2  sr“l  ym”1.  Shown  in  (b)  is  the  the  (2,1)  band 
fit  at  a  rotational  temperature  of  120  K.  The  spectrum  was  taken 
at  a  nitrogen  mole  fraction  of  1.2  x  10”^  and  has  a  maximum 
intensity  of  3.6  x  10”^°  W  cm”2  sr”1  ym-1. 

Figure  5.  Determined  populations  for  W2au  (v*  =  3)  (•)  and  w^u  (v*  =  0) 

(o)  versus  discharge  nitrogen  mole  fraction.  Absolute  number 
densities  pertain  to  the  3  mtorr  interaction  region  in  the  center 
of  the  reaction  volume. 

Figure  6.  a)  Relative  vibrational  population  distribution  for  w3au  at  a 

discharge  nitrogen  mole  fraction  of  1.2  x  10”1  (•)  normalized  to 
v'  =  2.  Also  shown  are  the  data  back-corrected  for  a  500  ys 
radiative  decay  (o).  b)  Relative  vibrational  population 
distribution  for  w^Au  at  a  discharge  nitrogen  mole  fraction  of 
1.2  x  10_1  (•)  normalized  to  v*  =  0.  Also  shown  are  the  data 
back-corrected  for  a  500  ys  radiative  decay  (o). 
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